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CHAPTER  1 
INTRODUCTION 


The  printed  circuit  transmission  line  has  become  an  integral  part 
of  many  high  frequency  circuit  applications,  creating  an  ever 
increasing  demand  for  accurate  theoretical  predictions  of  circuit 
performance.  While  approximate  formulations  were  acceptable  in  the 
past,  rising  performance  expectations  as  well  as  the  desire  for 
circuits  that  don't  require  post-production  tuning  have  made  the  use  of 
more  exact  formulations  necessary.  This  report  uses  exact  formulations 
to  address  four  major  areas  of  concern  in  high  frequency  printed 
circuit  design;  complex  microstrip  structures  (multiple  lines  and/or 
dielectric  layers) ,  coupling  between  lines,  distortion  of  non-periodic 
signals  on  complex  structures,  and  a  new  method  to  control  coupling  on 
multilayer  structures,  as  well  as  presenting  numerical  results  for  each 
of  these  areas.  This  chapter  discusses  the  uses  of  microstrips, 
discusses  the  various  types  of  approaches  used  to  solve  the  above 
mentioned  problems,  and  outlines  the  solution  methods  and  numerical 
results  presented  in  this  report. 

1.1  History  and  Applications  of  Microstrips 

While  the  use  of  printed  circuit  transmission  lines  is  a 
relatively  recent  development,  today  its  technology  is  considered  to  be 


mature.  The  first  printed  circuit  transmission  line,  the  strip  line, 
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was  developed  in  World  War  II  while  the  microstrip  was  first  proposed 
and  investigated  in  the  early  fifties.  Yet,  by  1955,  the  IRE 
Transactions  on  Microwave  Theory  and  Techniques  [1]  devoted  an  entire 
issue  to  the  research  and  design  of  this  new  form  of  microwave 
circuitry.  Early  research  in  this  field  was  concerned  with  calculating 
the  low  frequency  parameters  of  the  microstrip  using  a  quasi-static  TEM 
approach.  Demands  for  higher  operating  frequencies,  shorter  pulse 
widths,  and  smaller  circuit  dimensions  required  that  researchers 
address  the  frequency  dependent  nature  of  the  microstrip  parameters. 
Today,  the  research  field  is  concerned  with  obtaining  more  rigorous 
solutions  of  complex  microstrip  structures  that  include: 

1.  Multiple  dielectric  layers. 

2.  Multiple  conductors  on  different  layers. 

3.  Anisotropic  dielectrics. 

4.  Operating  frequencies  that  approach  the  optical  range. 

5.  Discontinuities  in  the  line. 

6.  Losses  in  the  conductors  and  the  dielectrics. 

7.  Losses  due  to  radiation. 

The  design  of  microwave  integrated  circuits  has  become  very  advanced, 
so  that  there  are  several  commercial  computer  aided  design  packages  on 
the  market  [2] .  Thus  microstrip  circuit  design  has  become  sufficiently 
advanced  to  be  considered  a  mature  technology. 

In  spite  of  the  microwave  integrated  circuit's  (MIC's)  low  power 
handling  capabilities,  it  enjoys  increasing  popularity  due  to  its  many 
advantages  over  waveguide  and  coaxial  lines.  MIC's  are  relatively 
inexpensive  to  make  and  the  fabrication  process  is  simple.  Because  of 


3 


its  planar  structure,  modern  surface  mounted  components,  such  as  GaAs 
FET's,  are  easily  included  in  the  circuit  structure,  giving  the  MIC 
circuit  a  low  profile  and  reduced  weight.  MIC' s  are  very  versatile  so 
that  most  of  the  common  passive  devices,  such  as  inductors,  capacitors, 
and  resistors,  can  be  realized  using  microstrips,  eliminating  the  need 
to  add  these  passive  devices  as  discrete  components.  In  addition, 

MIC's  can  be  used  to  make  phase  shifters,  mixers,  amplifiers,  filters, 
power  dividers,  power  combiners,  and  directional  couplers.  As  an 
antenna,  the  MIC  is  often  used  on  the  exterior  skin  of  aircraft  and 
missiles  because  its  low  profile  does  not  affect  the  aerodynamics  of 
the  vehicle. 

MIC's  are  currently  being  used  in  a  wide  range  of  high  frequency 
applications  that  affect  all  fields  of  electrical  engineering.  MIC 
components  have  been  designed  to  operate  in  the  UHF,  microwave,  and 
millimeter  wave  regions,  and  its  use  is  being  extended  into  the  low 
optical  frequencies.  Current  practical  applications  include  high  speed 
digital  transmission  and  optoelectronic  switches.  Microstrips  have 
even  been  applied  to  pulse  shaping  for  a  laser  used  to  start  a  fusion 
reaction  [3] . 

1.2  The  Effective  Dielectric  Constant  of  a  Single  Microstrip 

One  of  the  most  common  figures  of  merit  used  when  describing  a 

microstrip  line  is  the  frequency  dependent  effective  dielectric 

constant,  c  „^(f).  The  frequency  dependence  of  c  is  due  to  the 

reff  reff 

discontinuity  of  the  dielectric  at  the  conductor  interface,  prohibiting 
the  propagation  of  a  pure  TEM  mode.  This  frequency  dependence  is 


4 


responsible  for  the  microstrip's  dispersive  behavior  at  higher 
frequencies,  which,  in  general,  is  an  undesirable  property.  Although 
the  early  efforts  in  research  emphasized  the  use  of  a  TEM  mode  or 
quasi-static  approximation,  later  work  in  the  field  approaches  the 
problem  from  an  exact  modal  or  so  called  full  wave  solution. 

The  first  attempts  at  solving  for  e  ^  were  approximations  that 
attempted  to  model  the  microstrip  as  a  variation  of  a  simpler 
electromagnetic  problem.  Conformal  mapping  of  the  microstrip  structure 
was  used  [4],  which  transforms  the  geometry  of  the  microstrip  into  a 
simpler  coordinate  system  and  then  solves  the  quasi-static  boundary 
value  problem.  The  microstrip  was  also  modeled  as  a  ridged  waveguide 
and  the  solution  obtained  in  terms  of  longitudinal  side  electric  (LSE) 
modes  [5] .  More  recently  a  coupled  mode  approach  has  been  used,  where 
a  TEM  is  coupled  to  either  at  TMq  mode  [6],  [7]  or  a  TEq  mode  [8] -[10] 

to  solve  the  boundary  conditions. 

As  computers  became  more  inexpensive  and  prevalent,  methods  were 
introduced  that  solved,  numerically,  the  microstrip  boundary  value 
problem  in  spatial  variables.  These  methods,  known  as  spatial  domain 
techniques,  achieve  answers  in  the  form  of  a  set  of  simultaneous  linear 
equations  that  must  be  solved  with  a  matrix  inversion.  The 
formulations  are  exact  in  the  limit  as  the  number  of  equations  are 
increased,  but  large  amounts  of  computer  time  are  required  to  achieve 
high  accuracy.  One  spatial  domain  approach  is  to  formulate  a  singular 
integral  equation  for  the  microstrip  [11],  [12]  which  results  in  the 

system  of  linear  equations.  Another  approach  is  to  use  the  finite 
difference  method  [13],  [14]  which  uses  Maxwell's  equations  in  its 


5 


first  order  differential  form  to  specify  the  fields  at  a  set  of 
discrete  points.  The  first  order  derivatives  are  calculated  using  a 
difference  approximation  to  create  the  set  of  linear  equations. 

Dekelva  [15]  used  a  moment  method  approach,  enforcing  the  boundary 
conditions  at  a  discrete  number  of  points  to  achieve  results  for  c 

ref  f 

Mode  matching  [16],  [17]  splits  the  microstrip  into  different  regions 

where  the  boundary  conditions  are  uniform.  The  fields  in  each  region 
are  expanded  into  an  infinite  set  of  modes  that  satisfy  the  boundary 
conditions  of  the  particular  region.  Then  the  continuity  of  the  fields 
between  regions  is  enforced  by  equating  the  modal  sets  at  the 
boundaries  of  each  region. 

In  an  effort  to  reduce  the  amount  of  computer  time  used  by  the 

spatial  domain  techniques  while  retaining  the  accuracy  of  an  exact 

formulation,  the  Spectral  Domain  Approach  (SDA)  was  introduced  [18] . 

The  SDA  begins  the  solution  of  the  boundary  value  problem  by  Fourier 

transforming  the  spatial  domain  fields  into  the  spectral  domain.  This 

transformation  greatly  simplifies  the  formulation  of  the  boundary  value 

problem  and  allows  accurate  results  with  minimal  computer  time.  Early 

z  z 

work  with  the  SDA  approached  the  problem  assuming  TE  and  TM  modal 
configurations  for  the  fields  [19] -[21],  While  this  approach  yields 
the  correct  answers,  as  guaranteed  by  the  uniqueness  theorem,  its 

.  .  y 

derivation  is  lengthy  and  complex.  To  simplify  this  derivation,  TE 
(LSE)  and  TMy  (LSM)  modal  configurations  were  also  used,  [22] ,  [23] , 

since  the  microstrip  structure  closely  resembles  the  partially  filled 
dielectric  waveguide.  Using  this  modal  configuration,  the  derivation 


is  simplified  and  there  is  less  chance  for  error.  The  SDA  has  become 
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perhaps  the  most  popular  technique  for  solving  microstrip  problems  and 
in  1985,  IEEE  Transactions  on  Microwave  Theory  and  Techniques  published 
an  invited  paper  [24]  which  contains  a  relatively  complete  history  and 
description  of  the  SDA  as  well  as  an  extensive  bibliography. 

1.3  Coupled  Microstrips 

While  research  of  single  isolated  microstrips  is  important,  the 
effects  of  neighboring  conductors  must  also  be  considered  in  practical 
circuit  design.  The  simplest  structure  that  can  be  used  to  study  these 
effects  is  the  symmetric  coupled  microstrip,  where  the  presence  of  the 
second  conductor  adds  another  degree  of  complexity  to  the  solution  of 
the  boundary  value  problem.  However,  it  is  possible  to  split  this 
problem  into  two  simpler  ones  that  can  easily  be  solved  by  utilizing 
the  concept  of  even  and  odd  modes,  proposed  by  Wheeler  in  1956  [25] . 
This  method  is  described  in  detail  in  Chapter  2.  When  the  structure  is 
not  symmetric,  the  problem  can  be  split  into  a  similar  pair  of  modes, 
usually  referred  to  as  c  and  n  modes  [26] . 

As  with  the  single  microstrip,  early  research  centered  on 
obtaining  easily  computed  approximate  equations  to  describe  the 
frequency  dependence  of  e  Many  of  the  approaches  used  for  single 

microstrips  have  also  been  used  for  the  computation  of  coupled 
microstrips,  such  as  conformal  mapping  [27],  an  LSE  capacitance  model 
[28],  and  a  parallel  plate  waveguide  model  [29].  In  addition,  some 
curve  fitted  formulas  have  also  been  presented,  using  both  experimental 
data  [30]  and  results  computed  with  the  spectral  domain  approach  [31]. 
Although  these  approximate  formulas  are  useful,  their  application  is 
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limited  to  a  few  specific  cases.  To  be  able  to  consider  more  general 
structures  with  high  accuracy,  full  wave  methods  are  used,  as  in  the 
case  of  the  single  microstrip. 

Like  the  approximate  modelings,  exact  formulations  for  coupled 
microstrips  closely  parallel  the  derivations  for  single  microstrips. 
Also  like  the  single  microstrip  formulations,  the  boundary  value 
problem  can  be  solved  in  either  the  spatial  or  spectral  domain. 
Solutions  for  the  spatial  domain  are  given  by  Krage  and  Haddad  [32]  by 
expanding  the  fields  into  an  infinite  set  of  modes  in  each  medium  and 
enforcing,  numerically,  the  vanishing  of  the  electric  field  on  the 
center  conductors.  The  spectral  domain,  however,  is  the  easiest  of  all 
methods  to  extend  the  solution  from  a  single  conductor  to  multiple 
conductors  because  it  uses  a  Green's  function  approach.  This  approach 
formulates  the  boundary  value  problem  assuming  only  a  point  source  at 
first.  The  solution  for  the  actual  structure  is  obtained  by  using 
superposition  of  the  point  sources  to  represent  the  actual  current 
densities.  Thus  the  Green's  function  derived  for  the  single  microstrip 
can  be  used  without  modification  for  the  coupled  microstrip  problem. 

The  spectral  domain  has  been  applied  to  single  [33]  and  multiple 
substrate  [34]  coupled  microstrip  structures. 

1.4  Pulse  Dispersion 

For  single  frequency  (CW)  operation  on  a  microstrip,  it  is 

sufficient  to  know  only  a  single  value  of  e  which  can  then  be  used 

ref  f 

to  compute  most  of  the  other  parameters  of  interest.  However,  if 
either  non-periodic  or  wide  band  signals  are  to  be  used  on  the 
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microstrip,  then  it  is  important  to  study  the  effects  of  the  frequency 
dependence  of  creff  on  the  signal.  Since  the  microstrip,  as  a 
transmission  line,  is  dispersive  (the  phase  velocity  is  a  function  of 
frequency) ,  its  response  must  be  analyzed  in  the  frequency  domain, 
whereas  a  component  whose  response  in  time  is  non-linear,  such  as 
transistors  or  diodes,  must  be  analyzed  in  the  time  domain.  There  are 
three  different  methods  for  determining  the  time  domain  response  of  a 
dispersive  transmission  line; 

1)  Represent  the  transmission  line  as  a  set  of  matrix  parameters, 
such  as  ABCD  or  S  parameters. 

2)  Compute  the  inverse  Fourier  transform  of  the  product  of  the 
frequency  domain  transfer  function  of  the  line  and  the  Fourier 
transform  of  the  signal. 

3)  Model  the  structure  using  the  time  domain  finite  difference 
method. 

The  transmission  line  parameter  approach  has  been  used  in 
conjunction  with  the  microstrip  [35],  [36],  where  the  matrix  parameters 
were  determined  using  approximate  formulas.  Although  this  approach  is 
useful  for  cases  where  the  transmission  line  has  unmatched  loads, 
accurate  computation  of  the  matrix  parameters  requires  additional 
effort  that  is  not  necessary  with  the  Fourier  transform  method.  The 
time  domain  finite  difference  method  [14]  is  similar  to  the  finite 
difference  approach  used  to  find  the  effective  dielectric  constant, 
except  that  the  time  derivatives  approximated  as  well  as  the  spatial 
derivatives.  This  additional  degree  of  freedom  greatly  increases  the 
complexity  of  the  derivation,  as  well  as  the  numerical  expense. 
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The  Fourier  transform  method,  using  both  approximate  formulas  and 

full  wave  solutions  to  determine  e  „  ff),  has  been  widely  used  for 

ref  f 

determining  the  time  domain  response  of  microstrip  and  coplanar 

waveguide  (CPW)  structures  to  non-periodic  signal  inputs.  Using 

approximate  formulas  for  e  „„(f)  sacrifices  the  accuracy  of  the  full 

reff 

wave  solution,  but  drastically  reduces  the  required  computer  time. 
Graphs  of  distorted  pulses  using  approximate  formulas  for  the  single 
microstrip  have  been  presented  [3],  [37] -[41],  with  experimental 

results  shown  in  [38],  [42].  CPW' s  have  also  been  examined  with 

approximate  formulas  [43] .  Numerical  predictions  of  pulse  distortion 
using  full  wave  solutions  have  been  investigated  for  both  the  single 
microstrip  [44] -[46]  and  for  CPW' s  [47],  [48].  However  pulse 

distortion  on  coupled  microstrips  has  not  been  examined  from  a  full 
wave  perspective,  or  for  very  general  structures. 

1 . 5  Scope  of  Research 

Decreasing  computational  costs  coupled  with  increasing 
computational  speed  have  made  the  rigorous  analysis  of  microstrip 
boundary  value  problems  a  practical  choice  for  calculating  highly 
accurate  solutions.  This  report  formulates  full  wave  solutions  of  the 
boundary  value  problem  to  achieve  expressions  that  can  be  solved  for 
the  frequency  dependent  parameters  of  very  general  microstrip 
structures.  The  versatility  of  the  derivations  is  emphasized,  allowing 
many  different  microstrip  circuit  problems  to  be  solved.  Since  the 
accuracy  of  the  solutions  was  the  prime  concern,  extensive  optimization 
of  the  computer  programs  in  order  to  reduce  run  time  was  not  attempted. 
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While  there  are  many  single  frequency  applications  for 
microstrips,  the  propagation  of  non-periodic  signals  on  printed  circuit 
lines  is  also  an  important  concern,  and  so  the  degradation  of  pulsed 
signals  is  considered  in  this  report.  The  use  of  the  time  domain 
Fourier  transform  for  the  calculation  of  transient  response  is 
discussed,  with  the  emphasis  on  application  to  non-periodic  signals. 

The  Fourier  transform  method  is  then  applied  in  conjunction  with 
results  from  other  sections  to  predict  pulse  distortion  on  complex 
microstrip  structures.  The  types  of  complex  microstrip  structures 
dealt  with  include; 

1)  Multiple  dielectric  substrates  and  superstrates . 

2)  Symmetric  coupled  conductors. 

3)  Shielded  and  open  configurations. 

The  derivations  for  pulse  distortion  and  other  parameters  are  done  in  a 
efficient  and  general  manner,  so  that  the  computations  are  both 
flexible  and  expedient. 

The  coupling  between  adjacent  center  conductors  on  the  microstrip 
structure  is  analyzed,  since  such  results  are  vital  for  any  realistic 
microstrip  circuit  design.  The  concept  of  even  and  odd  modes  is 
introduced  and  is  then  extended  so  that  any  number  of  symmetrically 
center  conductors  can  be  analyzed.  This  approach  emphasizes  the 
versatility  of  the  spectral  domain  approach,  because  no  modification  of 
the  Green's  function  is  necessary  for  computation.  Results  for  pulse 
dispersion  on  multiple  conductors  are  presented,  showing  not  only  the 
signal  degradation  due  to  coupling,  but  also  the  effect  of  crosstalk  on 


adjacent  lines. 
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Since  modern  microwave  circuit  design  involves  structures  that  are 
much  more  complicated  than  the  simple  single  conductor,  single  layer, 
open  microstrip,  the  Green's  function  is  formulated  in  a  very  general 
manner  so  that  it  is  applicable  to  any  planar  geometry.  First,  the 
spectral  domain  approach  is  discussed  for  a  simple  structure  to 
illustrate  how  the  Green's  function  is  used  to  calculate  the  microstrip 
parameters.  Next,  The  Green's  function  for  an  arbitrary  planar 
microstrip  structure  is  derived  using  a  rigorous  full  wave  solution, 
giving  a  simple  recurrence  formula  for  its  calculation.  Using  this 
Green's  function,  calculations  of  the  effective  dielectric  constant  for 
multiple  dielectric  structures  are  given,  showing  the  effect  of  cover 
layers  and  substrate  layers  on  the  frequency  dependence  of  both 
multiple  and  single  lines.  These  results  will  lead  to  the  design  of 
multi-conductor  structures  that  have  small  or  no  coupling  between 
lines.  Finally,  pulse  dispersion  on  these  lines  is  presented,  showing 
the  decrease  in  both  signal  distortion  and  crosstalk. 

Using  exact  formulations  whenever  possible,  highly  accurate 
results  for  microstrip  parameters  and  pulse  distortion  are  computed. 

The  time  domain  Fourier  transform  is  used  to  calculate  pulse  distortion 
on  complex  microstrip  structures.  The  coupling  between  adjacent  center 
conductors  is  considered  using  an  even/odd  mode  formulation.  Using  a 
rigorous  boundary  value  method,  a  generalized  Green's  function  is 
derived  for  arbitrary  planar  structures.  Numerical  results  for  pulse 
distortion  on  complex  microstrip  structures  are  presented  including 
shielded  and  open  configurations,  multiple  dielectric  layers,  and 


coupling  with  adjacent  center  conductors. 


CHAPTER  2 

FULL  WAVE  ANALYSIS  OF  MICROSTRIP 


This  chapter  considers  the  microstrip  structure  in  an  exact 
formulation.  First  the  various  types  of  lines  that  can  be  treated 
using  this  approach  are  described,  and  the  geometry  used  for  the 
derivation  is  presented.  Next,  the  Fourier  transform  used  in  the 
derivation  is  introduced  and  the  various  terms  and  conventions  used 
throughout  the  report  are  defined.  Then  the  spectral  domain  approach 
is  applied  to  a  simple  structure,  the  shielded  microstrip,  to 
illustrate  the  use  of  this  method.  This  derivation  is  then  modified 
to  better  handle  slot  line  type  structures.  Finally,  a  simple 
recursive  formulation  for  the  Green's  function  of  a  generalized  planar 
structure  is  presented. 

2.1  The  Microstrip  Structure 

There  are  many  different  types  of  printed  circuit  transmission 
lines,  each  appropriate  for  different  applications,  all  of  which  can  be 
solved  using  the  Spectral-Domain  Approach  (SDA) .  Four  types  of  these 
transmission  lines  are  shown  in  Fig.  2.1;  microstrip,  coupled 
stripline,  coplanar  waveguide  (CPW) ,  and  slot  line.  Of  the  four, 
microstrips  have  the  lowest  dispersion,  high  Q,  and  highest  power 
handling  capabilities,  and  while  series  mounting  of  devices  is  easy, 
shunt  mounting  is  very  difficult.  On  the  other  hand  the  CPW  can  mount 


(a)  Microstrip 


(b)  Coplanar  Waveguide 


(c)  Slot  line 


(d)  Coupled  Striplines 


Fig.  2.1.  Common  planar  waveguide  structures 
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both  series  and  shunt  elements  easily  and  can  have  a  high  Q  [49], 
however  it  has  average  power  handling,  and  dispersion  characteristics. 
The  slot  line  allows  easy  mounting  of  shunt  elements,  but  has 
difficulty  with  series  components.  In  addition,  it  suffers  from  high 
radiation  losses,  low  Q,  low  power  capabilities,  and  is  highly 
dispersive.  Finally,  the  coupled  strip  line  offers  a  wide  range  of 
characteristic  impedances,  but  has  a  low  Q,  average  power  handling,  and 
average  dispersion  characteristics. 

The  microstrip  geometry  used  for  the  derivation  of  the  SDA's 
Green's  function,  as  well  as  for  the  computation  of  the  current  density 
expansion  functions,  is  shown  in  Fig  2.2.  A  standard  right-handed 
cartesian  coordinate  system  is  used  (the  z  axis  is  out  of  the  paper) , 
and  the  structure  is  considered  to  be  uniform  in  the  z  direction, 
reducing  it  to  a  two-dimensional  problem.  There  are  N  planar 
dielectric  layers  below  the  center  conductor  interface  and  m  layers 

above  it,  where  the  height  of  the  1th  layer,  upper  or  lower,  is  denoted 

,  ,  th 

by  h  or  h  respectively.  The  electrical  parameters  of  each  i 
Ui  Li 

layer  are  denoted  by  e  and  p  ,  representing  the  permitivity  and 
permeability,  respectively.  The  dielectric  materials  are  assumed  to  be 
isotropic,  so  both  p  and  c  are  scalars,  although  p  is  not 
necessarily  equal  to  the  free  space  value.  The  dielectric  constant, 
e,,  is  considered  initially  to  be  pure  real  (lossless  case),  but  it 
could  be  generalized  to  be  complex  (lossy  case) .  The  microstrip  is 
surrounded  on  all  four  sides  by  perfect  electric  conductors  (PEC's),  as 
in  a  shielded  configuration.  If  the  particular  microstrip  structure  is 


Fig.  2.2.  Geometry  for  general 
microstrip  structure. 
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to  have  no  side  walls,  then  a  - >  oo,  while  if  it  has  no  cover  sheet  or 

no  ground  plane,  then  h  or  h]  ^  - >  oo,  respectively. 

Because  the  SDA  is  formulated  for  a  unit  impulse  current  source  at 
an  arbitrary  x  location,  any  arrangement  of  conductors  that  are  uniform 

in  the  z  direction  may  be  considered  through  superposition.  In  this 

.  .  th  ,  ... 

report,  the  position  of  the  i  signal  conductor  is  determined,  as  in 

Figure  2.2,  by  the  position  of  the  center,  x  ,  and  the  width,  w  .  In 

addition,  each  of  the  center  conductors  is  assumed  to  be  of  zero 

thickness  so  that  the  formulation  of  the  problem  doesn't  become 

excessively  complex.  Although  there  may  be  any  finite  number  of  these 

center  conductors  on  the  interface,  for  clarity,  only  one  is  shown. 


2.2  The  Fourier  Transform 

Throughout  the  paper  an  e+"|Wt  time  convention  is  used,  and  the 
complex  wave  number  in  the  ith  medium  and/or  direction  is  designated  by 

y  ,  with  a  real  part  a  and  an  imaginary  part  (3  .  The  two-dimensional 

i  i  i 

Fourier  transform  pair  used  with  the  SDA  is  defined  as: 


G(/3x,y,/3z) 


jg(x,y,  z)  j  = 

-  00  _ 

-  J-  f- 

1|G(Px,y,Pz)|  =  g(x,y,z) 


j  (0  x+/3  z) 

g (x, y , z )  e  dx  dz 


(2.1) 


(2n) 


-j  (/3  x+p  z) 

G<P  ,y,P  )  e  *  Z  dp  d/3 

00  X  Z  X  z 


(2.21 


where  8  and  8  are  the  wave  numbers  in  the  x  and  z  directions.  If  the 
x  z 


structure  under  consideration  is  finite  in  the  x  direction,  i.e.  a  *  m 
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then  the  wave  numbers  in  the  x  direction  are  discrete  instead  of 
continuous,  and  a  slightly  different  transform  pair  is  used; 


G(Pxn'^z' 


fjgU,y,z)|  = 

„  00  „ 

-  U  l 

&  |  =  g(x»y»z) 


j(^xnX+PzZ) 

g (x,y, z)  e  dx  dz 


(2.3) 


(2  n)  n=-oo'-J 


m  g(/3  ,y,/3  )e 

oo  xn  z 


-j/3zz 


where 


3. 


xn 


(2n-l )  71 
2a 


A/3  =  n/a 

xn 


-j  13  x 

dMe  xn^x 

(2.4) 


(2.5) 

(2.6) 


This  transform  is  equivalent  to  using  a  Fourier  series  expansion  in  the 
x  direction  instead  of  a  Fourier  transform.  All  Fourier  transform 
quantities  are  denoted  with  a  ~  superscript  while  vectors  are 
underlined. 

This  definition  of  the  Fourier  transform  is  not  the  same  as  the 
one  that  is  normally  found  in  the  literature  since  the  exponent  of  the 
exponential  in  (2.1)  is  positive  instead  of  negative.  It  is,  however, 
equally  valid  and  allows  the  transform  of  the  partial  derivative  with 

-D/3zz 

respect  to  z  to  be  equivalent  to  assuming  an  e  variation,  as  is 


conventionally  used;  i.e., 
rd  f(z) 


f-§H  -  -j 


jP  F(P  )  < - > 

z  z 


M- 


■jfiz 


-j|3  z 

=  -j/3  e  (2.7) 

z 


While  this  new  definition  changes  some  of  properties  of  the  Fourier 
transform,  as  well  as  some  of  the  transform  pairs,  there  is  a  simple 
relation  between  the  two  definitions.  If  the  transform  of  f (z)  using 
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the  more  common  definition  is  given  by  F(/3  ),  then  the  transform  using 

z 

(2.1)  is  given  by  F  (-(3^)  . 

2.3  The  Spectral-Domain  Approach 

The  Spectral-Domain  Approach  (SDA)  is  preferred  over  other  full 
wave  solutions  for  planar  circuits,  especially  spatial  domain  methods, 
because  of  its  simple  formulation,  high  versatility,  low  memory 
requirements,  short  computation  time,  and  accurate  results.  The 
derivation  of  the  SDA  is  simpler  in  many  ways  than  spatial  techniques. 
For  example,  for  the  geometry  of  Figure  2.2,  the  partial  derivatives 
with  respect  to  x  and  z  in  spatial  coordinates  become  multiplications 
in  the  spectral  domain.  The  Green's  function  approach  used  in  the  SDA 
makes  it  applicable  to  many  different  types  of  planar  circuits 
including  the  following: 

a.  An  unlimited  number  of  signal  conductors  confined  to  one  plane. 

b.  An  unlimited  number  of  dielectric  layers  above  and  below  the 
interface  containing  the  conductors. 

c.  Top  (cover)  and/or  bottom  (ground  plane)  conductor  (either 
present  or  absent) . 

d.  Conductor  side  walls  (either  present  or  absent) . 

Unlike  spatial  domain  formulations,  large  matrices  are  not 

required  to  obtain  accurate  results  with  the  SDA.  In  addition,  large 
structures,  such  as  open  microstrips,  are  handled  as  easily  as  smaller 
ones,  without  the  huge  increase  in  memory  requirements  as  there  is  with 
spatial  methods.  Because  the  matrices  used  by  the  SDA  are  much  smaller 
and  fewer  numerical  integrations  are  needed,  computer  run  time  is  much 
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less  than  that  of  other  formulations.  Since  the  SDA  is  a  full  wave 
formulation,  the  results  are  very  accurate,  and  it  has  even  been  used 
to  produce  a  reference  standard  [50] . 


2.3.1  The  Wave  Equations  in  the  Spectral  Domain 

The  SDA  begins  by  transforming  the  fields  and  boundary  conditions 
from  the  spatial  domain  to  the  spectral  domain  in  order  to  simplify  the 
derivation  process.  Using  a  transform  to  solve  the  second  order 
differential  equation  is  equivalent  to  using  an  eigenfunction 
expansion,  but  the  formulation  is  more  straight  forward.  As  indicated 
previously  in  Section  2.2,  the  Fourier  transform  changes  both  the  z  and 
x  partial  derivatives  into  multiplications,  simplifying  the  electric 
and  magnetic  wave  equations  as  well  as  the  vector  potential  approach. 
Transforming,  as  in  (2.1),  the  time-harmonic  electric  and  magnetic  wave 
equations  for  lossless  media  with  only  a  current  source  present  gives; 


V2E  +  p2 E  =  -jw/iJ  < - » 


V2H  +  £2H  =  -VXJ  < - » 


—  +  (32 

dy2  y 


+  P2 

dy2  y 


E  =  -jWfiJ 


(2.8) 


+  a  j 

y 


p  j  -  p  j 

z— X  z— z 


H  =  -a 


+  a 


j0  J  +  -jt  J 
z-  y  dy  —  z 


jp  J  +  -Q-  J 
x-y  fly  —  x 


(2.9) 


Since  the  microstrip  conductor  is  assumed  to  be  infinitely  thin,  there 
is  no  current  density  in  the  y  direction  and  the  partial  derivatives 
with  respect  to  y  of  the  x  and  z  current  densities  are  zero.  Using 
this,  the  transform  of  the  magnetic  wave  equation  can  be  written  more 


concisely  as: 
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—  +  /32 

aY2  Y 


H  =  j  a 


[/3  J  “  P  J  1 

'  I  z— x  z— zl 


(2.10) 


The  transformed  wave  equations  are  second  order  partial 

differential  equations  of  only  a  single  variable,  y,  and  so  they  give 

no  information  or  constraint  on  the  x  or  z  variations.  The  x  and  z 

dependence  of  the  problem  is  now  contained  implicitly  in  the  transforms 

of  either  the  current  densities  or  the  electric  fields  through  the 

functional  dependence  on  (3  and  p  .  To  illustrate  the  use  of  the  SDA, 

a  simple  structure,  the  shielded  microstrip  is  considered.  First,  the 

y  y 

derivation  for  TM  mode  is  done,  and  the  TE  follows. 


2.3.2  TM^  Modes  for  the  Shielded  Microstrip 

The  geometry  of  the  shielded  microstrip,  shown  in  Figure  2.3,  is 

the  equivalent  to  the  geometry  found  in  Figure  2.2  with  n  =  M  =  1, 

h  =h,h  =h,x=0,  and  w  =  w.  Since  the  structure  has  walls 
ui  2  LI  1  i  i 

at  a  finite  distance,  the  second  Fourier  transform  is  used  as  stated  in 
(2.3)  through  (2.6).  This  derivation  can  also  be  generalized  to  the 
open  microstrip  structure  by  letting  a  — »  oo  and  h^  — >  co  and  using  the 
first  Fourier  transform  pair  of  (2.1)  and  (2.2).  Transforming  the 
fields  and  the  vector  potential.  A,  into  the  spectral  domain  gives: 


A  = 

a  A  (x,y, z) 

y  y 

— > 

IV  .  IV 

A  =  a  A  (13 
-  y  y  x 

,y,Pz> 

(2.11) 

E™ 

.  i 

32  A 

13 

g™  _  x 

(3  iv 

_  A 

(2.12a) 

X 

dxdy  i 

/ 

x  w/ic 

ay  y 

E™ 

y 

.  1 

4*; 

py 

A 

y 

< — 

- 

~TM  .1 

E  =  -j - 

y  cope 

'-V 

py2 

A~ 

y 

(2.12b) 

E™ 

_ .  1 

a2  A 

13 

~TM  Z 

E  —  — 

d  ~ 

_  A 

(2.12c) 

z 

■'wpe 

dzdy  } 

/ 

< 

> 

Z  W/1C 

Sy  y 
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TM 


^  A 
M  dz  y 


P 

~T  M  ,  Z  ~ 

-»  H  =  3  -  A 

x  mv 


H™  =  0 

y 


~TM 

< - >  H  =  0 

y 


h’m  -  !  3§  » 

z  n  ox  y 


~TM  ,  X  ~ 

< - »  H  =  -  j  -  A 

z  My 


rk;  t  h 

The  scalar  potential  A  (6  ,y,fi  )  for  the  i  layer, 

yi  x  z 


(2 . 12d) 

(2 . 12e) 

(2. 12f) 
1  or  2 , 


is  chosen  as  follows: 

REGION  1:  (0  s  y  s  h  ) 

Ayi  =  Bi(Px^z)  cosh(ayiy)  (2.13) 

REGION  2:  (h  <  y  <  h  ) 

1  J  2 

A  =  BO  ,0  )  cosh  [a  (h  +h  -y)]  (2.14) 

y2  2  x  z  y2  2  1 

B  O  ,|3  )  and  B  (/3  ,  f3  )  are  unknown  functions  of  the  transform 

1  X  Z  2  X  Z 

variables,  which  are  to  be  determined  from  the  boundary  conditions. 

They  are  abbreviated  as  B  and  B  throughout.  A  and  A  are  chosen 

based  on;  1)  the  tangential  E  fields  vanishing  at  y  =  0  (at  the  ground 

plane)  and  2)  the  tangential  E  fields  vanishing  at  y  =  h^+h^.  The 

hyperbolic  cosine  is  chosen  since  the  structure  should  have  standing 

waves  in  the  y  direction, 

y 

If  the  TM  mode  is  an  independent  solution  of  the  boundary  value 
problem,  then  either  the  x  or  the  z  directed  fields  can  be  used  to 

enforce  the  boundary  conditions.  On  the  other  hand,  if  using  the  x  and 

.  •  .  y 

z  directed  fields  gives  different  answers,  then  the  TM  mode  is  not  a 

solution  and  the  problem  must  be  solved  as  a  superposition  of  two 

orthogonal  modes.  Using  the  scalar  potentials  in  (2.13)  and  (2.14) 

with  the  equations  for  the  fields  in  (2.12a-f),  the  electric  and 

magnetic  fields  in  each  region  are  calculated. 
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REGION  1: 


0  a 

~  ~  x  vi 

E  =  -b  - i-  sinh  (a  y) 

xl  1  yl 

0  a 

iv  i'v  Z  V  l 

E  =  -B  -  sinh(a  y) 

zl  l  w/^e  yl 

/3z 

H  =  j  B  —  cosh (a  y) 
xl  J  l  yl 

Px 

H  =  -j  B  —  cosh (a  y) 
zl  J  l  ^  yr 


where 


2  2  2  2  2 

0  =  0  +  0  -  a  =  w  V  c 

*1  x  '  z  yl  *1  1 


REGION  2: 


0  a 

~  ~  X  V2 

E  =  B  - — —  sinh  [a  (h  +h  -y)  ] 

x2  2  y2  1  2 

0  a 

~  ~  z  v2 

E  =  B  - —  sinh  [a  (h  +h  -y)  ] 

z2  2  uM2c2  y2  1  2 

0Z 

H  =  j  B  —  cosh [a  (h  +h  -y) ] 
X2  J  2  y2  1  2 

0 

H  =  -j  B  —  cosh[a  (h  +h  -y)  ] 
Z2  J  2  n2  y2  1  2  J 


where 


_2  _  2  2  2  2 

0  =  0  +  0  -  a  =  w  n  c 

^2  '  x  '  z  y2  ^2  2 


(2.15a) 

(2 ,15b) 

(2.15c) 

(2 . 15d) 

(2.16) 

(2.17a) 

(2.17b) 

(2.17c) 

(2. 17d) 

(2.18) 


Now  that  the  expression  for  the  fields  have  been  specified,  the 


boundary  conditions  can  be  applied  to  eliminate  B  and  B^ . 


1) 


xl 


u  =  E 

y=h  x2 
J  1 

0  « 


y=h 


_  0  a 

3  x  yl  sinh  (a  h  )  =  B  x  y2-  sinh  (a  h  ) 
l  Ufi  c  yl  l  2  Ufi  c  y2  2 


l  l 


2  2 


a  c  ii  sinh  (a  h  ) 

b  =  -  b  y2  rl  rl  y2  2 

l  2  a  e  n  sinh  (a  h  ) 


yl  r2  r2 


yi  i 


(2.19) 

(2.20) 


The  z  directed  electric  fields  are  now  considered. 
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2)  E 


zl 


=  E 

y=hj  z2 


y=h. 


Pa  pa 

~  Z  vl  .  ~  Z  V2 

-B  - —  sinh  (a  h  )  =  B  - —  sinh(a  h  ) 

l  u(i  e  yl  l  2  a>fi  e  y2  2 


l  l 


2  2 


a  c  u  sinh  (a  h  ) 

B  =  -  B  y2  rl  rl  y2  2 

l  2  a  e  u  sinh (a  h  ) 

yl  r2^r2  yl  1 

This  gives  the  same  result  as  the  x  fields,  as  desired. 


(2.21) 


(2.22) 


Next,  the  discontinuity  of  the  H  fields  at  the  center  conductor 


interface  due  to  the  presence  of  a  current  source  is  enforced. 


3)  H 


xi 


y=hl 

P_ 


X2 


~TM 

y=h  z 

J  1 


P 

.  -  Z  .  ~  Z  ~TM 

i  B  —  cosh  (a  h  )  -  -j  B  —  cosh  (a  h  )  =  J 
l  P1  yi  1  2  n2  y2  2  z 


(2.23) 


Substituting  in  for  B^  using  (2.20) 

p  a  c  P 

.  ~  Z  V2  rl  .  .  ~  Z  ~TM 

-j  B  -  coth  (a  h  )  sinh  (a  h  )  -j  B  —  cosh  (a  h  )  =  J 

J  2  u  a  e  yil  y2  2  J  2  u  y2  2  z 

r-2  yl  r2  2 

(2.24) 


e  coth (a  h  )  c  coth (a  h  ) 
rl  yl  l  r2  y2  2 


a 


yi 


Let 


a 


y2 


.  C  r2^  2  ~TM 

=  3  - 5-  J  csch  (a  h  )  (2.25) 

a  P  z  y2  2 

y2  z  J 


c  coth (a  h  ) 
-IH  rl  yl  l 


si 


a 


yi 


c  coth (a  h  ) 

~TM  r2  y2  2 

s2  a 

y2 

This  simplifies  (2.25)  to 
c  n 

~  r2  2  ~TM 

B  =  ]  - 5—  J  csch  (a  h  ) 

2  J  a  p  z  y2  2 

y2  z  J 


f~TM  +  ~Tm1  -1 
L  Si  S2j 


(2.26a) 


(2.26b) 


(2.27) 


Next,  the  z  directed  magnetic  fields  are  checked  to  see  if  they 
give  the  same  answer. 

»v 

4)  H 


zl 


v.  -  H 

y=h  z2 

J  l 


~TM 

—  -  J 

y=h  x 

J  1 


P  P 

,  /v  x  ,  ^  X  a/'I’M 

-i  B  —  cosh  (a  h  )  +  i  B  —  cosh (a  h  )  =  -J 
1  (ij  yl  1  2  n2  y2  2  x 


(2.28) 
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Using  (2.20)  for  B  as  before, 

0  a  e  0 

.  ~  X  V2  rl  .  .  ~  X  ~TM 

-3  B  - - -  coth(a  h  )sinh(a  h  )  -j  B  —  cosh  (a  h  )  =  -J 

J  2  u  a  c  vl  1  y2  2  J  2  a  y2  2  x 

f-2  yl  r2  J  J  J 

(2.29) 

Simplifying  this  equation  using  (2.26a)  and  (2.26b)  leads  to 
c  u  r  -1  -1 

~  r2  2  ~TM  l~TM  ~TM  I 

B  =  j  V  j  csch  (a  h  )  Y  +  Y  (2.30) 

2  J  ay20x  x  y2  2  [  S1  s2J 


In  order  for  the  TM*  mode  to  be  an  independent  solution  in  the 
spectral  domain,  (2.27)  should  be  equivalent  to  (2.30).  This  occurs  if 
the  x  and  z  directed  modal  currents  are  related  as  follows; 


0  J™  =  0  J™ 

Z  X  X  z 


(2.31) 


To  show  that  this  relation  is  indeed  true,  the  magnetic  wave 
equation  with  only  an  electric  current  source  present,  as  in  (2.10)  is 


used: 


—  +  02 

ay2  y 


;  %  -  i  [1 


A/  ^ 

0  J  -  0  J 
z  x  z  z 


(2.32) 


Since  the  TM*  mode  has  no  H  ,  the  left  side  of  (2.32)  is  zero  and  so 

y 


0  j™  =  0  j™ 

Z  X  X  z 


(2.33) 


Therefore  the  Thr  mode  satisfies  all  the  boundary  conditions  in 
the  spectral  domain.  To  complete  the  solution,  the  tangential  fields 

/V 

at  the  center  conductor  interface  are  required.  Replacing  B^  in 
(2.17a)  and  (2.17b)  with  (2.27),  and  evaluating  them  at  y  *»  h  leads  to 


~TM  ,  l  x  ~TM  ~TM 

E  =  j  -  75—  Z  J 

X2  coc  0  z 

0  z 


(2.34) 


~TM  ,  1  ~TM  ~TM 

E  •“  n  -  Z  J 

z2  we  z 

0 


(2.35) 
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~TM  ,  , 

where  Z  is  the  modal  input  impedance,  defined  as 

-1 


~TM 

z 


PvTM  +  ~TM' 
[sl  S2 


(2.36) 


2.3.3  TE ^  Modes  for  the  Shielded  Microstrip 

y  y 

Having  solved  the  TM  mode,  the  TE  mode  is  considered. 


F  =  a  F  (x,y, z) 

y  y 


F  =  a  F  (13  ,y,0  ) 
—  y  y  x  z 


(2.37) 


TE 


I-»F 
e  3z  y 


-TE 


±  F 


(2.38a) 


TE 

E  =0 

y 


TE 


i  — —  f 

e  3x  y 


ute  ■  1  o 

x  Ju >pc  dxdy  y 


~TE 

E  =0 

y 


ft 


~TE  .  '  X  ~ 

E  =  j  -  F 

z  c  y 


hte 

X 


_fx  a_  ~ 

w/ie  3y  y 


(2.38b) 

(2.38c) 

(2 . 38d) 


hte 

y 

,  i 

4«2 

V 

F 

y 

- 

- 

hte  = 
y 

,  i 

W/1E 

i  „2 

py 

F 

y 

(2 . 38e) 

hte 

z 

.  1 

s2  F 

hte  = 

Z 

^Z 

_  XT’ 

(2 . 38f ) 

3z3y  i 

/ 

dy  y 

The 

t 

scalar  potential  F^ 

(Px' y'Pz> 

for  the 

•  th  1 

i  layer,  i 

=  l 

or  2, 

is  chosen  as  follows: 


REGION  1:  (0  s  y  s  h  ) 

F  =  A  (|3  ,/3  )  sinh  (a  y) 
yl  1  x  z  yl 

REGION  2:  (h  <  y  <  h  ) 

1  J  2 

F  =  A  (P  ,/3  )  sinh  [a  (h  +h  -y) ] 
y2  2  x  z  y2  2  1 


(2.39) 


(2.40) 


A  O  ,£  )  and  A  (f3  ,  (3  )  are  unknown  functions  of  the  transform 

1  X  Z  2  X  Z 

variables,  which  are  to  be  determined  from  the  boundary  conditions. 

They  are  abbreviated  as  A  and  A  throughout.  F  and  F  are  chosen 

l  2  yl  y2 


based  on;  1)  the  tangential  E  fields  vanishing  at  y  =  0  (at  the  ground 
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plane)  and  2)  the  tangential  E  fields  vanishing  at  y  =  h^+h^.  The 
hyperbolic  sine  is  chosen  since  the  structure  should  have  standing 
waves  in  the  y  direction. 

As  with  the  TM  mode,  if  the  TE  mode  is  an  independent  solution 

of  the  boundary  value  problem,  then  either  the  x  or  the  z  directed 

fields  can  be  used  to  enforce  the  boundary  conditions.  Likewise,  if 
,  y 

they  are  not  equivalent  then  the  TE  mode  isn't  an  independent  solution 
and  the  problem  needs  to  be  solved  with  a  superposition  of  two 
orthogonal  modes.  Using  the  scalar  potentials  in  (2.39)  and  (2.40) 
with  the  equations  for  the  fields  in  (2.38a-f),  the  electric  and 
magnetic  fields  in  each  region  are  calculated. 

REGION  1: 


E  =  -j  A  —  sinh (a  y) 
xl  J  l  c]  yi/ 

/3 

A/  ,  A/  X  , 

E  =  i  A  —  sinh (a  y) 
zl  J  l  e  yl 

l 

13  a 

~  ~  x  vi 

H  =  -A  -  cosh  (a  y) 

xl  l  WM1e1  yl 

/3  a 

~  ~  z  vl 

H  =  -A  -  cosh  (a  y) 

zi  i  yi 


where 


2  2  2  2  2 

0  =  P  +  0  -  a  -  an  c 

'l  'x  *z  yl  ^l  l 


REGION  2: 


0Z 

E  =  -j  A  —  sinh [a  (h  +h  -y) ] 

X2  J  2  C2  y2  1  2 

Px 

E  =  j  A  —  sinh[a  (h  +h  -y) ) 
Z2  J  2  C  y2  1  2 

2 

f 3  a 

~  ~  x  v2 

H  =  A  - —  cosh  [a  (h  +h  -y) ] 

X2  2  W/l 2C2  y2  1  2  1 


(2.41a) 

(2.41b) 

(2.41c) 

(2.4  Id) 

(2.41) 

(2.42a) 

(2.42b) 

(2.42c) 
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/3  a 

/v  z  V2 

H  =  A  - —  cosh  [a  (h  +h  -y)  ] 

z2  2  y2  l  2  J 


(2 . 42d) 


where 


2  2  2 

0=0  +  (3 

2  X  Z 


2  2 
a  =  w  u  e 

y2  2  2 


(2.42) 

Now  that  the  fields  have  been  determined,  the  boundary  conditions 


can  be  applied  to  eliminate  A^  and  A 


1) 


xl 


w  =  E 
y=hi  x2 


y=h. 


ez  Pz 

-i  A  —  sinh(a  h  )  =  -j  A  —  sinh  (a  h  ) 

J  l  e  yl  l  J  2  e  y2  2 

l  r  2  * 

c  sinh (a  h  ) 

~  ~  ri  y2  2 

p.  =  _  _ 

l  2  c  sinh (a  h  ) 

r2  yi  l 

The  z  directed  electric  fields  are  now  considered. 

2)  El  =  E 

zl | y=h^  z2 | 

0, 


(2.43) 


(2.44) 


y=h 


0x 

-j  A  —  sinh (a  h  )  =  -j  A  —  sinh (a  h  ) 
Jie  yl  l  2€  y2  2 

1  jt  2  r 

e  sinh (a  h  ) 

~  ~  rl  y2  2 

A  =  A  _ 

l  2  e  sinh (a  h  ) 

r2  yl  l 


(2.45) 


(2.46) 


which  is  the  same  result  as  using  the  x  fields,  as  required. 

Next,  the  discontinuity  of  the  H  fields  at  the  center  conductor 
interface  due  to  the  presence  of  a  current  source  is  enforced. 

3)  H 


xl 


u  -  H 
y=h  x2 

3  1 


y=h 


~TE 

J 

z 


0  a  0  a 

«  x  yl  ~  x  y2  ~te 

-A  - —  cosh  (a  h  )  -  A  - —  cosh  (a  h  )  =  J 

l  up  c  yl  l  2  wp  c  y2  2  z 


(2.47) 


l  l 


2  2 


Substituting  in  for  A^  using  (2.44)  leads  to 

fi  a  0  a 

~  x  yl  ,  ~  x  y2 

-A  -  coth  (a  h  )  sinh  (a  h  )  -  A  - =—  cosh  (a  h  )  =  J 

2  UfJL  C  yl  1  y2  2  2  W/J  c  y2  2  z 


TE 


1  2 


a  coth  (a  h  )  a  coth  (a  h  ) 
yl  yl  1  y2  y2  2 


2  2 


(2.48) 


rl 


r2 


Oli  c 

0  2  ~TE 

— ^ —  J  csch  (a  h  ) 
P  z  y2  2 


(2.49) 
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a  coth (a  h  ) 
'TE  _  yl _  yl  1 


(2.50a) 


a  coth (a  h  ) 
-TE  _  y2  y2  2 


(2 . 50b) 


This  simplifies  (2.49)  to 


wu  e 

A  =  -  — -  JTE  csch  (a 
2  pz  y2 

x 


V  [' 


~TE  +  ~TE 
Si  32 


(2.51) 


Next,  the  z  directed  magnetic  fields  are  checked  to  see  if  they 
give  the  same  answer. 


4)  H 


zi  y=h  z2  y=h 


~TE 
—  —  J 


Pa  pa 

~  X  Vl  ~  X  V2  ~TE 

-A  - —  cosh  (a  h  )  -  A  - i-  cosh  (a  h  )  =  -  J 

l  yl  l  2  y2  2  x 


(2.52) 


Using  (2.44)  for  A^  as  before, 

Pa  Pa 

«  X  Vl  .  ~  X  v2  ~T£ 

-A  - —  coth  (a  h  )smh(a  h  )  -  A  - —  cosh  (a  h  )  =  -J 

2  COfI  £  yl  1  y2  2  2  lOfl  c  y2  2  x 


Simplifying  this  equation  using  (2.50a)  and  (2.50b)  leads  to 


wu  e  r 

~  02  ~TE  ~TE  ~TE 

A  =  — J  csch  (a  h  )  Y  +  Y 
2  P  x  y2  2  si  s2 

z  L 


(2.53) 


(2.54) 


In  order  for  the  TEy  mode  to  be  an  independent  solution  in  the 
spectral  domain,  (2.51)  should  be  equivalent  to  (2.54).  This  occurs  if 
the  x  and  z  directed  modal  currents  are  related  as  follows; 


P  JTE  =  -  P  JTE 

Z  Z  XX 


(2.55) 


To  show  that  this  relation  is  indeed  true,  the  continuity  equation 


is  used  with  a  time  varying  volume  charge  density,  q(x,z;t),  equal  to 


,  v  jwt  .  , 

p^(x, z)e  ,  giving 


_  d  jwt 

Voj  =  -  _  ^pv(x,z)e 


(2.56) 
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Computing  the  time  derivative  and  suppressing  the  time  dependence 
reduces  (2.56)  to 

(2.57) 

Fourier  transforming  the  (2.57)  twice,  with  respect  to  x  and  z  as  in 


a  a  a 

_  J  +  —  J  +  —  J  =  -  JW  p  (x,z) 
ox  x  a  y  y  a  z  z  v 


(2.1)  and  letting  =  0  (due  to  zero  strip  thickness)  leads  to 
.  „  ~TE  .  „  ~TE  .  ~ 

-  ]P  J  -  j/3  J  =  -  jo >  p  (2.58) 

X  X  Z  Z  V 

y 

For  a  TE  field,  the  normal  E  field  (E^)  is  zero  everywhere.  The 
boundary  condition  of  n  «  eE  =  then  implies  that  p  (x, z)  =  0 
everywhere,  and  so  p  (|3  ,  0  )  =0  also.  Therefore  (2.58)  reduces  to 

V  X  z 


(2.59) 


Thus  the  TE*  mode  independently  satisfies  all  the  boundary 

conditions  in  the  spectral  domain.  To  complete  the  solution,  the 

tangential  fields  at  the  center  conductor  interface  are  required. 

Replacing  in  (2.42a)  and  (2.42b)  with  (2.51),  and  evaluating  them  at 

y  =  h  leads  to 
l 


~TE  .  Z  ~TE  ~TE 

E  =  ]  toil  Z  J 

x2  J  *0  8  z 

x 


~TE  ,  ~TE  ~TE 

E  =  -]  UU  Z  J 
Z2  0  Z 


(2.60) 


(2.61) 


~te  ,  y 

where  Z  is  the  TE  modal  input  impedance,  defined  as: 


zTE  =  |VE  +  ^TE1  1 

[si  S2j 


(2.62) 


The  tangential  electric  fields  at  the  center  conductor  interface 
have  been  determined  for  each  mode  in  terms  of  the  modal  current 
densities.  However,  since  solution  of  the  microstrip  problem  requires 
the  total  fields  at  the  interface,  the  modal  fields  of  (2.34),  (2.35), 
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(2.60)  and  (2.61)  are  added  together  to  form  the  x  and  z  directed 
fields  giving 


E 

X2 


/V 

E 


Z2 


~TM  .  ~TE 
E  +  E 
X2  x2 


~TM  ,  ~TE 
E  +  E 
Z2  z2 


(2.63a) 

(2.63b) 


2.3.3  Current  Density  Relations 

V  V 

To  combine  the  TE  and  TM  modal  solutions  there  needs  to  be  an 

.  .  ,  .  ~TE  ~TE  ~TM  ~TM 

equation  relating  the  modal  current  densities  J  ,  J  ,  J  ,  and  J  , 

z  x  z  x 

to  the  total  current  densities  J  and  J  .  Each  modal  solution,  TEy  or 

x  z 

y  .  . 

TM  ,  contains  only  one  independent  current  component,  the  other  is 
related  to  it  by  either  (2.33)  or  (2.59).  Therefore  the  combined  modal 
solution  for  the  microstrip  must  contain  two  independent  current 
distributions.  Since  the  total  fields  are  a  combination  of  the  modal 
fields,  the  total  current  densities  can  be  written  as  the  sum  of  the 
modal  current  densities;  that  is 


~TM  /vTE  ~ 

J  +  J  —  J 

z  z  z 

/vTM  /vTE  'v 

J  "i"  vJ  =  vJ 

XXX 


Substituting  (2.33)  and  (2.59)  into  (2.65)  leads  to 

P  P 

~  _  ~TM  X  ~TE  Z 

x  z  6  z  B 

z  x 

Multiplying  (2.64)  by  PZ^PX  reduces  to 

P  P  P 

~  Z  _  ~TM  Z  ~TE  Z 

z  P  JZ  f3  +  JZ  P 

XXX 


Adding  (2.67)  to  (2.66)  leads  to 

0 


~TM 

u 

Z 


^z  ^x 

p  +  p 

X  Z 


7T  +  J 

Z 


(2.64) 

(2.65) 


(2.66) 


(2.67) 


x 


(2.68) 
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Multiplying  both  sides  of  (2.68)  by  (3  /3  /  (p2  +  f3? )  ,  the  current  density 

X  Z  X  z 

~TM 

J  can  be  written  as 
z 


3t„  _  VA^A1 

z  O2+02) 

X  z 


(2.69) 


To  get  the  relationship  for  the  TEy  current  density,  multiply 


(2.66)  by  -p  //3  reduces  to 
z  x 

0 

~  Z  _  ~T  M  ~TE 

-J  -5—  -  ~ J  +  J 

X  B  Z  Z 

X 


T  2 


(2.70) 


Adding  (2.70)  to  (2.64)  leads  to 

.2n 


~TE 

J 

z 


1  + 


0. 


A/  /V  Z 

=  J  “  J  — 

Z  X  p 

X 


(2.71) 


2  2  2 

Multiplying  both  sides  of  (2.71)  by  0  /  (0  +  0  ),  the  current  density 

XX  z 

~TE 

J  can  be  written  as 
z 


A/  Ay 


~rE .  ^IJA-JA' 

Z  (0Z+02) 

X  z 


(2.72) 


2.3.5  Green's  Functions 

Having  obtained  relationships  between  the  modal  and  total  current 
densities,  the  total  electric  fields  at  the  interface  can  be  expressed 
in  terms  of  the  total  current  densities.  Using  (2.72)  and  (2.69)  to 
express  the  modal  currents  in  terms  of  the  total  currents  in  (2.63a) 
and  (2.63b),  the  total  fields  at  the  center  conductor  interface  can  be 


written  as 
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X2 


22 


2  2 

0+0 

x  z 


2  2 

0  +0 

X  z 


0  0 

X  Z  T  M  0  0  TE 

- Z  +  Wli  0  0  Z 

we  Mopzpx 


+J 


0 

X  TM  _2  TE 

- Z  -  Wli  0  Z 

we  'o  z 

0 


z 

V.  L 


’0 

Z  TM  _  2  TE 

- Z  -  Wli  B  Z 

we  o  x 

o 


+J 


0  0 

X1  Z  TM,  0  „  TE 

- z  +  wu  0  0  z 

we  o  z  x 


where 


E 

z2 

=  —  r  j  g  +j  g  ] 

WE  L  Z  ZZ  X  XzJ 

0 

=  —  [  J  G  +  J  G  1 
we  L  z  zx  x  xxJ 

0 

E 

X2 

r;  — 

02Z™-/320VE 

z  x'  0 

VJ 

Z  Z 

2  2 

0+0 

X  z 

rz  — 

„2  TM  _2  _2  TE 

0  Z  -0  0  Z 

X  z  0 

VJ 

XX 

2  2 

0+0 
x  z 

J  = 

XZ 

TM  „2  TE 

Z  +13  Z 

G  =00  - - - 

ZX  X  Z  02+02 

X  z 

(2.73a) 

(2.73b) 

(2.74a) 

(2.74b) 


(2.75a) 


(2.75b) 


(2.75c) 


Equations  (2.74a)  and  (2.74b),  for  a  given  frequency  and 

structural  parameters,  have  two  unknown  variables;  0  and  0  ,  and  four 

x  z 

unknown  functions;  J  (0  )  ,  J  (0  ) ,  E  (0  ,0  ),  and  E  ({3  , 13  )  . 

z  x  x  x  z2  x  z  x2  x  z 

Furthermore,  the  magnetic  and  electric  wave  equations,  after  being 
transformed  twice,  give  no  information  or  constraints  on  the  form  of 
the  four  functions,  as  it  does  in  the  spatial  domain.  It  is 
advantageous,  then  to  use  a  method  which  can  eliminate  the  four 
functions  and  one  variable. 


2.3.6  Galerkin's  Method 

Galerkin's  method  is  a  method  of  moments  technique  that  uses  the 
same  expansion  and  testing  functions  to  solve  an  integral  equation.  To 
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apply  this  method,  the  current  densities  are  expanded  into  complete 


sets  of  basis  functions 
oo 

J  =  Y  a  J 
z  ,  n  zn 

n=l 


(2.76a) 


j  =  y  b  j 

x  ,  m  xm 

m=l 


(2.76b) 


where  a  and  b  are  unknown  constants  that  specify  the  amplitudes  of 
n  m 

each  current  expansion  function. 

The  x  and  z  current  densities  in  the  electric  fields  of  (2.74a) 
and  (2.74b)  are  now  replaced  by  the  series  expansions  for  the  currents 
of  (2.76a)  and  (2.76b).  This  replaces  the  unknown  x  and  z  current 
densities  with  known  basis  functions  of  unknown  amplitude,  resulting 


E  =  y  a  J  G  +  r  b  J  G 
z2  .  n  zn  zz  ,  m  xm  zx 

n=l  m=l 


(2.77a) 


E  =  r  a  J  G  +  £  b  J  G 

x2  ,  n  zn  zx  ,  m  xm  xx 

n=l  m=l 


(2 . 77b) 


Now  (2.77a)  is  multiplied  by  the  complex  conjugate  of  (2.76a)  and 
(2.77b)  is  multiplied  by  the  complex  conjugate  of  (2.76b)  which  gives 


n  zn  z2 


m  xm  x2 


00 

00 

_ 

^  a  J 
n  zn 

l 

n=l 

n=l 

00 

00 

_ 

y  b  j 

^  m  xm 

I 

m=l 

n=l 

«V  /V/ 


<V  /V 


)  J  G 
m  xm  zx 


<v  /x/ 

3  J  G 
m  xm  xx 


(2.78a) 


(2.78b) 


Integrating  both  sides  of  (2.78a)  and  (2.78b)  from  minus  infinity  to 


plus  infinity  with  respect  to  gives 


n=l  -oo 

oo  r  oo 


J  E  d/3 
zn  z2  x 


00  00 


y  a  V  a  K  +  y  b  K 

,  n  ,  n  nn  ,  m  mn 

n=l  n=l  m=l 


n=l  -co 


00  00 


j  e  d/?  =  r  b  r  a  k  +  r  b  K 

xm  x2  x  ^ '  L~‘  r'rn  "> 


,  m  .  n  nm  ,  m  mm 

m=l  n=l  m=l 


(2.79a) 


(2.79b) 
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where  for  open  configuration  (a  =  oo) 


k  =  r 

nn  J 

K  =  f 

mn  J 

K  =  f 

nm  J 

K  =  [ 
mm  J 


J  (0  )J  (0  )  G  (0  ,0  )  d0 
zn  x  zn  x  zz  x  z  x 


-oo 

oo 


J  (0  )J  (0  )G  (0  ,0  )  d0 
zn  x  xm  x  zx  x  z  x 


J  (0  )J  (0  )G  (0  ,0  )  d0 
xm  x  zn  x  zx  x  z  x 


J  (0  )J  (0  )G  (0  ,0  )  d0 
xm  x  xm  x  xx  x  z  x 


(2.80a) 


(2.80b) 


(2.80c) 


(2. 80d) 


For  shielded  configurations  (a  *  oo)  the  integrations  are  replaced  by 

summations  since  0  is  a  discrete  instead  of  a  continuous  variable  and 
x 

d0^  is  replaced  by  A0^  as  in  (2.5)  and  (2.6).  Using  these  changes  in 
the  above  four  equations  gives 


00  * 

K  =  £  J  (0  .)J  (0  . ) G  (0  .,0  )  A0 

nn  .  ^  zn  xi  zn  xi  zz  xi  z  x 

1  =-oo 

00  * 

K  =  £  J  (0  ,)J  (0  . ) G  (0  .,0  )  A 0 

mn  ,  zn  xi  xm  xi  zx  xi  z  x 

1  =-00 

00  * 

K  =  £  J  (0  .) J  (0  ,)G  (0  .,0  )  A0 

nm  .  xm  xi  zn  xi  zx  xi  z  x 

1  =-00 


(2.81a) 


(2.81b) 


(2.81c) 


mm 


£  J  (0  , ) J  (0  ,)G  (0  .,0  )  A0 

u  xm  xi  xm  xi  xx  xi  z  x 


(2 . 81d) 


l  =-oo 

In  practice,  the  integrations  in  (2.80)  will  be  done  numerically, 
so  (2.81)  with  finite  limits  for  the  summations,  are  more 
representative  of  the  actual  equations  that  need  to  be  computed. 

The  left  hand  sides  of  (2.79a)  and  (2.79b)  can  be  shown  to  be  zero 
through  the  application  of  Parseval's  theorem,  which  states  [51): 

Given  two  functions  f(t)  and  g(t)  that  are  Fourier  transformable 
with  transforms  represented  by 
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G  (w)  and  f  (t) 


(2.82a) 


g(t)f(t)  dt 


G  (to)  F  (w)  dw 


(2.82b) 


where  denotes  the  complex  conjugate.  Thus  the  integration  of  the 
left  hand  side  of  (2.79a)  is 


J  O  )E  (|3  ,0  )  d|3  =  2n  J  (x)E  (x,/3  )dx 

zn  x  z2  x  z  x  J  zn  z2  z 

-oo  J  -00 


(2.83) 


Note  that  the  current  density  exists  only  on  the  strip  and  the 
tangential  electric  fields  are  zero  there  because  it  is  a  perfect 
electric  conductor.  The  tangential  electric  fields,  on  the  other  hand, 
are  non-zero  only  where  there  is  not  a  conductor,  and  therefore  where 
there  is  no  current  density.  Thus  the  product  of  the  current  density 
and  electric  field  is  always  zero  since  neither  is  non-zero  where  the 
other  exists.  Therefore,  the  integrand  is  zero  for  all  x  and  so  the 
integration  in  the  spatial  domain  is  also  zero.  Since  the  integration 
in  the  spatial  domain  is  zero,  through  Parseval's  theorem,  the 
integration  in  the  spectral  domain  is  also  zero.  In  a  similar  fashion, 
the  left  side  of  (2.79b)  can  also  be  shown  to  be  zero.  Replacing  the 
infinite  limits  of  the  summations  with  finite  ones,  (2.79a)  and  (2.79b) 


become 


N  N  M 

Fa  FaK  +  F  b  K  =0 
,  n  .  n  nn  m  mn 

n=l  n=l  m=l 


M  N  M 

Fb  FaK  +  F  b  K  =0 
,  m  ,  n  nm  ,  m  mm 

m=l  n=l  m=l 


(2 . 84a) 


(2.84b) 


where  N  and  M  are  the  total  number  of  current  densities  used  in  the 


expansion . 
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These  equations  can  be  solved  numerically  for  /3  as  a  function  of 

z 

frequency  and/or  structural  parameters  by  using  a  nonlinear  equation 
solver  to  find  the  zero  of  the  determinant  of  the  resulting  matrix.  To 
obtain  an  exact  solution,  the  number  of  basis  functions  used  must  be 
infinite.  However,  if  the  basis  functions  are  chosen  properly,  so  that 
they  closely  model  the  actual  current  distributions,  then  few  terms  are 
necessary  to  obtain  very  good  results.  For  frequencies  well  into  the 
dispersive  region  of  the  structure,  only  a  two  by  two  matrix  (i.e.  m=1, 
n=1)  is  necessary  to  get  very  good  results. 


2.3.7  Slot  Lines  and  Coplanar  Waveguides 

For  slot  line  and  CPW  type  structures,  the  electric  fields  exist 
over  a  smaller  portion  of  the  center  conductor  interface  than  the 
current  densities.  Therefore  it  is  usually  more  convenient  to  specify 
the  electric  fields  in  the  spaces  as  opposed  to  the  currents  on  the 
conductors.  To  solve  these  types  of  problems,  it  is  advantageous  to 
reformulate  (2.74a)  and  (2.74b)  so  that  the  currents  are  expressed  in 
terms  of  the  electric  fields. 

The  analysis  begins  by  simplifying  (2.74a-b)for  the  fields  on  the 

conductor  interface.  Both  sides  of  these  equations  are  multiplied  by 
-  j  2  2 

—=4— (8  +B  )  and  they  can  then  be  written  as 
u>Ho  x  z 


^-E  tf2+l32)  =  J  |py2Z™-,32ZTE 

z2  X  z  z  L  z  o  X 


]/\/  r  (  -  2~TM  ~TE  I 
+Jx[«{0o  Z  +Z 


-  ■)  'v  2  -  2 

— E  O  +£  )  =  J 
U)U  X2  X  z  z 

0 


(2.85a) 

(2.85b) 
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These  equations  are  now  in  the  form  of  a  matrix  equation,  i.e. 


^oV) 

X  z 


E 

a  b 

"  J 

Z2 

= 

X 

E 
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J 
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L  J 

L  ZJ 
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ig'V\zTEl 

|h0 


- 


o  =  ey2z™-/3VE 

pxp0  pz 

The  matrix  equation  in  (2.86)  can  be  solved  for  and  by 
multiplying  both  sides  by  the  inverse  of  the  square  matrix;  i.e 


J 

*-  zJ 


=  y  yv  > 

W/i  X  z 
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a  b 
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E 
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L  x3J 
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e  f 
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Z3 

f  9 

E 

L  X3J 

where 


e  =  c/ (ca-b  ) 

f  =  -b/(ca-b2) 

2 

g  =  a/ (ca-b  ) 

Begin  by  finding  the  denominator  of  the  matrix  elements 
(ca-b2) 


*2 - 2  TM 

_2  TE" 

*2-2  TM 

_2  TE- 

P  P  z 

-  p  z 
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X  0 

z 

(ca-b2)  =  pV^z™  yyyyyyyy” 
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Z  '  X  z 
2 
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(2.86) 

(2.87a) 

(2.87b) 

(2.87c) 

r 

(2.88) 

(2.89a) 

(2.89b) 

(2.89c) 

(2.90) 

(2.91) 

(2.94) 


(2.95) 
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Y™  =  l/z™  (2.96a) 

YTE  =  l/zTE  (2.96b) 

Using  these  substitutions  with  the  equation  for  the  denominator  in 


yte  =  i/zte 


(2.95),  the  elements  of  the  inverted  matrix  are; 
e  =  c/(ca-b2)  =  [eV Y™  -  62YTE1  \b2+B2] 
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„2  TE 
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(2.97a) 

(2.97b) 

(2.97c) 


Substituting  back  into  the  matrix  and  writing  out  the  equations 


j  =  Teh  +e  h  ] 

z  wp  *-  z2  zz  x2  xzJ 


j  =  Teh  +e  h 

x  we  L  z2  zx  x2  xxJ 
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(2.98a) 

(2.98b) 


where 
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(2.99a) 


(2.99b) 


(2.99c) 


Equations  (2.98a)  and  (2.98b)  can  then  be  solved  in  the  same 

manner  as  (2.76a)  and  (2.76b)  for  ft  by  expanding  the  electric  fields 

z 

instead  of  the  currents.  The  expansion  functions  for  the  electric 

field  are  chosen  to  be  the  same  as  the  current  expansion  functions, 

except  that  the  J  functions  are  used  for  E  expansions  and  the  J 

zn  xm  xm 


functions  are  used  for  the  E  expansions. 

zn 
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2.3.8  Green's  Function  for  Multilayer  Structures 

With  the  increasing  emphasis  on  MMIC  technology,  the  need  to 
compute  microstrip  parameters  for  multilayer  structures  has  become  more 
important.  The  introduction  of  additional  substrate  layers,  even 
relatively  thin  ones,  has  a  strong  effect  on  the  effective  dielectric 
constant  of  the  structure  [52] .  In  general,  the  addition  of 
superstrates  does  not  significantly  affect  the  parameters  of  planar 
circuits.  However,  in  certain  cases,  especially  when  the  circuit  is 
shielded,  these  layers  also  can  have  a  large  impact  on  the  propagation 
characteristics  of  the  structure.  Although  the  effect  of  multiple 
layers  is  very  important,  almost  all  of  the  approximate  formulations 
have  been  done  for  a  single  dielectric  substrate,  and  they  cannot  be 
modified  to  handle  more  general  cases.  However,  because  of  their  exact 
formulation,  most  full  wave  solutions  can  treat  the  multilayer  problem. 

One  drawback  with  using  full  wave  solutions,  such  as  the  SDA,  to 

solve  multilayer  problems  is  that  they  normally  require  a  different 

computer  program  for  each  structure  that  would  potentially  be  studied. 

In  the  case  of  the  SDA,  a  new  Green's  function  must  be  derived  for  each 

possible  configuration  by  solving  the  appropriate  boundary  value 

problem.  As  the  number  of  layers  increases,  the  number  of  unknown 

constants  and  equations  in  the  boundary  value  problem  quadruples.  Thus 

to  obtain  the  Green's  function  of  a  structure  with  six  dielectric 

X  X  z  z 

layers  using  a  combination  of  TE  ,  TM  ,  TE  ,  or  TM  modes  would  require 

simultaneously  solving  twenty  four  equations  for  twenty  four  unknowns. 

y  y  , 

If  TE  and  TM  modes  are  used,  the  problem  simplifies  slightly  to  two 
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separate  sets  of  twelve  equations  and  twelve  unknowns.  Thus  the 

development  of  an  easily  computed,  generalized  Green's  function  for 

planar  printed  circuit  lines  will  greatly  simplify  both  the  theoretical 

development  as  well  as  the  computational  process. 

In  this  section,  a  general  Green's  function  is  presented  that  is 

easily  computed  through  a  recursive  formulation.  This  Green's  function 

is  applicable  to  structures  with  any  number  of  layers  above  or  below 

the  center  conductor  interface.  Since  the  information  about  the  y 

variations  in  the  structure  are  contained  in  the  modal  input 

impedances,  they  are  the  only  part  of  the  Green's  function  that  must  be 

recalculated  for  multi-layer  structures.  The  input  impedance  for 
.  V  V 

either  the  TE  or  TM  mode  can  be  thought  of  as  a  parallel  combination 
of  the  admittances  seen  above  and  below  the  center  conductor  interface. 


i.e. 


■  ( i) 


~(i)  ~(i) 
Y  +Y 
LN  UM 


(2.100) 


where  i  is  either  te  or  tm,  l  and  u  indicate  the  lower  and  upper  layers 

respectively,  and  n  and  m  are  the  total  number  of  planar  layers  below 

and  above  the  interface,  as  in  Fig.  2.2.  The  admittances  Y<x)  are 

(j) 

determined  using  a  recursive  formulation  for  either  the  upper  or  lower 
layers,  beginning  with  j  =  Ll  or  ui  through  j  =  ln  or  um,  using 


yTM 

(j) 


~TE 

(J) 


~TM  ~TM  ,  2  .2 

Y  Y  +a  /ll 

s  ( j )  (j-l)  y  ( j)  'r  ( j) 

~TM  +  ~TM 
S  ( j)  (j-l) 

~TE  ~TE  2  .2 

Y  Y  +E  /a 

s(j)  (j-l)  r  ( j)  y  ( j) 

~TE  +  ~TE 

S(j)  (j-l) 


(2.101) 


(2.102) 
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yTE 
S  ( j) 


*TM 
3  { j) 


a  .cothTa  h 

y ( j )  L  y < j)  (j)J 

11 

r  <  j ) 

c  .cothTa  h 

r(j)  L  vMl 


y(j)  (D- 


a 


(2.103) 


(2.104) 


y  ( j) 


r(i> 


Yj~J  can  be  thought  of  as  the  input  admittance  seen  looking  outward 
through  the  jth  layer  (away  from  the  conductor  interface)  and  is 

the  self  admittance  of  the  layer  for  the  particular  mode  configuration. 

For  j  =  in  or  LI,  then  the  (j-l)  layer,  uo  or  lo,  is  a  perfect 
electric  conductor  (i.e.,  a  ground  plane  or  cover  sheet),  and 

/n/  ( i )  , 

Y  — >  oo,  for  k  =  u  or  l,  since  a  perfect  conductor  has  zero 

( k)  0  ^ 


resistance.  Taking  the  limit  of  either  (2.101)  or  (2.102b)  as  Y 

~(i)  rsj  (  i  ) 

goes  to  infinity,  then  Y  =  Y 

( k)  1  S (k )  1 

If  the  (k )  1th  layer  is  dielectric  of  infinite  height  (i.e., a 


(i) 
(k)  0 


structure  with  no  ground  plane  or  cover  sheet) ,  then  h 


(k)l 


oo  and 


therefore  coth (a  h  )  — >  1,  respectively.  The  conditions  for 
y(k)i  (k) l 

~(i)  ~(i) 

Y  and  Y  are  the  same  as  above,  as  though  there  is  a  conductor 
(k) 0  ( k) 1 

at  an  infinite  distance.  This  is  a  valid  assumption  since  the  outward 
radiation  condition  specifies  that  E  — »  0  as  y  — >  ±oo.  This  particular 
configuration,  no  ground  plane  and/or  cover  sheet,  results  in 

(2.105) 

(2.106) 

Note  that  these  equations  are  the  same  for  layers  both  above  and 
below  the  conductor  interface  because  the  uniqueness  theorem  states 
that  solution  of  the  boundary  value  problem  is  invariant  to  the  method 
of  solution  (if  done  correctly) .  Certainly  it  is  not  possible  for  the 
structure  to  'know'  which  layers  are  'above'  or  'below' .  This  symmetry 


~te  ~TE  . 

Y  =  Y  =  a  /u 

(k)  1  s  (k)  1  y  (k)  1  (k)  1 

~TM  ~TM  , 

Y  =  Y  =  c  /a 

( k)  1  S  (k)  1  r(k)l  y  (k)  1 
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in  the  equations  is  valid  only  for  strips  of  vanishingly  small 
thickness  and  becomes  invalid  when  the  strip  thickness  approaches  the 
height  of  the  dielectric  layer  that  contains  it. 

Thus  the  Spectral  Domain  Approach  is  a  simple,  fast  and  efficient 

method  for  determining  the  propagation  constants  of  planar  transmission 

lines.  Starting  by  Fourier  transforming  of  all  the  field  quantities, 

....  y  y 

the  derivation  of  the  Green's  function  is  simplified.  TE  and  TM 
modes  are  then  used  because  they  can  independently  satisfy  the  boundary 
condition  in  the  spectral  domain.  The  total  electric  fields  at  the 
interface  can  then  be  obtained  by  relating  the  modal  currents  to  the 
total  current  densities.  After  the  total  electric  fields  are  found, 
the  total  currents  are  expanded  into  a  complete  set  of  basis  functions 
so  that  the  problem  can  be  solved  using  Galerkin's  method.  Finally  the 
equations  are  solved  numerically  to  obtain  the  propagation  constants  of 


the  structure  under  consideration. 


CHAPTER  3 

TRANSIENT  SIGNAL  ANALYSIS 


While  many  applications  of  printed  circuit  transmission  lines 
require  only  single  frequency  operation,  there  also  is  a  demand  for 
circuits  that  use  finite  time  domain  signals.  This  chapter  presents 
the  formulations  necessary  to  analyze  the  distortion  of  transient 
signals  on  isolated,  as  well  as  coupled  transmission  lines.  Since 
planar  printed  circuits  are  inherently  dispersive,  i.e.  the  phase 
velocity  is  a  function  of  frequency,  the  analysis  of  transient  signal 
distortion  on  these  lines  begins  with  a  frequency  domain  approach.  In 
order  to  consider  the  coupling  between  adjacent  lines  in  a  rigorous 
manner,  the  even/odd  mode  formulation  is  then  presented.  Finally,  this 
formulation  is  used  to  derive  the  appropriate  equations  for  the 
propagation  of  finite  time  domain  signals  on  dispersive  and  lossy 
symmetric  coupled  transmission  lines. 

3.1  The  Frequency  Domain  Approach 

The  analysis  of  the  time  dependent  response  of  a  transient  signal 
may  be  formulated  in  either  the  time  domain  or  the  frequency  domain. 

For  a  non-dispersive  transmission  line,  it  is  easiest  to  analyze  signal 
propagation  in  the  time  domain.  On  the  other  hand,  if  the 
transmission  line  is  dispersive,  then  the  analysis  in  the  time  domain, 
while  possible  [14],  is  much  more  difficult  and  time  consuming.  The 
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simplest  and  most  efficient  method  for  examining  pulse  propagation  on 
these  types  of  structures  is  to  consider  the  pulse  in  the  frequency 
domain. 


The  analysis  of  the  system  response  begins  by  defining  a  Fourier 
transform  pair  that  relates  the  time  domain  response  to  the  frequency 
domain  spectrum.  Using  a  time  variable  t  and  the  corresponding 


transform  variable  w,  the  Fourier  transform  pair  is  defined  as 


V(w,z)  =  ?jv(t,z)  1  =  f  v(t,z)  e  dt 

v  '  -T/2 

v(t,z)  =  ?  ^VtWjZ)!  =  J  V(w,z)  e ^ 


(3.1) 


(3.2) 


where  T  is  the  total  time  over  which  the  pulse  exists.  Since  the 


signal  exists  only  over  a  finite  time  period,  its  spectrum  V  (w, z)  is 


continuous  and  exists  from  -co  <  w  <  oo. 


The  time  domain  response  of  a  linear  system  to  an  input  signal  is 
equal  to  the  convolution  integral  of  the  impulse  response  of  the  system 
and  the  input  signal.  If  the  input  signal  is  v(t,z)  and  the  impulse 
response  of  the  system  is  h(t,z),  then  the  system  response,  in  time,  is 
given  by 

00 

y(t,z)  =  T  h(r,z)  v(t-x,z)  dr  (3.3) 

**  -CO 

Since  the  convolution  integral  of  two  functions  in  the  time  domain  is 
equivalent  to  multiplication  of  their  Fourier  transforms  in  the 
frequency  domain,  (3.3)  can  also  be  written  as 

Y(u,  z)  =  H(w)  V(w,  z)  (3.4) 

where  Y(w, z),  H(w, z),  and  V(w, z)  are  the  Fourier  transforms  of  y(t,z), 
h(t,z),  and  v(t,z),  respectively.  To  find  the  system  response  in  time. 
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the  inverse  transform  of  (3.2)  is  applied  to  both  sides  of  (3.4)  giving 

00 

y  (t,  z)  =  ^  J  V(w,z)H(w)  e3^  du  (3.5) 

J  -00 

The  impulse  response  of  a  matched  transmission  line  is  simply  a 

phase  advance  equal  to  the  product  of  the  complex  propagation  constant 

and  the  distance  traveled,  i.e. 

-y  «a»  z  -[a  (w)  +  j/3  (w)  ]  z 

H(w,z)  =  e  Z  =  e  2  Z  (3.6) 

where  y  (a>)  is  the  complex  propagation  constant,  a  (u>)  is  the 
z  z 

attenuation  constant,  and  ft  (w)  is  the  phase  constant.  Using  (3.6)  for 

z 

H(w, z)  in  (3.5),  the  response  on  a  transmission  line  at  a  time  t  and  a 

position  zq  to  a  signal  v(t,z)  can  be  written  as 

oo^  -[a  (u)+j ft  (u>)  )  z  jut 

y(t,zo)  =  —  V(u,z)  e  Z  Z  e  du  (3.7) 

J  -00 

or 

roo  -a  (u)  z  j[wt-/3  (M)z  l 

i  r  ~  z  o  z  o 

y(t,zo)  =  —  V(u,z)  e  e  du  (3.8) 

-CO 

The  exponential  containing  a  in  the  integrand  causes  the  degradation 

z 

of  the  signal  due  to  attenuation  whereas  the  exponential  with  ft is 

responsible  for  the  phase  advance. 

Dispersion  distortion  happens  when  ft  (u)  is  a  nonlinear  function 

z 

of  frequency,  i.e.  the  phase  velocity  is  not  constant  with  varying 
frequencies.  The  effect  of  dispersion  distortion  can  be  visualized  by 
considering  the  pulse  as  a  summation  of  many  single  frequency  signals 
with  different  amplitudes,  i.e.  a  Fourier  series  representation.  At 
the  beginning  of  the  transmission  line,  all  the  signals  line  up 
properly  and  add  together  to  recreate  the  input  signal.  However,  as 
the  pulse  travels  down  the  line,  each  frequency  component  travels  at  a 
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slightly  different  speed,  since  the  phase  velocity  is  not  the  same  for 

each  frequency.  Thus,  each  frequency  component  of  the  pulse  arrives 

slightly  shifted  in  time  with  respect  to  the  other  components.  Since 

the  frequency  components  do  not  add  together  as  they  did  originally, 

the  pulse  shape  is  distorted.  Although  dispersion  distortion  does  not 

change  the  amplitudes  of  the  the  individual  frequency  components,  the 

amplitude  of  the  distorted  pulse  will  normally  be  less  than  the 

amplitude  of  the  original  pulse.  In  addition,  the  pulse  spreads  out  in 

time  since  some  frequency  components  arrive  earlier  and  some  later  than 

if  they  had  all  traveled  the  same  speed. 

If  a  (w)  is  a  constant  with  respect  to  frequency,  then  attenuation 
z 

will  only  reduce  the  amplitude  of  the  signal  and  not  distort  it.  On 

the  other  hand,  when  a  (W)  varies  with  frequency,  then  each  frequency 

z 

component  is  reduced  by  a  different  amount  and  it  changes  the 
reconstruction  of  the  signal.  However,  unlike  dispersion  distortion, 
attenuation  in  general  will  not  cause  widening  of  the  pulse,  only  a 
reduction  in  amplitude. 

3.2  Even/Odd  Mode  Formulation 

In  1956,  J.  Reed  and  H.  Wheeler  introduced  a  method  for  obtaining 
the  response  of  a  symmetric  four-port  network  called  the  even/odd  mode 
approach  [25]  .  The  even/odd  mode  approach,  also  known  as  the  symmetry 
approach,  simplifies  the  analysis  of  a  symmetric  four-port  network  by 
splitting  it  up  into  two  simpler,  two-port  networks.  Although  this 
method  of  analysis  applies  only  to  symmetric  networks,  asymmetric 
networks  can  be  analyzed  in  a  similar  fashion  [26] . 


48 


The  even/odd  mode  analysis  begins  by  assuming  that  a  signal  of 
unit  amplitude  is  applied  to  port  one  as  in  Fig.  3.1,  and  the  outputs 
are  measured  at  ports  two  and  three.  It  is  also  assumed  that  ports  two 
and  three  are  perfectly  matched,  so  that  there  are  no  reflections.  The 
input  signal  is  split  up  into  two  different  signals,  an  even  mode  and 
an  odd  mode  as  in  Fig.  3.1.  The  even  and  odd  mode  signals  are  chosen 
such  that  they  add  together  to  produce  a  unit  amplitude  signal  at  port 
one  and  no  signal  at  port  four,  recreating  the  original  signal.  The 
even  mode  consists  of  two  signals  of  one-half  amplitude  that  are  in 
phase  and  so  they  create  a  voltage  maximum  at  every  point  along  the 
line  of  symmetry.  This  is  equivalent  to  an  open  circuit  or  a  perfect 
magnetic  conductor  separating  the  two  lines.  The  odd  mode,  on  the 
other  hand,  has  two  signals  of  one-half  amplitude  that  are  180  degrees 
out  of  phase,  and  it  produces  a  voltage  minimum  at  every  point  along 
the  line  of  symmetry.  This  is  equivalent  to  a  short  circuit  or  having 
the  lines  separated  by  a  perfect  electric  conductor.  In  general,  the 
even  and  the  odd  modes  will  have  different  phase  velocities, 
attenuation  constants,  and  characteristic  impedances. 

Coupling  distortion  can  best  be  understood  by  considering  the 
response  on  both  lines  to  be  a  linear  combination  of  four  pulses,  two 
on  each  line,  as  in  Fig.  3.2.  In  general,  the  in  phase  pair  (even)  and 
the  out  of  phase  pair  (odd) ,  will  travel  down  the  line  at  two  different 
speeds,  due  to  the  differences  in  the  phase  velocities  of  the  even  and 
odd  modes.  The  even  and  odd  modes  pairs  of  pulses  add  constructively 
on  the  signal  line  and  destructively  on  the  sense  line.  To  isolate  the 
effects  of  even/odd  mode  distortion,  a  hypothetical  lossless  and 
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Fig.  3.2.  Coupling  distortion  on  lossless  and 
dispersionless  transmission  lines. 
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dispersionless  symmetric  two  conductor  transmission  line  is  considered. 
A  rectangular  pulse  is  used  for  simplicity,  and  it  is  assumed  that  the 
odd  mode  phase  velocity  is  higher  than  that  of  the  even  mode. 

When  the  signal  first  starts  out,  the  even  and  odd  pairs  have  not 
separated  very  much,  and  almost  completely  cancel  each  other  on  the 
sense  line  and  reproduce  the  signal  fairly  accurately  on  the  signal 
line,  as  shown  in  Fig.  3.2(b).  As  the  signal  travels  farther,  Fig. 
3.2(c),  the  even  and  odd  pairs  begin  to  separate.  Now  they  do  not 
cancel  completely  on  the  sense  line  or  reproduce  the  signal  very  well 
on  the  signal  line;  the  result  being  distortion  and  crosstalk.  The 
signal  line  response  has  also  been  widened  and  the  leading  and  trailing 
edges  have  dropped  off  in  amplitude.  After  the  signal  has  traveled  a 
very  long  distance,  Fig.  3.2(d),  the  even  and  odd  mode  pairs  separate 
completely,  so  that  there  are  two  pulses  of  one-half  amplitude  that  are 
in  phase  on  the  signal  line  and  two  pulses  of  one-half  amplitude  that 
are  180  degrees  out  of  phase  on  the  sense  line. 

In  Fig.  3.2(b) -(d),  the  leading  response  of  the  sense  line  is 
negative,  because  it  was  assumed  that  the  odd  mode  had  a  higher  phase 
velocity  than  the  even  mode.  If  this  is  not  the  case,  then  the  leading 
response  on  the  sense  line  is  positive,  because  the  even  mode  will 
arrive  before  the  odd  mode.  In  addition,  as  the  distance  increases, 
the  amplitude  of  the  response  on  both  lines  will  tend  to  approach 
one-half  of  the  value  of  the  response  on  a  single,  isolated  line. 

Since  pulse  spreading  is  due  to  the  difference  in  the  phase  velocities, 
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the  total  pulse  spread,  in  time,  due  to  only  even/odd  distortion  at  a 
distance  z,  can  be  written  as 


tg  =  |-^ - fH  (3.9) 

where  v  and  v  are  the  even  and  odd  mode  phase  velocities, 
e  o 

respectively.  Since  the  phase  velocities  are  related  to  the  effective 
dielectric  constants,  (3.9)  can  also  be  written  as 

(3.10) 

where  c  and  c  are  the  effective  dielectric  constants  for  the  even 


t  =  - 

Vc  -  Vc 

s  c 

re  ro 

re  ro 

and  odd  modes,  respectively,  and  c  is  the  speed  of  light  in  free  space 
Note  that  distortion  due  to  coupling  may  be  present  even  if  the  lines 
are  dispersive  and  lossy  as  well,  since  each  of  the  three  distortion 
mechanisms,  dispersion,  losses,  and  coupling,  are  independent  of  each 
other.  If  the  transmission  line  is  lossy  and  dispersive,  then  the 
analysis  is  the  same,  except  that  the  even  and  odd  mode  pairs  suffer 
distortion  due  to  dispersion  and  losses  as  they  travel  down  the  line, 
increasing  the  pulse  spread  and  reducing  further  the  amplitude  of  the 
signal . 


3.3  Symmetric  Coupled  Transmission  Lines 

In  this  section,  the  even/odd  mode  formulation  is  applied  with  the 
frequency  domain  approach  to  determine  the  response  on  symmetric 
coupled  lines.  To  compute  the  response  on  both  lines  due  to  an  input 
signal  on  only  one  line,  the  input  signal  is  split  into  the  even  and 
odd  mode  components.  The  response  on  both  lines  due  to  each  mode  is 
then  obtained  as  a  function  of  time  or  distance  for  a  specific  distance 
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or  time  respectively.  These  responses  are  then  added  together  to 
obtain  the  total  response  on  both  lines  at  a  particular  time  or 
distance . 

Using  this  approach,  the  response  on  both  lines  due  to  a  signal  on 
the  first  line,  the  signal  line,  with  no  excitation  on  the  second  line, 
the  sense  line,  can  be  written  using  the  superposition  of  the  even  and 
odd  mode  responses  as 

1 

(3.11a) 


v  (t,z)  =  i  v  (t,z)  +  v  ( t ,  z )  I 

1  2  L  e  °  J 

V2(t,z)  =  \  [ve(t,z)  "  Vo(t'z)] 


(3.11b) 


where  v  (t,z)  and  v  (t,z)  are  the  responses  of  the  even  and  odd  modes, 
e  o 

respectively,  to  the  input  signal  and  v^(t,z)  and  v2(t,z)  are  the 
voltages  on  lines  one  and  two  at  a  time  t  and  position  z.  Using  (3.8), 
the  even  and  odd  mode  responses  can  be  written  as 


1 

2n  . 

00 

V(U,  z) 

-y  (w)z 
ze 
e 

jut 

e 

du 

(3.12a) 

-00 

i 

r  00 

-y  (w) z 
zo 
e 

jut 

1 

“  2n  . 

V  (u>,  z) 

e 

du 

(3.12a) 

-00 

where  y  (u)  and  y  (u)  are  the  frequency  dependent  complex  propagation 
ze  zo 

constants  for  the  even  and  odd  modes  respectively.  Substituting  these 


two  equations  into  (3.11a)  and  (3.11b)  yields 

co 


Vl(t,Z)  =  4i  [  ^(U'Z) 

J  -00 

00 

Vt,z)  =  4i  J  ^(w'z) 


-T  z  jwt  -r  z  jut 
ze  zo 

e  e  +  e  e 


-y  z  jut  -y  z  jut 
ze  zo 

e  e  -  e  e 


dw 


dw 


(3.13a) 


(3.13a) 
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jwt  -z[y  +2T  ]/2 

Z  6  Z  O 

Factoring  out  e  e  in  (3.13a)  and  (3.13b)  leads  to 


'l(t'Z)=2iI  ^(W'Z) 

**  _  rv-i 


oo  jwt  -zfy  +y  ] /2 

ze  zo 

e  e  x 


1 

2 


z[y  -r  1/2  -z[y  -y  ] /2 
zo  ze  zo  ze 

e  +  e 


du> 


(3.14a) 


V2(t'Z>=2^J  ^(U'Z) 


00  jwt  -z[y  +y  ] /2 

ze  zo 

e  e  x 


1 

2 


]/2 


ze 


ze 


]  /  2 


dw 


(3.14b) 


These  equations  can  be  simplified  further  by  recognizing  that  the  terms 
in  the  parenthesis  in  (3.14a)  and  (3.14b)  are  the  hyperbolic  cosine  and 
sine  functions,  respectively.  Using  this  identity,  (3.14a)  and  (3.14b) 
can  be  expressed  as 


vi (t, z) 
V2(t,Z> 


sj. 

Hj 


jut-zt?  +y  ]/2 
ze  zo 

V(w,  z)  e  cosh 


jwt-z[y  ] /2 

~  ze  zo 

V(w, z)  e  sinh 


dw  (3.15a) 


dw  (3.15a) 


Equations  (3.15a)  and  (3.15b)  may  be  used  to  compute  the  pulse 
distortion  on  coupled,  lossy,  and  dispersive  transmission  lines  by 
evaluating  the  integrals  either  in  closed  form  or  numerically. 

However,  in  their  present  form  these  formulas  do  not  give  much  insight 
into  the  mechanisms  of  distortion  due  to  coupling,  losses,  or 
dispersion.  In  order  to  gain  an  understanding  of  these  coupling 
mechanisms,  equations  (3.15a-b)  need  to  be  simplified  to  separate  out 
each  of  the  three  mechanisms. 

The  first  step  is  to  expand  each  term  in  the  integrands  into  real 
and  imaginary  parts.  While  many  signals  of  interest  are  symmetric 
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about  some  specific  time  and  thus  have  a  completely  real  Fourier 
transform,  many  useful  signals  are  asymmetric  with  respect  to  time,  and 
therefore  their  Fourier  transforms  are  complex.  The  real  part  V(u,  z) 

t  •  /V 

is  designated  as  V  and  the  imaginary  part  as  V^.  In  addition,  to 
simplify  the  notation,  the  following  substitutions  are  made 


.  f  *zo-*ze  1 

Aa  =  Re^  j  = 

r  *Zo-*ze  I 

Ap  =  inj  —  j  = 


A0  =  Im 


a  =  Re 
av 


y  -y 

zo  ze 


y  +y 

zo  ze 


a  -a 
zo  ze 


0  -0 
zo  ze 


a  +a 
zo  ze 


(3.16a) 


(3.16b) 


(3.16c) 


0  =  Im 

av 


y  +y 

zo  ze 


0  +0 

zo  ze 


(3.1 6d) 


Using  these  replacements,  the  integrands  of  (3.15a)  and  (3.15b)  can  be 
rewritten  as 


I  =  V(w,z)  e 


jwt-z[y  +9c  ] /2 

ze  zo 


cosh  z 


y  -y 

zo  ze 


[fV  ^  I 

Vr  +  3  ViJ 


-aavz  j[«t-0avz]  r  I 

e  e  cosh  z (Aa  +  jA0) 


(3.17a) 


I  =  V(w,z)e 


jwt-zty  +  y  ]  / 2 
ze  zo 


sinh  z 


y  -y 

zo  ze 


[*N/  I 

Vr  +  3  ViJ 


-a  z  j [wt-0  z]  1 

e  e  sinh  z (Aa  +  jA0) 


(3.17b) 


Next  the  complex  exponential  in  (3.17a-b)  is  expanded  into  cosine  and 
sine  terms  giving 


I  =  e 
l 


-a  z 
av 


+  1  ’J  [ 

cosh  J^z  (Aa  + 


cos (wt -0  z)  +  i  si 
av 


n(wt-0  z) 
av  J 


jA0) 


(3.18a) 
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I  =  e 

2 


-a  z 
av  0 


*  1  ’J  [ 

sinh  j^z  (Aa  + 


cos  (wt-S  z)  +  j  sin  (wt-/3  z) 
av  av 


jA0) 


(3.18b) 


The  hyperbolic  sine  and  cosine  can  also  be  expanded  into  real  and 
imaginary  parts  with  the  following  formulas 


sinh(0+j0)  =  sinh (0)  cos  (0)  +  j  cosh  (0)  sin  (0) 


(3.19a) 


cosh(0+j0)  =  cosh (0) cos (0)  +  j  sinh (0) sin (0) 


(3.19b) 


Using  these  formulas  in  (3.18a-b)  leads  to 


I  =  e 
l 


K +  j  ~vJ  [ 


cos  (wt-/3  z)  +  j  sin(wt-p  z) 
av  av 


cosh(Aaz)  cos(Apz)  +  j  sinh(Aaz)  sin(A|3z) 


(3.20a) 


I  =  e 
2 


[vr  *  i  h]  [ 


cos (wt-6  z)  +  j  sin(wt-S  z)  x 
av  av 


sinh(Aaz)  cos  (A|3z)  +  j  cosh(Aaz)  sin(A|3z) 


(3.20b) 


Finally,  the  multiplications  in  (3.20a)  and  (3.20b)  are  carried  out  and 
the  results  are  separated  into  real  and  imaginary  parts.  The  real  part 


of  the  first  integrand  is 


Re(I  )  =  e 


-a  z 
av 


M 


cos (wt -B  z)  cos(Apz)  cosh(Aaz) 
av 


sin(wt-S  z)  sin(ASz)  sinh(Aaz) 
av 


-  V,  sin 
i 


(wt -B  z)  cos(A8z)  cosh(Aaz) 
av 


+  cos (wt-8  z)  sin(ASz)  sinh(Aaz) 
av 


(3.21a) 
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and  its  imaginary  part  is 


Imdj)  =  e 


-a  z 
av 


w 


sin(o)t -B  z)  cos(A/3z)  cosh(Aaz) 
av 


+  cos (wt z)  sin(Afiz)  sinh (Aaz) 
av 


<v 

V. 

1 


cos (wt -8  z)  cos(A8z)  cosh(Aaz) 
av 


-  sin(wt-8  z)  sin(A8z)  sinh(Aaz) 
av 


(3.21b) 


Similarly,  the  real  part  of  the  second  integrand  is  found  to  be 


Re  ( 1 2 )  =  e 


-a  z 
av 


w 


cos  (wt-/3  z)  cos(Apz)  sinh(Aaz) 
a  v 


-  sin(wt-/3  z)  sin(A8z)  cosh (Aaz) 
av 


/,  sii 
i 


V,  sin(wt-8  z)  cos(A8z)  sinh(Aaz) 
’  1  av 


+  cos  (cot-6  z)  sin(A/3z)  cosh(Aaz) 
av 


(3.21c) 


and  its  imaginary  part  is 


Im(I2>  =  e 


-a  z 
av 


sin(wt-/3  z) 

L  av 

cos  (A/3z) 

sinh (Aaz) 

+ 

cos(wt -8  z) 
av 

sin (Apz) 

cosh (Aaz) 

+  v. 

1 

cos (wt -8  z) 
av 

cos  (A/3z) 

sinh (Aaz) 

- 

sin(wt-8  z) 
av 

sin  (A/3z) 

cosh (Aaz) 

(3 . 21d) 


For  real  time  signals,  the  real  part  of  the  Fourier  transform,  V  , 

is  an  even  function  of  w,  and  the  imaginary  part,  V  ,  is  an  odd 

function  of  u)  [51]  .  In  addition,  while  f3  (w)  is  an  odd  function  of  0), 

z 


a  (G))  is  an  even  function  of  cj.  This  can  be  shown  by  considering  the 
z 


phase  propagation  equation  for  both  positive  and  negative  frequencies. 
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For  positive  frequencies,  the  phase  accumulation  of  a  propagating  wave 
is 

jwt-y  (w)  z  -a  (w)  z  j[wt-£  (a>)  z ] 

0(w)  =  e  Z  =eZ  e  Z  (3.22) 

and  for  negative  frequencies, 

-jwt-y  (-w)z  -a  (-w)z  - j [wt+/3  (-w)z] 

<f>  ( — w)  =  e  Z  =eZ  e  Z  (3.23) 

In  order  to  satisfy  the  outward  radiation  condition,  the  magnitude  of 

the  signal  cannot  increase  with  increasing  distance,  and  so  a  (-w)  must 

z 

be  equal  to  a  (w) .  Similarly,  in  order  for  phase  to  advance  as 
z 

distance  or  time  increases  then  /3  (-w)  must  be  equal  to  -[3  (w)  .  Thus 

z  z 

(w)  is  an  even  function  of  frequency  and  (w)  is  an  odd  function. 
Note  that  the  imaginary  part  of  the  integrands  in  (3.21b)  and 
(3.21d)  are  odd  functions  of  w  and  that  the  integrals  in  (3.15a)  and 
(3.15b)  are  over  symmetric  interval,  from  minus  infinity  to  plus 
infinity.  Since  the  integral  of  an  odd  function  over  a  symmetric 
interval  is  zero,  the  imaginary  components  of  the  integrand  do  not 
contribute  to  the  voltage  responses  and  can  be  ignored.  On  the  other 
hand,  the  real  part  of  the  integrands,  (3.21a)  and  (3.21c),  are  even 
functions  of  w,  and  so  the  integration  from  minus  infinity  to  plus 
infinity  will  be  simply  twice  the  integration  from  zero  to  either  plus 
or  minus  infinity.  This  is  the  anticipated  result,  since  if  real  time 
signal  is  input  to  a  real  system,  then  the  output  should  be  a  real 
function  as  well.  Eliminating  the  imaginary  parts,  the  two  integrands 


become 
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I  =  e 
l 


-a  z 
av 


w 


cos  (o)t-B  z)  cos(ASz)  cosh(Aaz) 
av 


sin(wt-0  z)  sin(A0z)  sinh(Actz) 

3  V 


S',  sii 

l 


V,  sin(wt-0  z)  cos (A0z)  cosh(Aaz) 
■  1  av 


+  cos (wt-0  z)  sin(A0z)  sinh(Aaz) 


]} 


(3.24a) 


I  =  e 
2 


-a  z 
av 


w 


cos (tot -3  z)  cos(Afiz)  sinh(Aaz) 
av 


-  sin  (tot -p  z)  sin(A/3z)  cosh(Aaz) 
av 


iV 

-  v,  si 

i 


in(wt-0  z)  cos(A0z)  sinh(Aaz) 
av 


+  cos (wt-0  z)  sin  (ASz)  cosh(Aaz) 
av 


]} 


(3.24b) 


av 


Equations  (3.24a)  and  (3.24b)  are  the  most  general  forms  of  the 
integrands  in  (3.15a)  and  (3.15b).  The  sine  and  cosine  terms  with  0 
represent  the  propagation  of  the  signals  at  an  average  phase  velocity, 
while  the  A0  sine  and  cosine  terms  account  for  the  distortion  of  the 
signal  due  to  even/odd  mode  coupling.  Similarly,  the  signals  can  be 
thought  of  as  having  and  average  attenuation  value  of  a^,  due  to  the 
real  exponential,  while  attenuation  coupling  between  the  lines  is  due 
to  the  hyperbolic  sine  and  cosine  terms  with  Act  arguments. 

Having  obtained  these  formulas,  some  special  cases  may  be 

considered.  If  the  attenuation  constants  are  the  same  for  each  mode, 

then  Aa  =  0  and  a  =  a  .  This  reduces  (3.21a-b)  to 
av  z 


-a  z 

I  =  e  Z  cos(A0z)  ■{ V  cos  (wt-0  z)  -  V,sin(wt-0  z) 

1  -  av  i  av 


I  = 

2 


-a  z 
2 

e  sin(A0z) 


{v 

•jv^sin 


(wt-0  z)  +  V,cos(wt-0  z) 
av  l  av 


(3.22a) 

(3.22b) 
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In  this  case,  the  signals  on  both  lines  are  attenuated  equally,  but 


there  is  coupling  due  to  A 8  and  dispersion  if  6  and  6  are  functions 

ze  zo 

of  frequency. 

If,  in  addition,  the  propagation  constants  for  the  even  and  odd 
modes  are  equal,  i.e.  Ap  =  0  and  |3  =  &zr  then  there  is  no  response  on 

the  sense  line,  and  the  signal  line  response  is  given  as 


I  =  e 
l 


-a  z 
z 


cos (wt-p  z)  -  V.sin(wt-p  z) 
Z  1  z 


(3.23) 


This  is  the  familiar  result  for  pulse  propagation  on  a  single, 
isolated,  lossy  transmission  line  [40],  [41],  [44],  [46]. 

However,  if  the  attenuation  constants  for  the  even  and  odd  modes 
are  different  but  the  propagation  constants  are  the  same,  then  the 


voltages  on  the  two  lines  can  be  expressed  as 
-a  z 

av  ,  ,  „ 

I  =  e  cosh(Aaz)  <V  cos (wt-8  z)  -  V,sin(a>t-6  z) 

l  I  r  av  i  av 


{v 

(v 


-a  z 

I  =  e  sinh(Aaz)  -(v  cos (wt -8  z)  -  V,sin(wt-/3  z) 

2  1  r  av  i  av 


(3.22a) 


(3.22b) 


Expanding  the  hyperbolic  sine  and  cosine  into  exponentials  and 
multiplying  the  result  with  the  exponential  in  front  of  (3.22a)  and 


(3.22b)  results  in 


■.-it 
‘  §[ 


-a  z  -a  z- 
ze  zo  | |~ 

e  +  e 


V  cos  (wt-/3  z)  -  V,  sin  (a)t -/3  z) 
r  av  l  av 


-a  z  -a  z 
ze  zo 

e  -  e 


cos  (wt-3  z)  -  V.  sin  (wt-p  z) 
av  i  av 


(3.22a) 


(3.22b) 


Since  the  exponentials  in  (3.22b)  are  not  equal,  then  I^  is 
non-zero  and  therefore  there  will  be  a  response  on  the  sense  line.  Thus 
even  if  the  phase  velocities  of  the  even  and  the  odd  mode  are  the  same, 
there  still  may  be  coupling  due  to  differences  in  the  attenuation 
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constants  for  the  two  modes.  The  response  on  the  sense  line  has  zero 


magnitude  at  zq  =  0  and  is  zero  again  as  zq  — >  oo.  The  maximum  response 
on  the  sense  line  occurs  when  the  derivative  with  respect  to  z  of  the 


exponentials  in  the  front  of  (3.22b)  vanishes,  i.e. 

0 


d_  r 

dz 


-a  z  -a  z 
zo  ze 

e  -  e 


Performing  the  indicated  derivative, 

-a  z  -a  z 

ze  zo  _ 

-a  e  +  a  e  =0 

ze  zo 

-a  z  -a  z 

ze  zo 

a  e  =  a  e 

ze  zo 


Both  sides  are  multiplied  by  exp (a  z) /a  ,  giving 

zo  ze 

-a  z+a  z  a 
ze  zo  zo 


(3.23) 

(3.24) 

(3.25) 


a 


(3.26) 


ze 


To  remove  the  exponential,  the  natural  logarithm  is  taken  of  both 
sides . 

-a  z+a  z  =  ln(a  /a  )  (3.27) 

ze  zo  zo  ze 

This  equation  can  now  be  easily  solved  for  z. 


z  = 


In (a  /a  ) 
zo  ze 


a  -a 


(3.28) 


zo  ze 

This  is  the  position  of  the  maximum  response  on  the  sense  line  due  to 
attenuation  coupling.  Note  that  if  the  attenuation  constants  are  a 
function  of  frequency,  which  is  usually  the  case,  then  location  of  the 
maximum  sense  line  response  is  different  for  each  frequency  component 
of  the  pulse.  Thus  it  may  not  be  possible  to  accurately  predict,  using 
(3.28)  the  actual  location  of  the  maximum  sense  line  response  due  to 
attenuation  coupling. 

Using  Fourier  transform  theory  and  the  even/odd  mode  approach,  the 
equations  for  signal  propagation  on  symmetric  coupled,  lossy,  and 
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dispersive  transmission  lines  were  derived.  The  Fourier  transform 
provides  a  simple  and  straightforward  way  to  represent  the  frequency 
dependent  characteristics  of  both  the  transmission  line  and  the  signal, 
allowing  quick  computation  of  the  time  domain  response  of  the  system. 
The  even/odd  mode  approach  reduces  the  symmetric  four-port  system  by 
representing  it  as  two  simple  two-port  networks.  Together,  these  two 
methods  were  used  to  derive  equations  that  relate  the  distortion 
mechanisms  of  dispersion,  even/odd  mode  coupling,  and  attenuation 
coupling. 


CHAPTER  4 
RESULTS 

This  chapter  uses  the  formulations  from  chapters  two  and  three  to 
compute  the  frequency  dependent  parameters  of  complex  microstrip 
structures  and  to  analyze  pulse  distortion  due  to  dispersion  and 
coupling  on  these  lines.  Initially,  a  very  simple  structure  is 
considered,  the  single  layer  open  symmetric  coupled  microstrip.  The 
effective  dielectric  constant  is  computed  versus  frequency  for 
different  substrate  materials  and  strip  spacings.  Using  this  data, 
pulse  distortion  on  these  structures  is  considered  for  different 
lengths  and  spacings.  Next,  the  generalized  Green's  function  is  used 
to  compute  the  effective  dielectric  constant  of  structures  with 
multiple  dielectric  layers.  Again  results  are  presented  for  different 
spacings,  materials,  and  heights  of  the  layers.  Finally  multilayer 
symmetric  coupled  microstrips  are  considered,  particularly  with  respect 
to  structures  that  exhibit  low  coupling  due  to  the  control  of  the 
electrical  properties  of  the  substrates  and  superstrates . 

4.1  Symmetric  Coupled  Microstrips 

Coupling  between  adjacent  transmission  lines  is  a  very  important 
consideration  in  circuit  design,  since  it  is  a  limiting  factor  in  the 
overall  size  and  operation  speed  of  the  circuit.  However,  due  to  the 
wideband  nature  of  finite  time  domain  signals,  no  single  parameter  or 
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rule  of  thumb  can  accurately  take  into  account  the  effects  of  coupling 
on  pulse  propagation,  especially  for  multilayer  structures.  For 
accurate  results,  the  analysis  needs  to  be  carried  out  in  a  more 
rigorous  manner,  considering  the  distance  the  signal  travels  and  the 
entire  spectrum  of  both  the  signal  and  the  impulse  response  of  the 
network . 

This  section  considers  open  symmetric  coupled  microstrip 

structures  with  a  single  dielectric  substrate.  While  this  is  not  a 

very  complex  structure,  it  is  very  useful  in  demonstrating  many  of  the 

principles  of  coupling  and  dispersion  distortion.  Four  different 

dielectric  substrates  are  used  throughout  this  chapter  as  typical 

examples  of  some  of  the  different  types  of  materials  that  are  used  in 

production:  1)  RT/duroid  5880,  =  2.2,  2)  beryllium  oxide,  =  6.8, 

3)  alumina,  e  =9.7,  and  4)  gallium  arsenide,  c  =12.2.  To 
r  r 

demonstrate  distortion,  a  Gaussian  pulse  is  used  which  has  time  and 

frequency  responses  given  by 

-In  (2)  (t/x) 2  _  | - - —  - (wt/2) 2/ln  (2) 

v  ( t )  =  Ae  < — )  V(w)  =  Ax  e  (4.1) 

where  A  is  the  amplitude  of  the  pulse,  t  is  the  voltage  half  width- 

half  amplitude  maximum,  and  the  pulse  is  centered  about  t  =  0. 

The  effective  dielectric  constant  of  the  even  mode,  e  ,  and  the 

re 

odd  mode,  e  ,  are  plotted  as  a  function  of  frequency  in  Figs.  4.1 
ro 

through  4.8  for  the  four  different  substrate  materials  with  six 
different  spacings .  The  even  mode  is  shown  in  Figs.  4.1,  4.3,  4.5,  and 
4.7,  while  the  odd  mode  is  graphed  in  Figs.  4.2,  4.4,  4.6,  and  4.8. 


The  smallest  spacing,  s  =  0.6  mm  is  very  close  and  would  normally  be 
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w  =  0.6mm,  h  =  0.6mm 


Frequency  in  GHz 

Fig.  4.2.  e  vs.  frequency  for  open  symmetric  coupled  microstrips  on 
RT/duroid  5880  substrate,  odd  mode. 
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Fig.  4.5.  e 
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e  ^  vs.  frequency  for  open  symmetric  coupled  microstrips  on 
alumina  substrate,  even  mode. 
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Frequency  in  GHz 

Fig.  4.6.  £  vs.  frequency  for  open  symmetric  coupled  microstrips  on 

alumina  substrate,  odd  mode. 
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Fig.  4.7.  e  vs.  frequency  for  open  symmetric  coupled  microstrips  on 
gallium  arsenide  substrate,  even  mode. 
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Fig.  4.8.  crefj  r  vs.  frequency  for  open  symmetric  coupled  microstrips  on 
gallium  arsenide  substrate,  odd  mode. 
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used  only  when  strong  coupling  between  the  lines  is  desired.  The 

values  for  s  =  5.0  mm,  both  even  and  odd  mode,  are  very  close  to  the 

values  for  s  — >  oo,  i.e.  a  single  isolated  microstrip,  and  so  e  “  e 

re  ro 

for  this  spacing.  In  microwave  circuit  design,  the  most  common  rule  of 

thumb  for  designing  uncoupled  lines  is  to  place  them  at  least  three 

center  conductor  widths  or  three  substrate  heights  apart.  For  the 

substrate  height  and  center  conductor  width  of  Figs.  4. 1-4. 8,  this 

spacing  is  represented  by  the  s  =  1.8  mm  graphs.  The  greatest 

difference  in  e  and  e  for  the  single  layer  configuration  occurs 
re  ro 

when  the  spacing  is  the  smallest.  As  the  spacing  increases,  e 

decreases  in  value,  e  increases,  and  both  asymptotically  approach  the 

value  of  c  ..  for  a  single  isolated  microstrip, 
reff 

As  can  be  seen  from  the  derivation  in  Chapter  3,  dispersion 
distortion  is  due  to  the  change  in  the  phase  velocity  with  frequency. 

In  particular,  dispersion  distortion  of  finite  time  domain  signals  is 
related  to  the  differences  in  the  inverse  of  the  phase  velocities  over 
the  band  of  frequencies  of  interest.  Since  the  inverse  of  the  phase 
velocity  is  proportional  to  the  square  root  of  the  effective  dielectric 


constant,  the  amount  of  distortion  due  to  dispersion  is  related  to 


l  Vc  ..(f)  -  v'c  ..(f) 
reff  l  reff  2 


(4.2) 


where  l  is  the  distance  traveled  down  the  line,  and  f^  and  f  are  the 


lower  and  upper  frequencies  that  define  the  band  of  interest.  Thus, 


the  amount  of  dispersion  that  is  acceptable  depends  on  three  factors; 

1)  how  e  varies  as  a  function  of  frequency,  2)  the  bandwidth  of  the 

reff 

system  and  signal,  and  3)  the  distance  the  pulse  must  travel.  The 
least  dispersive  structure  of  the  four  is  the  RT/duroid  5880  substrate, 
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whose  results  are  shown  in  Figs.  4.1  and  4.2,  since  its  effective 
dielectric  constant  changes  very  little,  only  about  0.23  from  0  to  125 
GHz.  The  beryllium  oxide.  Fig.  4.3  and  4.4,  and  alumina  substrates, 
Fig.  4.5  and  4.6,  are  a  little  more  dispersive  than  the  duroid  and  the 
gallium  arsenide  is  the  most  dispersive,  showing  a  change  in  of 

about  3.8  from  0  to  125  GHz. 

While  dispersion  is  an  important  consideration  in  the  design  of 

MMIC's,  the  effects  of  coupling  distortion  and  crosstalk  also  are  a 

major  concern.  Since  both  coupling  distortion  and  crosstalk  are  due  to 

the  difference  in  the  even  and  odd  mode  phase  velocities,  the  amount  of 

pulse  spread  due  to  even/odd  coupling  is  characteristic  of  the  amount 

of  pulse  distortion  due  to  coupling.  From  (3.24a-b)  the  amount  of 

pulse  spread  due  to  differences  in  c  and  e  can  be  written  as 

re  ro 


=  ir 

c 


t  (u)  =  -I  Ve  (wj~  -  Vc  (wj~ 


re 


ro 


(4.3) 


where  t  (w)  is  the  amount  of  pulse  spread  in  time  for  a  given 

frequency,  l  is  the  distance  traveled,  and  c  is  the  speed  of  light. 

Figs.  4.9  and  4.10  show  e  and  e  as  a  function  of  frequency  for  the 

re  ro 

four  different  substrates.  As  the  substrate  dielectric  constant 


increases,  the  separation  between  e  and  e  at  a  given  frequency 

increases  as  well.  Thus  the  gallium  arsenide  substrate  has  a  much 

tighter  coupling  between  the  lines  than  the  duroid  for  structures  with 

identical  dimensions.  For  the  single  layer  microstrip,  the  largest 

difference  in  c  and  c  occurs  at  low  frequencies,  in  the  quasi-TEM 
re  ro 


region,  although  this  is  not  the  case  for  all  structures,  as  will  be 


reff 
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shown  in  section  4.3.  As  frequency  increases,  e  ^  and  c  q  approach  the 

value  of  e  __  for  the  isolated  microstrip, 
ref  f 

This  behavior  is  characteristic  of  all  coupled  microstrip 
structures,  even  multilayer  ones,  because  as  frequency  increases,  the 
electrical  distance  separating  the  coupled  lines  increases,  decreasing 
the  coupling  and  enhancing  the  isolation  of  the  lines.  So  while 
dispersion  is  a  high  frequency  concern,  since  it  increases  as  the 
bandwidth  of  the  signal  increases,  even/odd  mode  coupling  is  normally  a 
low  frequency  phenomenon,  decreasing  as  frequency  increases.  Therefore, 
even  if  the  signal  is  narrowband  and  therefore  would  not  experience 
significant  dispersion  distortion,  even/odd  mode  distortion,  however, 
may  be  an  important  factor  in  the  design. 

Now  that  the  frequency  dependent  effective  dielectric  constant  has 

been  computed  for  the  structures,  pulse  distortion  can  be  considered 

using  the  Fourier  transform  approach.  Throughout  this  chapter  four 

different  designations  are  used  to  refer  to  the  time  domain  responses 

of  the  transmission  lines;  signal  line,  sense  line,  undistorted,  and 

isolated.  The  undistorted  designation  refers  to  pulses  that  are 

propagated  on  a  TEM,  lossless,  uncoupled  transmission  line  whose 

effective  dielectric  constant  is  invariant  with  respect  to  frequency 

and  it  is  equal  to  the  zero  frequency  value  of  the  structure  being 

considered,  i.e.  c  „{w)  =  e  „^(0)  for  all  w.  The  isolated  response 

reff  reff 

indicates  that  the  pulse  is  sent  down  a  dispersive,  lossless,  and 
uncoupled  transmission  line.  The  effective  dielectric  constant  for 
this  pulse  is  computed  using  the  SDA  with  just  a  single  center 
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conductor  whose  width  is  the  same  as  the  width  of  either  center 
conductor  of  the  symmetric  coupled  lines. 

When  the  signal  is  impressed  on  a  coupled  line,  there  are  two  time 
domain  responses  of  interest,  the  response  on  the  line  where  the  signal 
is  impressed,  called  the  signal  line,  and  the  response  on  the  line 
coupled  to  it,  referred  to  as  the  sense  line.  The  amount  of  distortion 
due  to  dispersion  is  shown  by  the  difference  between  the  undistorted 
and  the  isolated  response.  Distortion  due  to  even/odd  mode  coupling  is 
represented  by  the  difference  between  the  isolated  line  and  the  signal 
line  waveforms,  and  the  sense  line  response  gives  an  indication  of  the 
amount  of  crosstalk  present  in  the  structure. 

The  effects  of  dispersion  distortion  and  even/odd  mode  coupling 
distortion  on  pulse  propagation  are  shown  in  Figs.  4.11  through  4.18 
for  different  distances  and  substrate  materials.  A  Gaussian  pulse  with 
an  amplitude  of  5  and  a  t  equal  to  30  picoseconds  is  used  along  with  a 
spacing,  center  conductor  width,  and  substrate  height  of  0.6  mm.  Since 
the  pulse  travels  at  different  speeds  on  each  substrate,  the  distances 
were  scaled  appropriately  to  have  all  the  responses  arrive  at  similar 
times  and  to  more  clearly  illustrate  the  distortion  mechanisms.  Also, 
two  distances  were  chosen  for  each  pulse,  in  order  to  show  how 
increasing  distance  affects  both  coupling  and  dispersion  distortion. 

The  first  structure  considered,  shown  in  Figs.  4.11  and  4.12,  uses 
RT/duroid  as  the  substrate  material.  In  spite  of  the  large  distance 
that  the  input  signal  has  traveled,  it  shows  very  little  distortion  due 
to  dispersion.  Even  at  250  mm,  the  isolated  and  undistorted  responses 
agree  very  well.  Even/odd  mode  distortion,  on  the  other  hand,  has 
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Fig.  4.12.  Pulse  distortion  on  coupled  lines,  s  =  0.6mm,  1 
r  =  30ps,  RT/duroid  5880  substrate. 
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significantly  degraded  the  signal  line  response  at  just  125  mm,  and 
shows  a  large  response  on  the  sense  line  as  well.  At  250  mm  the  signal 
has  become  almost  unrecognizable  and  the  sense  line  response  has  risen 
to  almost  the  same  amplitude  as  the  signal  line. 

As  the  relative  dielectric  constant  of  the  substrate  increases, 

changing  from  2.2  to  6.8  and  9.7,  the  effects  of  dispersion  become  more 

pronounced  due  to  the  greater  change  in  the  phase  velocity  with 

frequency  for  these  structures.  The  beryllium  oxide  substrate.  Figs. 

4.13  and  4.14,  has  slightly  more  distortion  due  to  dispersion, 

particularly  at  1  =  150  mm,  and  even/odd  mode  distortion  affects  the 

pulse  even  more  than  the  duroid  structure  did.  At  150  mm,  the  signal 

line  response  begins  to  look  like  two  separate  in  phase  pulses  and  the 

sense  line  like  two  pulses,  180  degrees  out  of  phase.  Also,  like  all 

the  pulse  response  graphs  in  this  section,  the  leading  response  on  the 

sense  line  is  negative,  indicating  that  the  odd  mode  leads  the  even 

mode.  In  other  words,  the  odd  mode  phase  velocity  is  greater  than  the 

even  mode  and  hence  c  >  c  This  is  characteristic  of  all  single 

re  ro 

layer  microstrip  structures,  but  as  will  be  shown  in  section  4.2,  it  is 
not  necessarily  true  for  multilayer  structures.  The  structure  with 
alumina  substrate.  Figs.  4.15  and  4.16,  has  results  very  similar  to  the 
beryllium  oxide,  except  that  the  distances  have  decreased  from  75  and 
150  mm  to  62.5  and  125  mm  and  the  distortion  of  the  pulses  has 
increased  slightly. 

The  gallium  arsenide  structure,  as  anticipated,  has  the  greatest 
dispersion  and  coupling  distortion  of  the  four  substrate  materials.  In 
spite  of  having  scaled  the  distances  down  to  50  and  100  mm,  there  is 
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still  significant  distortion,  both  from  coupling  and  dispersion.  At 
50mm,  dispersion  has  shifted  the  pulse  in  time,  slowing  it  down 
slightly,  as  well  as  widening  it.  Coupling  had  reduced  the  signal  line 
to  70  percent  of  the  original  amplitude  and  has  a  sense  line  response 
that  is  almost  50  percent  of  the  undistorted  pulse.  By  100  mm,  both 
the  signal  and  sense  line  responses  have  separated  almost  completely 
into  the  even  and  odd  mode  pulses.  Unlike  the  other  substrates, 
dispersion  adds  significant  distortion  to  the  pulse,  reducing  the 
amplitude  by  ten  percent  and  significantly  increasing  the  pulse  width. 

The  effect  of  center  conductor  spacing  on  the  transient  response 
of  the  coupled  transmission  lines  is  as  important  as  the  distance 
considerations.  Figures  4.19  through  4.26  show  the  effects  of 
different  line  spacings  on  pulse  distortion  for  the  four  different 
substrates.  The  signal  line  responses  are  given  in  Figs.  4.19,  4.21, 
4.23,  and  4.25  with  the  corresponding  sense  line  values  in  Figs.  4.20, 
4.22,  4.24,  and  4.26.  The  isolated  response,  i.e.  s  =  oo,  is  not 
included,  but  is  closely  approximated  by  the  s  =  5  mm  response  on  the 
signal  line  and  by  a  zero  response  on  the  sense  line. 

As  before,  the  structure  with  RT/duroid  substrate,  with  results 
shown  in  Figs.  4.19  and  4.20,  is  not  significantly  affected  by  the 
dispersive  effects  of  the  line,  but  is  highly  susceptible  to  coupling 
distortion.  At  the  smallest  spacing,  s  =  0.6  mm,  the  signal  line 
amplitude  has  decreased  by  one-half  and  the  sense  line  has  risen  to  the 
same  magnitude  as  the  signal  line.  As  the  spacing  increases,  the 
amplitude  on  the  signal  line  increases  and  the  sense  line  decreases. 
When  s  =  1.8  mm,  the  signal  line  response  has  been  degraded  by  only  10 
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Fig.  4.22.  Pulse  distortion  on  coupled  lines  vs.  spacing,  sense  line  response 
1  =  150mm,  r  =  30ps,  beryllium  oxide  substrate. 
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percent,  but  the  sense  line  is  up  to  35  percent  of  the  initial 
amplitude  of  the  input  pulse.  Depending  on  the  system  specification, 
this  may  be  too  large  of  a  spurious  signal,  even  though  this  spacing 
satisfies  the  normal  design  criteria. 

Increasing  the  substrate  relative  dielectric  constant  to  6.8, 

Figs.  4.21  and  4.22,  and  9.7,  Figs.  4.23  and  4.24,  increases  dispersion 
and  coupling  distortion  even  as  the  distance  traveled  decreases.  At 
the  smaller  spacings,  s  =  0.6  mm  and  s  =  1.2  mm,  the  tighter  coupling 
on  this  structure  causes  both  the  signal  and  sense  line  responses  to 
separate  into  the  even  and  odd  mode  pulses,  i.e.  two  pulses  in  phase  on 
the  signal  line,  and  two  pulses  180  degrees  out  of  phase  on  the  sense 
line.  Due  to  the  distortion  of  the  pulse  from  dispersion,  the  maximum 
positive  sense  line  response,  which  is  greater  than  one-half  of  the 
amplitude  of  the  undistorted  signal,  occurs  when  the  spacing  is  1.8  mm 
and  is  slightly  less  for  the  smaller  spacings.  Normally,  it  is 
expected  that  the  sense  line  amplitude  would  not  exceed  one-half  of  the 
input  signal  because,  as  shown  in  section  3.2,  the  magnitude  of  both 
the  signal  and  sense  line  responses  tend  to  approach  one-half  the 
original  amplitude  as  the  even  and  odd  mode  pulses  separate  completely. 

The  larger  response  on  the  sense  line  is  a  result  of  the  increased 
dispersive  characteristics  of  the  line  due  to  the  higher  relative 
dielectric  constant.  This  can  be  visualized  by  approximating  the 
distortion  due  to  coupling,  as  in  Fig.  3.2,  as  being  a  result  only  of 
the  separation  of  the  even  and  odd  mode  pulse  pairs  over  distance  due 
to  differences  in  the  modal  phase  velocity.  To  include  the  effects  of 


dispersion,  the  undistorted  pulse  pairs  in  Fig.  3.2  are  replaced  with 
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one-half  of  the  isolated  line's  response.  This  is  only  approximate, 
because  the  difference  in  the  modal  phase  velocity  is  a  function  of 
frequency,  and  therefore  affects  each  frequency  component  of  the  pulse 
differently  and  also  because  the  dispersive  characteristics  of  the  even 
and  odd  modes,  while  being  very  similar,  are  not  quite  exactly  the 
same . 

In  Figs.  4.21  and  4.23,  the  isolated  line  response  is 
approximately  the  same  as  the  s  =  5.0mm  graphs.  In  both  graphs,  either 
=  6.8  or  =  9.7,  the  isolated  pulse  has  a  relatively  large 

negative  swing  on  the  trailing  edge  due  to  the  dispersion  distortion. 
Since  the  sense  line  response  is  equal  to  the  even  mode  response  minus 
the  odd  mode  response  (eq.  3.11a),  and  because  the  odd  mode  leads  the 
even  mode,  the  trailing  negative  swing  of  the  odd  mode  pulse  adds 
constructively  to  the  even  mode  pulse.  This  constructive  interference 
increases  the  sense  line  response  so  that  it  is  larger  than  one-half  of 
the  input  signal.  Therefore  the  maximum  possible  sense  line  response, 
for  a  given  spacing,  would  occur  when  the  maximum  of  the  negative  swing 
of  the  odd  mode  pulse  aligns  with  the  positive  maximum  of  the  even  mode 
pulse.  Since  the  odd  mode  leads  the  even  mode  in  this  structure,  there 
is  no  possibility  of  constructive  interference  in  the  sense  line' s 
negative  leading  response,  and  so  it  attains  a  maximum  when  the  line 
spacing  is  a  minimum.  Thus  in  order  to  obtain  accurate  predictions  of 
the  signal  and  sense  line  responses,  the  dispersive  characteristics  of 
the  line,  as  well  as  the  pulse  spectrum  must  be  accounted  for  in  a 
rigorous  manner,  such  as  the  Fourier  transform  approach  used  in  this 


report . 
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When  the  relative  dielectric  constant  of  the  substrate  is 
increased  to  12.2,  Figs.  4.25  and  4.26,  both  dispersion  and  coupling 
have  a  significant  effect  on  the  distortion  of  the  pulse.  The  signal 
and  sense  line  responses  on  the  gallium  arsenide  substrate  are  very 
similar  to  the  corresponding  pulse  shapes  for  the  beryllium  oxide  and 
alumina  substrates,  except  that  the  pulses  on  the  gallium  arsenide 
structure  are  affected  significantly  more  by  dispersion. 

The  s  =  5.0  mm  graph,  in  which  dispersion  distortion  is  the 
dominant  mechanism,  shows  a  12  percent  reduction  in  the  maximum 
amplitude  on  the  signal  line.  In  addition,  the  magnitude  of  the 
negative  trailing  swing  of  the  signal  line  has  risen  to  over  20  percent 
of  the  undistorted  amplitude.  Like  the  previous  two  substrate 
materials,  the  maximum  positive  sense  line  response  does  not  occur  when 
the  spacing  is  the  smallest  but  instead  happens  for  s  =  1.2  mm  or 
1 . 8  mm. 

When  s  =  1.8  mm,  the  signal  line  response  has  degraded  to  65 
percent  of  the  initial  amplitude  and  the  sense  line  response  has  risen 
to  a  maximum  magnitude  of  50  percent  of  the  undistorted  value.  For 
most  circuits,  these  would  be  unacceptable  values  for  coupling  loss  and 
the  amount  of  crosstalk,  even  though  this  spacing  is  within  the  three 
widths/heights  rule  of  thumb.  At  a  smaller  distance,  the  amount  of 
coupling  will  not  be  as  great,  but  clearly,  the  allowable  proximity  of 
the  transmission  lines  for  uncoupled  operation  depends  not  only  on  the 
spacing  of  the  center  conductors,  but  also  the  substrate  materials, 
structural  dimensions,  and  the  distance  over  which  the  lines  are 
coupled  to  each  other.  In  general,  for  single  layer  microstrips, 
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relatively  short  distances  and/or  small  dielectric  constants  allow 
smaller  spacings  while  longer  distances  and/or  larger  dielectric 
constants  require  larger  spacings. 

Using  the  Spectral  Domain  Approach  to  determine  the  line 
parameters  and  the  even/odd  mode  formulation  with  the  Fourier  transform 
approach  to  compute  the  time  domain  results,  the  transient  analysis  of 
coupled  lines  has  been  analyzed,  showing  the  effects  of  both  dispersion 
and  even/odd  mode  coupling  distortion.  First,  the  frequency  dependent 
Ereff  °f  couPlec*  lines  were  presented  for  both  the  even  and  the  odd 

modes  with  different  line  spacings  and  substrate  materials.  These 
results  were  then  used  to  examine  pulse  distortion  on  tightly  coupled 
lines  at  different  distances,  showing  how  both  dispersion  and  even/odd 
mode  coupling  distort  the  input  pulse  and  produce  a  spurious  response 
on  the  adjacent  line.  Next,  pulse  distortion  was  studied  as  a  function 
of  the  spacing  between  the  center  conductors.  This  analysis  showed  how 
the  dispersive  characteristics  of  the  lines  can  cause  constructive 
interferences  on  the  sense  line,  giving  a  larger  response  than 
anticipated.  In  addition,  these  results  showed  that  not  only  the 
separation  of  the  center  conductors  needs  to  be  considered  in  designing 
decoupled  lines,  but  equally  important  factors  are;  the  substrate 
material,  structural  dimensions,  and  the  distance  over  which  the  lines 
are  close  to  each  other. 

4.2  Multilayer  Structures 

The  desire  to  improve  the  performance  of  MMIC' s  and  the  increasing 


interest  in  them  has  created  a  demand  for  the  accurate  characterization 
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of  multi-layer  structures,  and  therefore  requires  the  use  of  rigorous 
analytical  techniques.  The  desire  to  equalize  the  even  and  odd  mode 
phase  velocities  in  coupled  lines,  as  well  to  limit  the  space 
requirements,  have  made  multilayer  structures  attractive  for  many 
applications  as  well.  A  full  wave  analysis  is  necessary  because  the 
addition  of  either  substrates  or  superstrates  produces  significant 
changes  in  the  frequency  dependent  parameters  that  cannot  be  predicted 
accurately  by  quasi-static  or  other  approximate  formulations. 

This  section  uses  the  generalized  Green's  function  in  conjunction 

with  the  SDA  to  compute  the  frequency  dependent  parameters  of  both 

single  and  coupled  lines  for  multilayer  structures.  First  the 

effective  dielectric  constant  is  computed  as  a  function  of  frequency 

for  different  combinations  of  substrate  materials  and  heights  to 

illustrate  the  effect  of  the  addition  of  substrates  on  the  line's 

performance.  Next  a  superstrate  of  varying  height  with  the  same 

relative  dielectric  constant  as  the  substrate  is  added  to  an  open 

coupled  line  structure  to  consider  the  effects  of  superstrates.  Then 

as  the  height  of  the  superstrate  layer  is  fixed,  an  upper  ground  plane 

is  introduced,  and  the  c  __  for  the  isolated,  even,  and  odd  mode  cases 

ref  f 

is  plotted  as  the  height  of  the  ground  plane  is  varied.  Finally,  an 

open  symmetric  coupled  line  structure  with  two  substrates  is  used  to 

show  how  exchanging  the  two  relative  dielectric  constants  of  the 

substrates  changes  both  c  and  pulse  propagation  when  all  other 

reff 

dimensions  are  held  constant. 

Increasing  the  number  of  layers,  either  substrates  or 
superstrates,  can  create  a  significant  difference  in  the  effective 
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dielectric  constant.  An  open  single  microstrip  with  two  substrates  is 

used  in  Figs.  4.27  through  4.30  to  illustrate  the  effects  of  these 

additional  substrates.  The  total  substrate  height  is  held  constant  and 

six  different  combinations  of  h  and  h  are  used.  The  h  =0.6 

L2  LI  L2 

mm  graph  represents  a  single  substrate  structure  with  a  relative 

dielectric  constant  of  e  .  Likewise  the  h  =0.0  mm  graph  also 

rL2  L2 

represents  a  single  substrate  line,  but  with  a  relative  dielectric 

constant  of  crL1  •  The  other  graphs  represent  two  substrate  structures 

where  the  heights  of  the  substrates  are  h  and  h  ,  with  relative 
y  L2  Ll 

dielectric  constants  of  e  and  e  ,  respectively. 

rL2  rLl  J 

The  effective  dielectric  constant  begins  at  the  single  layer  value 

for  e  when  h  =0.0  mm,  or  zero  percent  of  the  total  height.  Since 
rLl  L2 

the  dielectric  substrate  being  added  has  a  lower  relative  dielectric 

constant  than  the  substrate  already  present,  as  h  increases,  e 

L2  reff 

begins  to  decrease.  The  h  values  of  0.1,  0.25,  0.35,  and  0.5  mm 

L2 

represent,  respectively,  changes  of  16.6,  41.6,  58.3  and  83.3  percent. 


When  h  =  0.6  mm,  or  100  percent  of  the  total  height,  Ereff  is 

identical  to  the  single  layer  value  for  e  .  Thus  the  expectation  is 

that  e  __  will  decrease  monotonically  as  h  increases  and  the  total 
reff  L2 


change  in  Ereff  at  a  given  frequency  will  be  equal  to 

Ac  „(f)  =  eL*(f)  -  eL2„(f) 
reff  reff  reff 


(4.4) 


where  the  Ll  and  L2  superscripts  indicate  that  e  „„  is  evaluated  for  a 

reff 

single  substrate  microstrip  with  a  relative  dielectric  constant  of  e  ^ 
and  e  ,  respectively. 


The  first  structure,  with  its  results  shown  in  Fig.  4.27,  consists 


of  a  layer  of  RT/duroid  on  top  of  a  layer  of  beryllium  oxide.  At  low 
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Fig.  4.27.  e 


Frequency  m  GHz 

vs.  frequency  for  open  multilayer  single  microstrip  on 


beryllium  oxide  and  RT/duroid  5880  substrate. 


Frequency  in  GHz 

4.29.  e  vs.  frequency  for  open  multilayer  single  microstrip  on 
l*  cff 

gallium  arsenide  and  beryllium  oxide  substrate. 
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Frequency  in  GHz 

Fig.  4.30.  £reff  vs-  frequency  for  open  multilayer  single  microstrip  on 
gallium  arsenide  and  alumina  substrate. 
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frequencies,  the  smallest  duroid  layer,  h  =0.1  mm,  decreases  e 

L2  reff 

by  55  percent  of  Ae  ^  even  though  the  height  added  is  only  16.6 

percent  of  the  total  height.  When  h^  =  0.25  mm,  or  41.6  percent  of 

the  total  height,  the  decrease  in  e  is  almost  80  percent.  In  Fig. 

4.28,  which  represents  data  for  a  structure  with  duroid  on  alumina,  the 

difference  at  low  frequencies  is  even  larger;  that  is  the  16.6  percent 

height  change  gives  a  62  percent  change  in  e  and  the  41.6  percent 

height  change  results  in  a  difference  in  e  __  of  85  percent  of  Ac  „. 

reff  reff 

Thus  the  addition  of  a  relatively  thin  substrate  can  have  a  profound 

impact  on  the  effective  dielectric  constant  of  a  multilayer  structure, 

especially  for  the  lower  frequencies. 

The  third  structure,  gallium  arsenide  on  beryllium  oxide,  shown  in 

Fig. 4. 29,  gives  results  that  are  similar  to  the  previous  ones.  The 

addition  of  the  0.1  mm  layer  gives  a  44.7  percent  decrease  while  the 

change  of  h  from  0.5  mm  to  0.6  mm  results  in  only  a  change  of  only  6 

percent  of  Ac  „„  at  low  frequencies.  Therefore,  while  the  addition  of 
reff 

a  thin  upper  substrate  layer  with  a  lower  c^_  gives  large  changes  in 

e  ,  .  the  addition  of  an  equally  thin  lower  substrate  layer  with 
reff 

higher  c  results  in  relatively  small  changes  in  c  Also,  as  the 

r  reff 

beryllium  oxide  layer  is  increased  in  height,  the  dispersive 

characteristics  of  the  structure  decrease.  This  is  expected  because 

single  layer  structures  are  less  dispersive  when  the  relative 

dielectric  constant  of  the  substrate  is  less. 

The  final  two  substrate  configuration  considered  has  a  layer  of 

alumina  on  a  layer  of  gallium  arsenide.  The  results,  plotted  in  Fig. 

4.30,  do  not  have  as  drastic  of  changes  in  c  as  those  of  the 

reff 
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previous  graphs,  because  the  relative  dielectric  constants  of  the 

substrates  are  closer  to  each  other  than  those  of  the  previous 

structures.  At  low  frequencies,  the  0.1  mm  and  0.25  mm  graphs  have 

changes  of  only  38  and  64  percent  of  ^creff  respectively.  However, 

like  the  other  structures,  as  frequency  increases,  the  percent  changes 

in  eref£  become  closer  to  the  percent  changes  in  the  layer  height.  At 

100  GHz,  the  0.1  and  0.25  mm  structures  have  changes  in  c  ,,  of  16.3 

ref  f 

and  36.2  percent  of  Ae  which  are  much  closer  to  the  height  changes 

reff 

of  16.7  and  41.6  percent  of  the  total  height.  Additionally,  the 

increase  in  h  to  0.1  mm  (h  =  0.5  mm)  and  0.25  mm  (h  =  0.35  mm), 

gives  changes  of  16.8  and  46.7  percent  respectively,  which  are  much 

closer  to  the  corresponding  changes  for  h  .  Note  that  while  at  low 

frequencies  the  addition  of  a  thin  upper  substrate  gives  a  much  larger 

change  in  e  „„  than  the  addition  of  a  thin  lower  substrate,  as  the 
reff 

frequency  is  increased,  the  addition  of  a  lower  substrate  gives  about  the 
same  change  as  the  addition  of  an  upper  substrate. 

This  change  in  the  effect  of  upper  and  lower  substrates  as  the 
frequency  is  increased  can  be  attributed  to  changes  in  the  electric 
field  structure  as  the  frequency  increases  beyond  the  quasi-static 
region.  As  the  frequency  is  increased,  the  electric  field  lines  tend 
to  bunch  underneath  the  center  conductor,  connecting  directly  with  the 
ground  plane,  rather  than  fringing  out  into  the  upper  layers.  This 
results  in  a  field  configuration  that  is  similar  to  a  parallel  plate 
capacitor  filled  with  two  different  dielectric  slabs.  Since  the 
parameters  of  the  parallel  plate  capacitor  are  independent  of  the 


relative  placement  of  the  slabs,  then  the  microstrip  parameters,  in  the 
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high  frequency  region,  will  also  be  independent  of  the  substrate 
placement . 

Since  at  very  high  frequencies  the  microstrip  may  be  approximated 

by  a  parallel  plate  capacitor,  it  is  instructive  to  consider  the 

equivalent  dielectric  constant  of  a  two-layer  parallel  plate  capacitor. 

The  equivalent  dielectric  constant,  c  ,  of  the  dual  dielectric 

eqv 

parallel  plate  capacitor  is  the  dielectric  constant  of  a  single  layer 

of  height  h  +h  that  when  used  in  the  parallel  plate  capacitor  gives 

the  same  capacitance  as  one  filled  with  two  layers  that  have  dielectric 

constants  c  and  e  and  heights  h  and  h  .  Using  elementary 

circuit  theory  and  the  capacitance  formula  for  a  parallel  plate 

capacitor,  the  equivalent  dielectric  constant  is  found  to  be 
h  +h 


L2  Ll 


'eqv  h  /e  +h  /e 

L2  L2  Ll  Ll 


e  is  also  the  high  frequency  limit  of  e  i.e. 

eqv  reff 


lim  c  _,(f)  =  c 
,  reff  eqv 

f-*D 


(4.5) 


(4.6) 


Note  that  this  equivalence  applies  only  at  frequencies  that  are 
sufficiently  high  so  that  almost  all  the  electric  field  lines  lie 
directly  between  the  center  conductor  and  the  ground  plane  with  almost 
no  fringing.  This  concept  of  the  equivalent  dielectric  constant  can  be 
extended  to  include  more  dielectric  layers  by  characterizing  each  layer 
as  another  parallel  plate  capacitor  attached  in  series. 

Since  the  difference  in  the  even  and  odd  mode  phase  velocity 
degrades  the  performance  of  a  circuit,  it  is  advantageous  to  be  able  to 
equalize  the  modal  phase  velocities.  One  method  that  has  been  used  in 


the  past  to  accomplish  this  is  to  place  a  superstrate  layer  with  the 
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same  dielectric  constant  as  the  substrate  on  top  of  the  center 
conductors.  If  the  height  of  this  superstrate  is  adjusted  properly, 
then  the  even  and  odd  mode  phase  velocities  will  be  equalized.  The 
concept  behind  this  approach  is  that  striplines,  which  have  a 
homogeneous  dielectric  medium,  are  TEM  structures  where  the  even  and 
odd  mode  phase  velocities  are  equal.  Thus  the  placement  of  such  a 
layer  of  the  correct  thickness  would  make  the  microstrip  more  like 
the  stripline,  i.e.  a  more  homogeneous  medium,  and  hence  equalize  the 
even  and  odd  mode  phase  velocities . 

This  method  of  equalizing  the  modal  velocities  is  now  applied  to 
structures  with  four  different  substrates  and  the  results  are  shown  in 
Figs.  4.31  through  4.34.  The  configuration,  which  is  shown  in  each  of 
the  figures,  is  an  open  symmetric  coupled  microstrip  with  a  single 
substrate  and  two  superstrates  where  h  — >  oo.  All  four  structures  are 
considered  at  f  =  100  MHz,  which  is  in  the  quasi-static  region  for 
these  microstrips.  This  frequency  was  chosen  because  it  gives  very 
broadband  characteristics,  from  0  to  about  1  GHz,  and  because  many 
microstrip  designs  operate  in  or  near  the  quasi-static  region. 

In  each  of  the  four  graphs,  when  the  superstrate  layer  is  very 
thin,  there  is  not  much  effect  on  the  effective  dielectric  constant. 
However,  when  the  thickness  of  the  superstrate  is  increased  to  10 
percent  of  the  substrate  height,  then  e  ^  for  ail  three  cases, 
isolated,  even  mode,  and  odd  mode,  begins  to  increase  as  well.  The 
effect  of  the  superstrate  layer  is  more  pronounced,  however,  if  its 
dielectric  constant  is  much  larger  than  that  of  the  substrate.  The  odd 


mode  effective  dielectric  constant,  however,  increases  much  faster  in 
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Fig.  4.31.  c  vs.  cover  layer  height  for  multilayer  symmetric  coupled 
microstrips  with  RT/duroid  5880  superstrate  and  substrate. 
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Fig.  4.32.  c  vs.  cover  layer  height  for  multilayer  symmetric  coupled 
microstrips  with  beryllium  oxide  superstrate  and  substrate. 


reff 


Fig.  4.34.  £  vs.  cover  layer  height  for  multilayer  symmetric  coupled 

microstrips  with  gallium  arsinide  superstrate  and  substrate. 
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this  region  than  the  other  two  e  ,  ,'s.  When  h  equals  infinity,  then 

the  microstrips  are  imbedded  in  a  completely  homogeneous  medium  with  a 

relative  dielectric  constant  of  c  =  c  .  Therefore,  in  the  limit 

rLl  rU2 

as  h  goes  to  infinity,  the  effective  dielectric  constants  for  all 

three  cases  approach  the  relative  dielectric  constant  of  the  substrate. 

For  the  four  structures  used,  the  infinite  height  condition  is 

approximately  satisfied  when  h  is  about  10  times  h^,  i.e.  when 

h  =10  mm. 

U2 

Although  the  modal  phase  velocities  are  approximately  equal  when 

h  is  about  10  times  h  ,  this  is  not  a  very  efficient  design  due  to 
U2  LI 

the  extraordinarily  large  thickness  of  the  superstrate.  Aside  from  the 
expense  of  such  a  circuit,  the  large  dielectric  slab  could  support 
spurious  guided  wave  modes  at  relatively  low  frequencies  that  would  be 
excited  by  discontinuities  in  the  lines.  However,  because  the 
effective  dielectric  constant  of  the  odd  mode  increases  much  faster 
than  the  even  mode  for  each  of  the  four  structures,  at  one  other  h 

U2 

value  besides  infinity,  the  even  and  odd  mode  phase  velocities  are 

equal.  For  all  four  materials  being  considered,  this  occurs  when  h 

is  approximately  the  same  height  as  the  substrate.  The  even/odd  mode 

equalization  occurs  when  h  =  0.97  mm  for  the  duroid  and  moves 

U2 

slightly  higher  as  the  dielectric  constant  increases.  For  the 

beryllium  oxide  substrate  and  superstrate  combination,  the  equalization 

occurs  when  h  is  equal  to  1.06  mm  and  at  around  1.07  mm  for  both  the 
U2 

alumina  and  gallium  arsenide  structures. 

The  motivation  behind  adding  a  superstrate  layer  to  the  microstrip 
was  to  make  the  structure  more  TEM,  like  the  stripline,  by  making  the 
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dielectric  surrounding  the  center  conductors  seem  approximately 

homogeneous.  For  the  structures  considered,  this  approximation  to  the 

stripline  occurred  when  the  superstrate  layer  was  about  ten  times  the 

substrate  layer.  This  might  have  been  anticipated  because  a  shielded 

microstrip  has  approximately  the  same  characteristics  as  an  open 

microstrip.  This  is  more  evident  when  the  distance  to  the  upper  ground 

plane  and  side  walls  are  at  least  five  times  the  substrate  height  or 

center  conductor  width  [46].  However,  the  even/odd  mode  equalization 

also  occurred  for  a  much  smaller  h  value,  a  value  at  which  the 

U2 

structure  does  not  approximate  the  stripline  and  is  not  TEM.  Also  note 

that  when  h  is  greater  than  this  first  equalization  point,  the  odd 

mode's  e  __  is  greater  than  the  even  mode's  e  __ ,  something  that  did 
reff  reff 

not  occur  in  the  single  layer  structure.  This  indicates  that  a 

completely  different  mechanism  is  responsible  for  the  equalization  of 

the  even  and  odd  mode  phase  velocities. 

As  long  as  the  microstrip  is  operated  in  the  quasi-static  region, 

then  the  location  of  the  equalization  point  stays  constant,  because 

c  , ,  is  constant  for  the  frequencies  in  this  region.  However,  as  the 
reff 

frequency  increases  past  the  quasi-static  region  of  the  structure,  the 
location  of  this  equalization  point  moves  to  larger  values  of  the 
superstrate  height.  At  some  frequency  outside  the  quasi-static  region, 
the  equalization  point  changes  to  a  value  so  large  that  it  gives  no 
advantage  in  design,  essentially  eliminating  it  as  a  possibility.  This 
occurs  because  as  frequency  increases,  the  electric  field  lines  tend  to 
go  in  a  straight  line  from  the  center  conductor  to  the  ground  plane 
with  less  and  less  lines  in  the  superstrate  regions.  Thus  the 
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superstrate  has  a  decreasing  effect  on  cre^^  as  the  frequency 
increases . 

While  the  inability  to  find  a  realizable  superstrate  height  that 
matches  the  even  and  odd  mode  phase  velocities  at  high  frequencies  is  a 
significant  drawback  to  this  method,  many  microwave  circuits  operate  in 
the  quasi-static  region,  and  so  there  are  applications  that  can  use 
this  concept.  In  addition,  it  will  be  shown  in  section  4.3  that  this 
even/odd  mode  equalization  phenomenon  can  also  occur  when  multiple 
substrates  are  used  in  the  structure  instead  of  an  additional 
superstrate.  Since  the  substrates  lie  between  the  center  conductor  and 
the  ground  plane,  as  the  frequency  increases,  they  still  have  a 
significant  effect  on  the  electric  field  lines  of  the  structure, 
increasing  the  potential  bandwidth  of  the  design.  Also,  since  the  even 
and  odd  mode  phase  velocities  tend  to  approach  each  other  at  very  high 
frequencies,  it  will  be  possible  to  achieve  a  very  good  match  of  the 
even  and  odd  mode  phase  velocities  over  an  almost  infinite  bandwidth. 

The  necessity  to  control  electromagnetic  emissions,  the  desire  to 
cutoff  higher  order  modes,  as  well  as  the  need  to  reduce  radiation 
losses,  often  necessitate  shielding  the  microstrip  circuit.  To  be  able 
to  accurately  account  for  the  effects  of  shielding  on  the  transmission 
lines  requires  the  use  of  a  rigorous  analytical  approach.  Fortunately, 
due  to  the  general  manner  in  which  the  generalized  Green's  function  was 
derived,  as  well  as  the  flexibility  of  the  spectral  domain  approach,  it 
is  as  easy,  if  not  easier,  to  compute  the  frequency  dependent 
parameters  for  shielded  structures  as  for  open  ones.  Figures  4.35 


through  4.38  deal  with  the  presence  of  an  upper  ground  plane,  or  cover 
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sheet,  in  a  symmetric  coupled  microstrip  structure  with  one  substrate 

and  two  superstrates .  The  same  four  dielectric  materials  are  used  and 

the  effective  dielectric  constant  for  the  isolated  line  and  the  even 

and  odd  modes  of  the  coupled  lines  are  computed  as  the  height  of  the 

upper  ground  plane,  h^,  is  varied. 

As  h  gets  smaller,  then  the  structure  becomes  more  like  a 

stripline,  and  so  the  effective  dielectric  constant  for  all  three  cases 

approaches  the  relative  dielectric  constant  of  the  substrate.  On  the 

other  hand,  as  h  becomes  large,  then  the  structure  becomes  more  like 
ui 

an  open  microstrip,  and  the  effective  dielectric  constants  for  the 

isolated,  even  and  odd  modes  approach  approach  the  corresponding  values 

shown  in  Figs.  4.31  through  4.34  when  h  =  0.3  mm.  As  with  the 

previous  four  structures,  for  generality,  the  frequency  at  which  the 

parameters  are  computed  is  100  MHz. 

The  most  noticeable  characteristic  in  all  four  graphs  is  that  the 

effective  dielectric  constants  vary  over  a  large  range  of  values  as  the 

height  is  changed  instead  of  monotonically  decreasing  from  the 

stripline  value  to  the  open  microstrip  value.  ereff  begins  at  the 

maximum  value  when  h  is  the  smallest  and  begins  to  decrease  in 

ui  y 

magnitude  until  a  minimum  is  reached  when  h  is  somewhere  between  0.1 

and  0.4  mm,  depending  on  the  dielectric  material  used.  As  h 

continues  increasing  past  this  value,  e  __  begins  to  increase  as  well, 

ref  f 

approaching  the  open  microstrip  value  when  h  is  about  10  times  the 
substrate  height,  i.e.  when  h  “  10  mm. 
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Fig.  4.35.  t  vs.  cover  sheet  height  for  multilayer  symmetric  coupled 
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microstrips  with  RT/duroid  5880  superstrate  and  substrate. 
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Fig.  4.36.  £  ff  vs.  cover  sheet  height  for  multilayer  symmetric  coupled 
microstrips  with  beryllium  oxide  superstrate  and  substrate. 
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Fig.  4.37.  t  ft  vs.  cover  sheet  height  for  multilayer  symmetric  coupled 
microstrips  with  alumina  superstrate  and  substrate. 
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The  manner  in  which  e  changes  as  h  varies  can  be  understood 

reff  ui 

by  considering  the  field  structure  in  shielded  configuration.  When  h 
is  very  small,  then  the  paths  of  the  field  lines  exist  almost  entirely 
in  the  dielectric  material  with  only  very  short  distances  in  the  air, 
and  so  creff  is  approximately  equal  to  the  e^_  of  the  substrate.  As  the 
air  gap  increases  in  height,  then  the  field  lines  that  connect  with  the 
cover  layer  pass  through  an  air  layer  whose  thickness  is  no  longer 
negligible  and  therefore  has  greater  effect  on  the  fields.  Because  the 
dielectric  constant  of  the  air  is  lower  than  that  of  the  substrate  or 
superstrate,  cref£  begins  to  decrease  in  value.  At  some  point, 
however,  the  upper  ground  plane  is  sufficiently  distant  so  that  more  of 
the  electric  field  lines  go  through  the  substrate  to  ground  rather  than 
through  the  superstrates,  bypassing  the  air  gap  and  negating  its 
effects.  This  occurs  even  though  the  distance  from  the  center 
conductors  to  the  upper  ground  plane  is  shorter  than  the  distance  to 
the  lower  ground  plane.  But,  because  the  substrate  has  a  higher 
dielectric  constant  than  the  equivalent  combination  of  the  upper  two 
layers,  it  attracts  more  of  the  field  lines. 

As  the  cover  sheet  moves  farther  away  from  the  center  conductors, 
less  field  lines  go  from  the  center  conductors  through  the  air  layer  to 
the  cover  sheet,  thereby  increasing  e  .  Finally,  when  the  cover 
sheet  is  far  enough  away,  it  attracts  virtually  no  field  lines  and  no 
longer  has  any  noticeable  effect  on  the  structure.  When  this  happens, 
for  all  three  cases  approach  the  corresponding  values  for  the 


open  structure. 
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The  h^  value  corresponding  to  the  minimum  c  value  varies 

slightly  as  the  relative  dielectric  constant  of  the  substrate  and 
superstrate  changes.  Additionally,  while  the  minimum  point  is 
approximately  the  same  for  each  of  the  three  cases,  it  varies  slightly, 
with  the  odd  mode  having  the  lowest  h  value  and  the  even  mode  the 
highest.  For  the  RT/duroid  substrate  and  superstrate  combination,  Fig. 

4.35,  the  minima  occur  when  h  is  about  0.43  mm.  If  the  relative 

ui 

dielectric  constant  of  the  material  was  decreased,  it  would  be  expected 
that  the  minimum  point  would  approach  0.7  mm,  so  that  the  total 
superstrate  height  would  be  equal  to  the  height  of  the  substrate. 

Also,  the  minimum  points  of  each  of  the  three  cases  would  tend  to 
approach  each  other  more  closely  with  the  even  and  odd  mode  e 
approaching  the  isolated  value,  which  in  turn  would  be  approaching 
unity.  In  the  limit  as  the  relative  dielectric  constant  approached 
unity,  the  graphs  would  flatten  completely  and  no  distinct  minimum 
would  exists. 

As  the  dielectric  constant  increases,  however,  the  minima  points 

occur  at  lower  values  of  h  ,  and  the  locations  for  each  of  the  three 

ui 

cases  increase  in  separation.  For  the  beryllium  oxide  structure.  Fig. 

4.36,  the  minima  occur  for  h  near  0.21  mm,  while  for  the  alumina 

ui 

substrate  and  superstrate,  Fig.  4.37,  they  occur  near  0.19  mm.  When 
gallium  arsenide  is  used,  Fig.  4.38,  the  minima  occur  near 
h  =  0.16  mm,  but  differ  noticeably  for  each  case.  The  odd  mode's 
minimum  occurs  at  h  “  0.155  mm  while  the  even  mode's  happens  around 
0.175  mm.  As  would  be  expected,  the  minimum  for  the  isolated  case  ends 

up  between  that  of  the  even  and  the  odd  mode  at  h  =*  0.165  mm. 

ui 


125 


In  addition  to  the  locations  of  the  minima  moving,  increasing  the 

relative  dielectric  constant  of  the  material  also  makes  the  change  in 

e  __  much  greater  as  well.  For  the  duroid  structure,  the  introduction 
ref  f 

of  a  0.01  mm  air  gap  decreases  e  by  about  2  percent  and  the  0.1  mm 

reff 

gap  decreases  it  by  12  percent.  The  decrease  at  the  minimum  point  is 

17  percent  lower  than  the  stripline  value  and  7  percent  lower  than  the 

open  microstrip  value.  When  the  dielectric  is  changed  to  gallium 

arsenide,  the  changes  in  c  __  become  much  more  drastic.  The  0.01  mm 

reff 

and  0.1  mm  air  gaps  decrease  e  by  14  and  40  percent,  respectively. 

The  minimum  value  of  c  is  now  40  percent  lower  than  the  stripline 

reff 

value  and  23  percent  lower  than  the  open  microstrip  case. 

Thus  if  a  stripline  circuit  were  designed  with  similar  dimensions 
on  gallium  arsenide,  but  in  production  there  were  a  10  micron  gap 
between  the  upper  ground  plane  and  the  dielectric  material,  there  would 
be  and  8  percent  increase  in  the  wavelength  and  hence  an  equivalent 
decrease  in  the  electrical  dimensions  of  the  circuit,  since 
A  oc  1  /^e  If  the  air  gap  were  0.1  mm,  however,  the  wavelength 

would  be  increased  by  23  percent.  Changes  in  the  electrical  dimensions 
of  the  circuit  of  this  magnitude  would  certainly  have  an  adverse  effect 
on  the  circuit  performance.  If,  instead,  the  circuit  had  been  designed 
as  an  open  structure  without  accounting  for  the  shielding,  then, 
depending  on  the  placement  of  the  cover  sheet,  there  could  also  be 
significant  changes  in  the  parameters  of  the  lines.  Clearly  the 
presence  of  a  cover  sheet,  which  is  very  common  in  production  designs, 
must  be  taken  into  account  using  a  rigorous,  full  wave  technique  like 


the  SDA. 
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As  with  the  previous  dual  superstrate  structures,  it  is  possible 

to  adjust  the  dimensions  to  equalize  the  even  and  odd  mode  phase 

velocities.  For  the  four  structures  under  consideration,  a  cover  layer 

height  of  slightly  less  than  2  mm  makes  the  even  and  odd  mode  e  „„  the 

reff 

same.  However,  also  like  the  two  superstrate  open  structure,  as  the 

frequency  increases  above  the  quasi-static  region,  the  h  value 

required  to  maintain  the  even/odd  mode  equalization  decreases.  As  the 

frequency  increases,  the  electric  field  lines  tend  to  gather  along  the 

path  with  the  highest  capacitance.  In  order  to  maintain  enough  of  the 

field  lines  in  the  superstrates,  the  capacitance  of  the  upper  layers 

must  be  increased,  and  so  the  cover  sheet  height  must  be  reduced  to 

maintain  the  even/odd  mode  equalization.  Due  to  the  large  change  in 

the  required  cover  sheet  height  for  frequencies  above  the  quasi-static 

region,  this  method  of  equalizing  the  even/odd  mode  phase  velocities 

can  only  be  used  either  in  the  quasi-static  region  or  for  very 

narrowband  applications  at  frequencies  beyond  the  quasi-static  region. 

It  was  shown  previously  that  the  presence  of  multiple  substrates 

in  the  microstrip  structure  causes  significant  deviations  in  creff  from 

the  single  layer  value.  However,  the  effect  of  the  relative  placement 

of  the  substrates  has  not  been  dealt  with  yet.  At  very  high 

frequencies,  the  microstrip  can  be  approximated  by  a  parallel  plate 

capacitor,  and  so  c  „  (oo)  is  determined  by  the  heights  and  relative 

reff 

dielectric  constants  of  the  substrates,  and  not  by  the  relative 
locations  of  the  substrates.  At  lower  frequencies,  however,  the 
position  of  a  substrate  layer  relative  to  the  other  substrate  layers 


can  have  an  important  effect  on  the  line  parameters. 
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To  illustrate  the  effect  of  substrate  placement  on  creff/  an  open 

symmetric  coupled  microstrip  structure  is  used  with  two  substrates, 

alumina  and  RT/duroid  5880.  The  first  configuration,  designated 

structure  #1,  places  a  0.3  mm  duroid  layer  on  top  of  an  alumina  layer 

of  the  same  height.  The  next  configuration,  called  structure  #2, 

switches  the  substrates,  placing  the  alumina  on  top  of  the  duroid, 

again  with  both  layers  having  the  same  height.  The  other  dimensions 

are  the  same  for  both  structures  and  are  listed  in  Fig.  4.39.  The 

effective  dielectric  constants  for  the  two  configurations  are  plotted 

in  Fig.  4.39  as  a  function  of  frequency  for  the  isolated  case  and  for 

the  even  and  odd  modes  of  the  coupled  lines. 

Structure  #1  has  a  much  lower  e  than  structure  #2,  about  36 

reff 

percent  lower,  indicating  that  structure  #1  should  be  less  dispersive 

than  #2.  On  the  other  hand,  the  separation  in  the  even  and  odd  mode 

phase  velocities  is  much  greater  in  #1,  16  percent,  compared  with  the 

5.6  percent  separation  for  structure  #2.  Thus  structure  #2  will  have 

much  less  distortion  due  to  coupling  than  #1.  Also  while  c  __  for  the 

reff 

odd  mode  is  lower  than  the  even  mode  in  #1,  switching  the  substrate 

layers  as  in  #2,  results  in  the  even  mode  e  „„  being  much  lower  than 

reff 

the  odd  mode . 

Pulse  distortion  on  these  structures  is  shown  in  Figs.  4.40  and 
4.41  for  a  Gaussian  pulse  with  t  =  30  ps  at  a  distance  of  250  mm.  The 
signal  line  response  on  structure  #1  has  been  critically  distorted  and 
has  separated  into  the  even  and  odd  mode  pulses.  On  the  other  hand, 
the  signal  line  response  on  structure  #2,  while  greatly  reduced  in 


magnitude,  is  not  nearly  as  distorted  as  the  response  on  #1  because  of 


reff 


Magnitude  of  pulse 


Fig.  4.40.  Pulse  distortion  on  multilayer  coupled  lines,  1  =  250mm, 
r  =  30ps,  structure  #1. 
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the  smaller  separation  in  the  even  and  odd  mode  e  __  for  this 

ref  f 

configuration.  The  sense  line  responses  for  the  two  structures  are 
also  markedly  different.  The  most  important  difference  between  the  two 
responses  is  that  structure  #1  has  a  negative  leading  sense  line 
response,  whereas  structure  #2  has  a  positive  leading  response.  This 
difference  arises  because  the  odd  mode  is  faster  than  the  even  mode  on 
structure  #1,  while  the  even  mode  is  the  faster  of  the  two  on  structure 
#2.  Therefore,  changing  the  relative  placement  of  the  the  substrate 
layers  has  a  significant  effect  on  the  characteristics  of  the 
structure . 

Using  the  generalized  Green's  function  in  conjunction  with  the 

spectral  domain  approach  allowed  the  accurate  analysis  of  complex 

multilayer  microstrip  structures.  Initially,  it  was  shown  that  the 

addition  of  even  relatively  thin  substrate  layers  with  different 

dielectric  constants  resulted  in  very  large  changes  in  e  f  for  the 

isolated  line,  even  mode,  and  odd  mode.  Then  attention  was  turned  to 

the  addition  of  superstrate  layers,  particularly  with  respect  to  the 

use  of  additional  superstrates  to  equalize  the  even  and  odd  mode  phase 

velocities.  While  it  was  possible  to  match  the  even  and  odd  modes  in 

the  quasi-static  region,  outside  of  this  region  the  match  was  either 

very  narrowband  or  required  impracticably  large  superstrate  heights. 

Next  a  cover  layer  was  introduced  in  the  two  superstrate  structure  and 

it  was  shown  that  the  addition  of  this  shielding  layer  created  a  wide 

variation  in  c  __  for  each  of  the  three  cases.  Finally,  the  effect  of 
ref  f 

changing  the  positions  of  the  substrates  was  examined,  showing  that  the 
relative  position  of  the  substrates  affects  not  only  the  magnitude  of 
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the  isolated  case,  but  also  can  cause  the  even  mode  c  to  be  lower 

ref  f 

than  the  odd  mode,  something  not  possible  with  single  substrate 
structures . 

4.3  Zero  Coupling  Structures 

It  was  shown  in  the  previous  section  that  the  dimensions  of  a 
multiple  superstrate  structure  could  be  adjusted  to  equalize  the  even 
and  odd  mode  phase  velocities.  Unfortunately  this  type  of  design  has  a 
very  narrow  bandwidth  and  matches  may  not  be  practical  at  certain 
frequencies.  The  bandwidth  limitation  arises  because  the  field 
structure  that  is  responsible  for  the  match  lies  in  the  superstrates 
and,  as  frequency  increases,  most  of  the  fields  move  into  the  substrate 
layers.  Thus  if  the  even/odd  mode  matching  conditions  could  be 
obtained  using  the  substrate  layers,  then  it  should  be  possible  to 
obtain  a  more  broadband  match. 

In  Fig.  4.39,  the  odd  mode  e  __  was  higher  than  the  even  mode 

reff 

c  ..  for  structure  #2.  However,  if  the  same  structure  were  completely 
reff 

filled  with  either  of  the  two  dielectric  materials  instead  of  a 

combination  of  the  two,  then  the  odd  mode  c  „„  would  be  less  than  the 

reff 

even  mode,  as  is  the  case  with  all  single  layer  structures.  Therefore 

there  must  exist  a  combination  of  h  and  h  which  would  equalize  the 

Ll  L2 

even  and  odd  mode  phase  velocities.  The  existence  of  these  even/odd 
mode  equalization  structures,  or  zero  coupling  structures,  is  due  to  a 
mechanism  that  is  similar  to  the  one  that  was  found  in  the  multiple 
superstrate  structures.  In  order  for  this  match  to  exist,  the 
dielectric  constant  of  the  lower  substrate  must  be  less  than  that  of 
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the  upper  substrate.  In  addition,  it  will  be  shown  that  the  difference 

in  the  dielectric  constants  must  surpass  some  value  that  depends  in 

some  way  on  all  of  the  dimensional  parameters  of  the  microstrip 

structure.  Thus  even  if  e  is  less  than  c  it  may  not  be  possible 

to  equalize  the  even  and  odd  mode  phase  velocities. 

To  study  the  characteristics  of  the  zero  coupling  structures,  an 

open  symmetric  coupled  microstrip  is  used  with  two  substrates  below  and 

air  above.  The  total  substrate  height  of  the  structure,  designated 

h  ,  ,  is  held  constant  while  the  heights  of  the  two  layers  are 
total 

varied.  When  h  is  equal  to  zero,  then  h  is  equal  to  h  ,  and  the 
li  M  L2  M  total 

structure  has  a  single  substrate  with  a  relative  dielectric  constant  of 

e  .  If  h  is  increased  until  it  is  equal  to  h  , ,  then  h  is 
rL2  Ll  total  L2 

equal  to  zero  and  the  structure  once  again  has  only  a  single  substrate, 
this  time  with  a  relative  dielectric  constant  of  c 

rl 

In  Fig.  4.42  the  upper  substrate  is  chosen  to  be  beryllium  oxide 

and  the  lower  substrate  to  be  either  air  or  RT/duroid.  The  effective 

dielectric  constants  for  the  isolated,  even  mode,  and  odd  mode  are 

plotted  as  h  varies  from  zero  to  h  ,  =  0.6  mm.  When  e  is  equal 
Ll  total  rLl  ^ 

to  unity,  the  places  where  the  even  and  odd  modes  have  the  same  c 

reff 

referred  to  as  the  zero  coupling  points,  occur  at  h  -  0.1  mm  and 

h  =0.6  mm  =  h  The  second  zero  coupling  point  occurs  when 

Ll  total  ^ 

h  =  h  ,  because  this  structure  is  a  pure  TEM  line  since  the 
Ll  total  ^ 

dielectric  is  homogeneous  with  a  relative  dielectric  constant  of  unity. 

This  structure  is  just  two  conductors  suspended  in  air  above  a  ground 

plane.  Therefore,  whenever  the  lower  substrate  is  air,  then  one  of  the 

possible  zero  coupling  points  will  be  at  h  =  h 

Ll  total 
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Fig.  4.42.  e  ..  vs. 
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h  for  open  symmetric  coupled  multilayer  microstrips 
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on  beryllium  oxide,  with  RT/duroid  5000  and  air  below. 
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When  the  lower  substrate  is  changed  to  RT/duroid,  then  there  is  no 

value  of  h  that  makes  e  the  same  for  the  even  and  odd  modes. 

LI  reff 

However,  if  the  dielectric  constant  of  the  lower  substrate  is  reduced 

slightly  from  2.2  to  1.8,  then  zero  coupling  points  will  exist. 

Although  there  is  no  geometry  that  allows  the  even/odd  mode 

equalization  when  e  =2.2,  it  is  possible  to  significantly  minimize 

the  difference  in  the  modal  e  and  therefore  drastically  reduce  the 

reff 

effects  of  coupling  distortion.  A  relatively  good  match  in  c  and  c 

re  ro 

is  obtained  when  h  is  in  between  0.32  and  0.42  mm. 

LI 

The  next  two-substrate  structure  has  a  layer  of  gallium  arsenide 

on  top  of  either  an  RT/duroid  or  a  beryllium  oxide  lower  substrate. 

When  the  lower  substrate  is  RT/duroid,  there  are  two  zero  coupling 

points,  one  near  h  =  0.14  mm  and  another  near  h  =0.54  mm.  The 

existence  of  the  zero  coupling  points  depends  on  having  a  configuration 

where  e  is  less  than  e  and  since  at  each  of  the  dimensional 
re  ro 

boundaries,  i.e.  at  h  equal  zero  and  h  equal  h  ,  ,  e  is  greater 
'  Ll  M  Ll  M  total  re  a 

than  c  ,  there  must  be  at  least  two  points  where  e  =  e  The  other 

ro  re  ro 

possibility  is  that  is  greater  than  e  for  all  h^  except  for  one 

point  where  they  are  equal.  This  is  similar  to  quadratic  equations 
which  always  have  two  roots,  although  it  is  possible  that  the  two  roots 
will  have  the  same  value. 

As  the  relative  dielectric  constant  of  the  lower  substrate  is 
increased,  the  locations  of  the  zero  coupling  points  move  closer 
together.  At  some  value  of  e  the  two  points  will  have  moved 
directly  on  top  of  each  other,  giving  only  one  physical  zero  coupling 
structure.  If  e  is  increased  beyond  this  point,  then  the  graphs  for 
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the  even  and  the  odd  modes  no  longer  cross  and  there  are  no  zero 

coupling  points.  This  case  is  shown  in  Fig.  4.43  for  c  ^  =  6.8.  Thus 

with  the  gallium  arsenide  and  beryllium  oxide  combination  it  is  not 

possible  to  make  c  equal  to  c  nor  is  there  any  value  of  h  that 

re  ro  li 

significantly  minimizes  the  differences  between  c  and  c 

re  ro 

Note  that  the  condition  on  e  and  e  for  the  existence  of  the 

rLl  rL2 

zero  coupling  points  is  not  dependent  only  on  the  difference  between 

the  dielectric  constants.  A  combination  of  c  =1.8  and  c  =6.8, 

rLl  rL2 

which  has  zero  coupling  points,  has  a  smaller  change  in  the  relative 

dielectric  constants,  Ae^_  =  5.0,  than  the  gallium  arsenide  and 

beryllium  oxide  combination,  which  has  a  Ae  of  5.4.  Nor  is  their 

existence  dependent  solely  on  the  ratio  of  the  relative  dielectric 

constants.  A  structure  with  e  equal  to  unity  and  e  equal  to  1.7 

rLl  rL2 

has  zero  coupling  points  even  though  the  ratio  of  the  relative 
dielectric  constants,  1.7,  is  less  than  the  gallium  arsenide/beryllium 
oxide  ratio,  which  is  1.79. 

One  reason  that  there  is  not  a  simple  requirement  on  the 
dielectric  constants  of  the  substrates  is  that  the  existence  of  zero 
coupling  points,  especially  at  lower  frequencies,  will  depend  on  the 
parameters  of  the  superstrates  as  well  as  the  substrates.  For 
simplicity  only  open  microstrips  are  considered,  but  the  presence  of  a 
cover  sheet  or  additional  superstrates  could  also  be  treated  in  a 
similar  manner  to  find  the  zero  coupling  points,  if  they  exist,  for  any 
structure.  Clearly  the  complexity  of  the  microstrip  boundary  value 
problem  precludes  the  possibility  of  finding  a  simple  formula  that 
determines  if  the  even/odd  mode  equalization  is  possible.  Indeed,  none 
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of  the  approximate  formulas  even  predict  the  existence  of  these  zero 

coupling  structures,  let  alone  being  able  to  determine  the  dimensional 

requirements.  Therefore  the  analysis  of  these  types  of  structures  must 

be  carried  out  with  a  rigorous  solution  of  the  boundary  value  problem. 

Figures  4.42  and  4.43  show  that  zero  coupling  structures  can  exist 

at  a  single  frequency.  The  next  step  is  to  examine  the  frequency 

dependence  of  the  zero  coupling  points.  In  Figures  4.44  through  4.46, 

the  locations  of  the  zero  coupling  points,  i.e.  the  value  of  h  that 

makes  e  **  e  ,  are  plotted  versus  frequency,  with  the  frequency  on  a 

logarithmic  scale.  As  the  operating  frequency  is  increased,  e  and 

begin  to  approach  the  isolated  line  value  and  hence  the  difference 

between  them  decreases.  This  occurs  for  all  values  of  h  since  the 

Ll 

electrical  separation  of  the  conductors  increases  as  frequency 
increases . 

As  the  frequency  to  continues  increase,  a  point  is  reached  where 

it  is  not  possible  to  uniquely  determine  a  zero  coupling  point  because 

c  e  for  all  h  values.  Beyond  this  frequency,  any  value  of  h 

re  ro  Ll  J  Ll 

will  give  a  good  match  for  e  and  e  .  The  frequency  at  which  this 

re  ro 

occurs  depends  on  how  close  the  modal  effective  dielectric  constants 

need  to  be  for  the  particular  design.  A  larger  tolerance  in  the 

separation  of  the  modal  effective  dielectric  constants  decreases  this 

frequency  and  a  smaller  tolerance  increases  it.  Thus  the  graphs  in 

Figs.  4.44  through  4.46  are  given  only  up  to  a  certain  frequency,  after 

which  almost  any  value  of  h  results  in  less  than  a  tenth  of  a  percent 

separation  in  c  and  c 

re  ro. 
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The  first  zero  coupling  structure,  whose  results  are  shown  in  Fig. 

4.44,  uses  beryllium  oxide  as  the  upper  substrate  and  air  as  the  lower 

substrate.  RT/duroid  could  not  be  used  for  the  lower  substrate  in  this 

case  because,  as  shown  in  Fig.  4.42,  there  is  no  value  of  h  which 

makes  c  equal  to  c  .  Because  the  lower  substrate  is  air,  then  one 
re  ro 

of  the  possible  solutions  is  h  =  h  ^  , ,  and  since  this  represents  a 

Li  total 

pure  TEM  structure,  the  solution  is  constant  for  all  frequencies. 

While  this  solution  is  not  generally  practical  in  terms  of  actual 

circuit  design,  the  second  root  does  represent  a  possible 

configuration,  with  h  =  0.11  mm.  The  zero  coupling  point  is 

constant  versus  frequency  in  the  quasi-static  region  which  ends  near  1 

GHz.  Above  this  frequency,  the  h  value  rises  slightly,  reaching  a 

maximum  around  17  GHz,  that  is  10  percent  higher  than  the  low  frequency 

value.  As  the  frequency  increases  further,  the  second  root  begins 

decreasing  slowly  until  the  root  location  can  no  longer  be  determined 

uniquely.  Although  there  is  a  fairly  rapid  change  in  the  root  location 

from  17  to  100  GHz,  at  these  higher  frequencies  the  separation  between 

e  and  e  is  a  relatively  weak  function  of  h  .  Thus  in  this  region, 
re  ro  Ll 

increasing  the  frequency  decreases  the  importance  of  the  location  of 
the  root. 

The  next  structure,  shown  in  Fig.  4.45,  uses  alumina  as  the  upper 
substrate  with  either  air  or  RT/duroid  below.  Increasing  the  relative 
dielectric  constant  of  the  upper  substrate  from  6.8  to  9.7  decreases 
the  second  zero  coupling  point  for  the  air  substrate  from  0.11  to 
0.068  mm.  The  required  h  achieves  a  maximum  near  17  GHz,  as  before, 
but  with  a  14  percent  deviation  from  the  low  frequency  location  of  the 


143 


root.  However,  the  change  is  only  1.6  percent  of  the  total  height  of 

the  substrate  compared  with  a  2  percent  change  when  beryllium  oxide  was 

the  upper  substrate.  Thus  increasing  the  relative  dielectric  constant 

decreases  the  required  height  of  the  lower  substrate  and  stabilizes  its 

value  with  respect  to  changes  in  frequency. 

When  the  dielectric  constant  of  the  lower  substrate  is  increased, 

the  required  values  move  closer  together  and  change  more  rapidly  with 

frequency  outside  the  quasi-static  region.  In  the  low  frequency 

region,  the  first  root  decreases  from  0.6  to  0.51  mm  while  the  second 

root  increases  from  0.068  to  0.2  mm.  As  c  is  increased  further, 

rLl 

then  the  root  locations  would  move  closer  and  closer  to  each  other 

until  they  reach  the  same  value.  If  e  is  increased  past  this  point, 

then  the  zero  coupling  points  would  no  longer  exist. 

As  the  frequency  is  increased  above  the  quasi-static  region,  then 

the  zero  coupling  points  begin  to  change,  with  the  first  root 

decreasing  slightly  and  the  second  root  increasing  in  value.  The 

minimum  of  the  first  occurs  near  17  GHz  and  is  a  decrease  of  3.4 

percent  of  h  ,  while  the  second  root  reaches  a  maximum  at  15  GHz 
total 

that  is  5.7  percent  of  •  Since  the  changes  in  the  required  h^ 

are  relatively  small,  even  with  the  duroid  substrate,  using  the 
substrates  instead  of  the  superstrates  to  equalize  the  modal  phase 
velocities  gives  very  wideband  results. 

The  final  structure,  shown  in  Fig.  4.46,  replaces  the  alumina 
layer  with  a  layer  of  gallium  arsenide  and  uses  the  same  two  materials 
for  the  lower  substrates.  With  the  air  as  the  lower  substrate,  the 
quasi-static  value  for  the  roots  are  at  0.6  and  0.055  mm  while  the 
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duroid  lower  substrate  give  values  of  0.54  and  0.15  mm.  Additionally, 
the  deviations  in  the  zero  coupling  points  versus  frequency  have  been 
reduced  from  the  previous  configuration.  With  air  as  the  lower 
substrate,  the  maximum  change  in  h  is  1.4  percent  of  h  ,  and  with 
the  duroid  substrate  the  first  root  varies  by  1.4  percent  and  the 
second  by  3.8  percent.  Again  these  very  small  deviations  give  this 
configuration  a  very  wide  bandwidth. 

Realistic  circuit  designs  require  many  different  center  conductor 
spacings,  and  so  it  is  important  to  consider  the  effects  that  different 
spacings  have  on  the  location  of  the  zero  coupling  points.  Figures 
4.47  and  4.48  show  the  zero  coupling  points  for  two  configurations  with 
RT/duroid  as  the  lower  substrate  and  either  alumina  or  gallium  arsenide 
as  the  upper  substrate.  Both  roots  are  computed  for  three  different 
center  conductor  spacings;  0.6,  1.2,  and  1.5  mm. 

As  the  spacing  increases,  there  are  two  major  effects  on  the  zero 

coupling  points;  1)  the  frequency  above  which  the  roots  cannot  be 

uniquely  defined  decreases  and  2)  the  zero  coupling  points  vary  more 

widely  with  frequency.  As  the  physical  spacing  of  the  center 

conductors  increases,  the  frequency  at  which  the  electrical  spacing  is 

sufficiently  large  to  isolate  the  conductors  decreases,  causing  the 

first  effect.  The  second  effect,  which  might  seriously  reduce  the 

bandwidth  of  the  design,  is  offset  by  the  decreasing  importance  of  the 

root  location  with  frequency  in  the  regions  where  it  is  changing 

rapidly.  In  addition,  larger  spacings  have  less  separation  between  the 

even  and  odd  mode  e  and  so  it  becomes  less  critical  to  obtain  and 

ref  f 

exact  match. 
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Fig.  4.47.  Zero  coupling  points  vs.  frequency  for  coupled  microstrips  with 
alumina  and  RT/duroid  5880  substrate,  different  spacings. 
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Using  alumina  as  the  upper  substrate,  as  shown  in  Fig.  4.47, 

results  in  a  larger  change  in  h  as  both  the  spacing  and  the  frequency 

Ll 

change  as  compared  to  using  gallium  arsenide  as  in  Fig.  4.48.  Changing 

the  center  conductor  spacing  from  0.6  to  1.2  mm  gives  an  8  percent 

decrease  in  the  first  root  and  an  11  percent  increase  in  the  second 

with  the  alumina  substrate,  and  a  5  percent  decrease  in  the  first  root 

and  a  7  percent  increase  in  the  second  root  when  the  upper  substrate  is 

gallium  arsenide.  Increasing  the  spacing  on  the  alumina  substrate  to 

1.5  mm  gives  14  and  18  percent  changes  in  the  first  and  second  roots, 

respectively,  compared  with  8  and  11  percent  changes  for  the  gallium 

arsenide  substrate.  The  zero  coupling  points  on  the  alumina  substrate 

also  vary  much  more  with  frequency  than  those  for  the  gallium  arsenide. 

For  both  upper  substrates,  data  at  the  higher  frequencies,  i.e.  around 

80  to  100  GHz,  are  in  the  region  where  e  “  c  for  wide  ranges  of 

re  ro 

h^  ^ ,  and  so  the  rapid  changes  on  the  roots  in  this  region  are  not  as 
relevant . 

The  equalization  of  the  even  and  odd  mode  phase  velocities  has 
many  single  frequency  advantages;  for  example  it  increases  the 
directivity  and  hence  the  isolation  of  microstrip  couplers.  However, 
the  ability  to  achieve  a  good  wideband  match  of  the  modal  phase 
velocities  also  offers  the  possibility  of  eliminating  almost  completely 
the  distortion  of  pulses  due  to  coupling.  Since  it  is  impossible  to 
get  a  perfect  even/odd  match  at  all  frequencies,  the  frequency  spectrum 
of  the  pulse  must  be  considered  in  order  to  determine  which  band  of 
frequencies  is  the  most  important  in  determining  the  h^ ^  value.  If  the 
pulse  spectrum  has  no  significant  frequency  components  outside  the 
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quasi-static  region,  then  the  low  frequency  value  of  h^ ^  could  be  used. 

However,  if  the  pulse  spectrum  is  very  wide,  as  is  the  case  with  the 

pulses  that  have  been  considered  so  far,  then  a  value  of  h  must  be 

Ll 

chosen  that  minimizes  the  difference  between  c  and  c  over  the 

re  ro 

bandwidth  of  interest.  In  general,  this  value  of  h  will  be  somewhere 

between  the  low  frequency  value  and  the  minimum  value,  if  the  first 

root  is  used,  or  the  maximum  value,  if  the  second  root  is  used. 

To  illustrate  how  the  zero  coupling  points  can  be  used  to 

eliminate  distortion  due  to  even/odd  mode  coupling,  the  single  layer 

open  symmetric  coupled  microstrip  on  alumina  substrate  in  Fig.  4.10  is 

redesigned  to  lower  the  coupling  between  the  lines.  The  center 

conductor  widths  and  spacing  are  retained  and  the  total  height  of  the 

substrate  is  maintained  as  well,  but  a  layer  of  duroid  is  added  below 

the  alumina  to  equalized  the  even  and  odd  mode  phase  velocities.  This 

example  will  use  the  second  root,  as  shown  in  Fig.  4.45,  although  it 

would  be  equally  valid  to  use  the  first  root.  The  low  frequency  value 

of  h  is  0.1991  mm  and  the  maximum  value,  which  is  at  15  GHz,  is 
Ll 

0.2334  mm.  The  Gaussian  pulse  has  a  half  width,  half  maximum  of  30 
picoseconds,  and  at  15  GHz,  the  spectrum  has  fallen  off  to  13.5  percent 
of  the  maximum  value.  Also,  since  the  spectrum  of  the  pulse  falls  off 
so  rapidly,  equation  (4.1),  the  lower  frequencies  are  the  most 
important  and  frequencies  higher  than  15  GHz  do  not  critically  affect 
the  pulse.  Using  this  information,  h  is  chosen  to  be  0.22  mm. 

Ll 

Using  the  chosen  dimensions,  c  __  is  computed  and  is  shown  in 

ref  f 

Fig.  4.49  as  a  function  of  frequency  both  as  a  preliminary  to  the 
computation  of  pulse  distortion  as  well  as  to  check  the  validity  of  the 
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4.49.  e  vs.  frequency  for  low  coupling  structure  using  alumina  and 
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RT/duroid  5880  substrates. 
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matching  process.  In  the  quasi-static  region  as  well  as  for 

frequencies  above  25  GHz,  e  is  greater  than  e  ,  which  is  a  reversal 

ro  re 

of  the  single  layer  situation,  while  from  5  to  25  GHz,  e  becomes 

slightly  larger  than  This  design  achieves  a  very  good  match  of 

the  even  and  odd  mode  e  at  the  lower  frequencies,  with  less  than  a 

one  percent  difference  in  the  values  up  to  40  GHz.  The  separation 

between  e  and  e  is  the  largest  in  frequencies  from  40  to  85  GHz  and 
re  ro 

reaches  a  maximum  separation  of  less  than  1.6  percent  at  60  GHz.  When 
the  frequency  gets  larger  than  100  GHz,  the  separation  between  and 

becomes  negligible  due  to  the  increasing  electrical  separation  of 
the  center  conductors.  Since  neither  the  even  nor  the  odd  mode  c 

ref  f 

is  consistently  larger  than  the  other,  the  distortion  of  the  pulses 
will  not  appear  the  same  as  the  previous  cases,  and  it  is  not 
immediately  evident  whether  the  leading  sense  line  response  will  be 
positive  or  negative.  However,  since  the  even  mode  is  faster  in  the 
quasi-static  region,  where  the  majority  of  the  pulse  spectrum  lies,  it 
is  probable  that  the  leading  response  will  be  due  mostly  to  even  mode 
components  and  therefore  the  sense  line  should  have  a  positive  leading 
response . 

In  addition  to  reducing  the  separation  of  c  and  c  ,  the 

re  ro 

addition  of  the  substrate  layer  drastically  changes  the  dispersive 
characteristics  of  the  lines.  Comparing  Fig.  4.49  with  Figs.  4.9  and 
4.10,  shows  that  the  new  structure  will  be  less  dispersive  than  the 
single  layer  alumina  structure,  but  slightly  more  dispersive  than  the 
configuration  with  the  duroid  substrate. 
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Having  computed  the  effective  dielectric  constant  for  the 
frequency  band  of  interest,  the  propagation  of  the  pulse  on  the  coupled 
lines  may  be  computed.  Three  different  distances  are  used,  62.5,  125, 
and  250  mm,  for  comparison  with  the  previous  pulse  distortion  results 
for  both  the  duroid  and  alumina  substrates.  The  62.5  mm  distance  used 
in  Fig.  4.50  corresponds  to  the  graphs  in  Fig.  4.15  where  the  substrate 
is  a  single  layer  of  alumina.  The  125  mm  length  used  for  Fig.  4.51  is 
comparable  to  both  Fig.  4.16,  alumina  substrate,  and  Fig.  4.11,  the 
duroid  substrate.  Finally,  Fig.  4.52  with  a  distance  of  250  mm, 
corresponds  to  the  results  in  Fig.  4.12  with  a  duroid  substrate.  The 
isolated  line  response  is  not  graphed  in  Figs.  4.50  through  4.52 
because  there  is  no  visible  difference  between  it  and  the  signal  line 
response . 

At  62.5  mm  the  signal  line  response  is  very  similar  to  the 
undistorted  pulse,  the  only  difference  being  that  it  is  shifted 
slightly  with  respect  to  it.  The  sense  line  response  is  negligible, 
having  a  maximum  amplitude  of  less  than  one  percent  of  the  undistorted 
pulse.  At  this  distance,  the  lines  on  the  alumina  substrate  showed  a 
significant  reduction  in  the  amplitude  of  the  signal  line,  about  40 
percent,  and  had  a  sense  line  response  that  approached  50  percent  of 
the  undistorted  pulse. 

As  the  pulse  travels  farther  on  the  low  coupling  structure, 
reaching  125  mm,  dispersion  distortion  begins  to  affect  the  pulse, 
while  coupling  distortion  is  unnoticeable .  The  signal  line  is  still 
unaffected  by  coupling,  and  the  maximum  sense  line  response  has  risen 
to  only  1.9  percent  of  the  undistorted  pulse.  Contrasting  this  is  the 
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Fig.  4.50.  Pulse  distortion  on  low  coupling  structure,  1  =  62.5mm, 
r  =  30ps. 


152 


-  Signal  line 

-  Sense  line 

.  Undistorted  / 


CO 


CrUl 

l-0,  hyj 

£  . 

=  9.7,  h 

rL2 

’  L2 

£  . 

=  2.2.  h. 

rLl 

Ll 

w  = 

0.6mm,  s 

0.38mm 

0.22mm 

1.2mm 


155 


responses  shown  in  Figs.  4.11  and  4.16,  where  the  signal  line  response 
on  the  alumina  substrate  had  almost  completely  split  into  the  even  and 
odd  mode  components  and  on  the  duroid  substrate  it  had  suffered  a  20 
percent  degradation  in  amplitude. 

At  the  final  distance,  250  mm,  dispersion  distortion  has 
noticeably  changed  the  signal  line  response,  but  coupling  distortion 
has  been  eliminated  as  a  limiting  factor  for  the  structure.  Even  at 
this  distance,  there  is  no  visible  difference  between  the  graphs  of  the 
signal  line  and  the  isolated  case,  indicating  that  coupling  has  not 
affected  the  intended  signal.  The  sense  line  shows  a  maximum  response 
that  is  just  a  little  over  four  percent  of  the  undistorted  pulse,  which 
indicates  that  there  is  a  negligible  amount  of  crosstalk  between  the 
lines.  The  previous  results  for  the  duroid  substrate,  shown  in  Fig. 
4.12,  indicated  that  coupling  was  by  far  the  dominant  mechanism, 
producing  a  very  large  spurious  response  and  seriously  degrading  the 
signal  line  pulse. 

Thus  through  the  use  of  the  zero  coupling  points,  even/odd  mode 
coupling  was  essentially  eliminated  as  a  distortion  mechanism  on 
microstrip  lines  even  though  they  were  in  close  proximity.  Using  the 
correct  combination  of  the  two  substrate  materials  enabled  the  design 
of  a  structure  that  exhibited  better  isolation  between  the  lines  than 
could  be  obtained  by  using  either  of  the  substrate  materials  alone.  In 
spite  of  the  extreme  distances  traveled,  crosstalk  was  held  to 
negligible  levels  and  the  intended  signal  suffered  no  noticeable 


distortion  due  to  coupling. 


CHAPTER  5 
CONCLUSIONS 


Solving  the  microstrip  boundary  value  problem  in  a  rigorous  manner 
gave  a  simplified  representation  of  multilayer  coupled  line  structures 
and  enabled  the  accurate  computation  of  the  line  parameters  as  well  as 
pulse  propagation  for  these  structures.  First,  a  generalized  Green's 
function  was  presented  as  a  simple  recurrence  relation,  that  is 
applicable  for  microstrip  structure  with  any  number  of  superstrates  or 
substrates.  Next  the  Fourier  transform  and  the  even/odd  mode  approach 
were  used  to  obtain  the  formulas  necessary  to  analyze  the  transient 
response  on  coupled  transmission  line,  including  the  effects  of 
dispersion,  losses  and  coupling.  Numerical  results  were  presented  for 
both  single  and  multilayer  structures,  showing  that  a  full  wave 
analysis  in  necessary  for  accurate  results.  Finally,  a  new  method  for 
equalizing  the  even  and  odd  mode  phase  velocities  through  the  control 
of  the  electrical  characteristics  of  the  substrates  was  presented. 

The  rigorous  solution  of  the  multilayer,  multiconductor  microstrip 
boundary  value  problem  with  the  spectral  domain  approach  resulted  in  a 

relatively  simple  formulation  for  the  microstrip  parameters.  The 

.  •  ,  y  y 

derivation  began  by  choosing  TE  and  TM  modes,  allowing  the  problem  to 

be  solved  separately  for  each  mode,  since  the  modes  satisfy  the 

boundary  conditions  independently  in  the  spectral  domain.  Next,  the 

SDA  was  used  to  obtain  the  Green's  function,  initially  for  a  simple 
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single  layer  shielded  structure.  To  illustrate  the  use  of  the  SDA,  a 
six-layer  structure  was  considered  and  the  derivation,  resulted  in  a 
simple  recurrence  relation  for  the  Green's  function  of  multilayer 
configurations.  The  recurrence  relation  allows  the  Green's  function  to 
be  easily  computed  with  a  single  subroutine  for  any  arbitrary 
multilayer  structure,  rather  than  requiring  a  different  subroutine  for 
each  possible  configuration.  Finally,  the  Green's  function  was  also 
expressed  in  a  form  suitable  for  the  calculation  of  the  parameters  of 
slot  line  and  coplanar  waveguide  type  structures. 

The  even/odd  mode  analysis  and  the  frequency  domain  approach  were 
used  to  derive  the  formulas  for  the  propagation  of  finite  time  domain 
signals  on  coupled  transmission  lines.  Using  the  frequency  domain 
approach  simplified  the  transient  analysis  of  the  coupled  line  system, 
allowing  fast,  accurate  computation  of  the  output  waveforms.  The 
even/odd  mode  approach  made  it  possible  to  split  the  symmetric 
four-port  system  into  two  simpler  two-port  networks.  Combining  the 
even/odd  mode  analysis  with  the  frequency  domain  approach  gave  simple 
expressions  for  the  time  domain  response  of  lossy  and  dispersive 
coupled  transmission  lines.  These  expressions  clearly  illustrated  how 
losses,  coupling,  and  dispersion  degrade  the  input  signal  and  how 
differences  in  the  complex  propagation  constants  for  the  even  and  odd 
modes  produce  a  spurious  response  on  the  sense  line. 

Numerical  results  for  the  frequency  dependent  effective  dielectric 
constant  and  pulse  dispersion  for  a  single  layer  symmetric  coupled 
microstrip  were  presented,  for  different  substrate  materials,  center 
conductor  spacings,  and  distances.  Coupling  between  the  lines 
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increased  when  the  center  conductor  spacing  was  increased,  the  distance 
increased,  and/or  the  dielectric  constant  of  the  substrate  increased. 
Examining  the  distortion  of  pulses  at  different  distances  showed  that 
coupling  distortion  separates  the  input  signal  into  even  and  odd  mode 
pulses  as  the  distance  gets  very  large.  By  studying  the  effects  of 
center  conductor  spacing,  it  was  demonstrated  that  the  conditions  for 
uncoupled  operation  of  transmission  lines  depended  not  only  on  the 
center  conductor  spacing  or  substrate  height,  but  also  on  the  substrate 
materials,  the  distance  over  which  the  lines  are  coupled,  as  well  as 
the  other  dimensions  of  the  structures. 

The  generalized  Green's  function  was  used  with  the  SDA  to  consider 

the  effects  of  multiple  substrates  and  superstrates  on  the 

characteristics  of  coupled  lines.  Adding  even  a  thin  lower  substrate 

layer  with  a  different  dielectric  constant  caused  significant  changes 

in  the  c  of  the  isolated  case  as  well  as  for  the  even  and  odd 

ref  f 

modes.  Placing  a  superstrate  layer  with  the  same  dielectric  constant 

as  the  substrate  on  top  of  the  center  conductor  interface  also  had  an 

important  effect  on  the  microstrip  parameters.  By  adjusting  the  height 

of  this  superstrate  layer,  it  was  possible  to  equalize  the  even/odd 

phase  velocities  for  the  structure,  although  the  equalization  was  good 

only  in  the  quasi-static  region  or  over  a  very  narrow  band  of 

frequencies.  When  a  cover  sheet  was  introduced  to  the  two  superstrate 

structure,  e  ^  changed  over  a  wide  range  of  values  as  the  height  of 

the  cover  layer  was  varied.  Finally,  it  was  shown  that  the  relative 

positions  of  the  substrate  layers  changes  not  only  the  value  of  e  ^ 

for  the  isolated  case,  but  can  cause  the  even  and  odd  mode  e  __  to 

ref  f 
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switch  in  relation  to  each  other. 

A  new  method  of  equalizing  the  even  and  odd  mode  phase  velocities 
was  presented  that  could  obtain  a  match  over  a  very  wide  band  of 
frequencies.  To  overcome  the  difficulties  encountered  from  using  the 
superstrate  equalization  scheme,  the  modal  phase  velocities  were 
equalized  by  choosing  substrates  with  the  appropriate  dimensions  and 
materials.  Since  for  all  frequencies  most  of  the  fields  lie  in  the 
substrate  layers,  this  matching  procedure  gives  good  results  not  only 
in  the  quasi-static  region,  but  also  for  frequencies  well  above  it. 

The  appropriate  substrate  heights  for  even/odd  mode  equalization, 
called  zero  coupling  points,  were  presented  as  functions  of  frequency 
and  center  conductor  spacing  for  different  substrate  materials.  Using 
these  zero  coupling  points,  a  structure  was  designed  that  exhibited 
better  isolation  between  the  adjacent  lines  than  could  be  obtained  by 
using  either  of  the  substrate  materials  alone.  When  pulse  distortion 
was  considered  for  this  structure,  it  was  shown  that  the  crosstalk  was 
almost  non-existent  and  the  input  signal  had  suffered  no  distortion  due 


to  coupling. 


CHAPTER  6 

RECOMMENDATIONS 


Future  research  in  this  area  would  be  concerned  with  increasing 

the  speed  and  accuracy  of  the  computations  and  considering  both  losses 

and  discontinuities  in  a  rigorous  manner.  To  increase  the  accuracy  of 

the  Spectral  Domain  Approach,  more  current  expansion  functions  should 

be  used.  It  is  very  important  to  choose  'good'  expansion  functions 

since  a  poor  choice  will  necessitate  a  very  large  matrix  which  may  not 

even  converge.  Since  the  SDA  is  used  in  a  very  general  manner,  it  is 

possible  to  include  the  effects  of  lossy  dielectrics  by  using  a  complex 

dielectric  constant.  This  substitution  would  make  the  Green's  function 

complex,  resulting  in  two  equations  (the  real  and  imaginary  parts)  to 

be  solved  for  two  unknowns,  a  and  6  .  The  results  from  this  research 

z  z 

could  be  used  to  verify  approximate  formulas  used  for  dielectric  losses 
and  to  specify  the  range  of  validity  of  the  approximations. 

It  was  shown  in  Chapter  3  that  differences  in  the  even  and  odd 
mode  attenuation  constants  cause  coupling  and  distortion  on  coupled 
lines.  Currently  there  are  no  approximate  formulas  or  results  from 
rigorous  solutions  available  for  the  modal  attenuation  constants.  The 
attenuation  constants  for  the  even  and  odd  modes  could  be  determined  by 
using  the  SDA  for  dielectric  losses  and  by  modifying  approximate 
formulations  for  the  conductor  losses.  Once  the  modal  attenuation 
constants  have  been  determined  as  a  function  of  frequency,  pulse 
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distortion  on  coupled  lines  due  to  attenuation  distortion  and 
attenuation  coupling  could  be  investigated. 

To  accurately  model  complete  MMIC  circuits,  the  effects  of 
discontinuities  must  be  considered  in  a  rigorous  manner.  The  SDA  can 
be  used  to  consider  these  problems  since  the  derivation  transformed  in 
both  the  z  and  x  directions.  By  using  the  appropriate  current 
expansion  functions,  which  now  must  be  functions  of  (3^  as  well  as 
it  will  be  possible  to  characterize  discontinuities  in  the  z  direction 
of  the  structure.  Discontinuities  that  might  be  considered  include: 
abrupt  changes  in  the  center  conductor  width,  gaps  in  the  center 
conductor,  bends,  crossing  lines,  etc.  These  results  could  then  be 
applied  to  consider  the  transient  response  of  very  complex  circuits. 
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GREEN’S  FUNCTION  FOR  GENERALIZED  PLANAR  STRUCTURE 


Efficient  solution  of  printed  circuit  transmission  lines  with  the 

Spectral  Domain  Approach  (SDA)  requires  a  closed  form  expression  for 
,  ,  ~TE  ~TM 

the  modal  input  impedances,  Z  and  Z  .  The  solution  begins  by 

solving  the  boundary  value  problem  for  a  six  layer  structure,  first  for 
y  V 

the  TM  mode  and  then  for  the  TE  mode.  This  derivation  shows  a 
pattern  that  leads  to  the  recurrence  relation  presented  in  section 


2.3.8. 

y  y  . 

Since  TE  and  TM  modes  are  being  used,  the  boundary  conditions 

can  be  enforced  separately  for  each  mode.  In  order  to  simplify  the 

notation,  a  different  geometry  is  used  for  the  derivation,  shown  in 

figure  A.l.  The  structure  is  surrounded  perfect  electric  conductors  on 

all  four  sides,  i.e.  at  x  =  a,  x  =  -a,  y  =  h  ,  y  =  h  +h  +h  +h  +h  +h  . 

Non-shielded  structures  can  be  considered  by  letting  a  — >  oo,  h^  — >  oo, 

and/or  h  — >  oo  as  is  appropriate  for  the  structure.  For  clarity,  only 
6 

one  center  conductor  is  shown,  but  any  finite  number  could  be 

.  y  . 

considered.  Beginning  with  the  TM  mode  first,  the  appropriate  vector 

potentials,  A  (0  ,  y,0  ),  for  each  of  the  six  regions  are  chosen  as 
yi  x  z 

follows : 

REGION  1 :  (0  <  y  <  h  ) 


<V  A/ 


A  =B  (0,0) cosh (a  y) 
yl  1  x  z  yl 

REGION  2:  (h  £  y  £  d  ,  d  =h+h) 

A  =B  (0  ,0  ) cosh [a  (d  -y) ] +B  (0  ,0  )sinh[a  (d  -y) ] 
y2  2A  x  z  y2  2  2B  x  z  y2  2 


(A. la) 


(A. lb) 
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REGION  3:  (d  s  y  £  d  ,  d  =  h  +h  +h  ) 

-  2  3  3  1  2  3 

A  =B  {p  ,p  )cosh[a  (d  -y)]+B  (P  ,/3  )sinh[a  (d  -y)  ]  (A.lc) 

y3  3A  '  x"  z  y3  3  *  3B  1  x"z  y3  3 

REGION  4:  (d  s  y  s  d  ,  d  =  h  +h  +h  +h  ) 

-  3  4  4  1  2  3  4 

A  =B  (p  ,P  ) cosh [a  (y-d  )]+B  (p  ,P  )sinh[a  (y-d  ) ]  (A. Id) 

y4  4A  x  z  y4  3  4B  x  z  y4  J  3 

REGION  5:  (d  ^  y  <  d  ,  d  =  h  +h  +h  +h  +h  ) 

-  4  5  5  1234  5 

A  =B  (p  ,P  ) cosh [a  (y-d  )]+B  {p  ,p  )sinh[a  (y-d  ) ]  (A.le) 

y5  5A  X  z  y5  4  5B  x  z  y5  4 

REGION  6:  (d  £  y  <  d  ,  d  =  h  +h  +h  +h  +h  +h  ) 

-  5  6  6  12345  6 

A  -B  AP  ,  P  )  cosh  [a  (d  -y)  ]  (A. If) 

y6  6  x  z  y6  6 

Where  the  B  (S  ,p  )'s  are  unknown  functions  of  the  transform 
n  x  z 

variables,  which  are  to  be  determined  from  the  boundary  conditions. 

They  are  abbreviated  as  B  throughout.  A  and  A  are  chosen  based  on; 

n  yl  y6 

1)  the  tangential  E  fields  vanishing  at  y  =  0  (at  ground  plane)  and  2) 

the  tangential  E  fields  vanishing  at  y  =  d  (at  the  cover  sheet) .  If 

6 

the  structure  doesn't  have  a  cover  sheet,  i.e.  it  is  an  open  structure, 

then  d  - >  oo  because  h  - >  oo  and  the  fields  at  infinity  must  also 

6  6 

,  ,  y 

vanish  due  to  the  outward  radiation  condition.  Because  the  TM  mode 
allows  either  the  x  or  z  directed  fields  to  be  used  to  enforce  the 
boundary  conditions,  only  the  x  directed  electric  and  magnetic  fields 
will  be  given.  The  six  vector  potentials  in  equations  (A.la-f)  are 
used  with  equations  (2.12a-f)  to  calculate  the  fields  in  each  region. 


REGION  1 : 


P  a 

tm  x  yl  ,  ,  , 

5  =  -  B  - —  sinh(a  y) 

xl  l  u)i  e  yl 


tm  .  z 

H  =  i  B  —  cosh (a  y) 
xl  1  yi 

2  2  2  2  2 

P  =  P  +  P  -  a  =  U  H  C 
*1  '  X  '  z  yl  Ml  1 


(A. 2a) 


(A. 2b) 


(A. 3) 
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REGION  2: 


„TM 

JX2 


TM 

H 

x2 

REGION  3: 


„TM 


TM 


P  a 

X  ' 

Wfi 


1  a  r  "l 

X  y2  B  sinh [a  (d  -y)]+B„  cosh[a  (d  -y) ] 

C  [  2A  y2  2  2B  y2  2  J 


3  ~  |B2Acosh[ay2(d2-y)]+B2BSinh[ay2(d2'y)1 

2  2  2  2 
p  +  p  -  a  „  =  0)  |!  e 
Mx  1  z  y2  2  2 


Pa  r 

JLVi  , 

*3  *3  L 

^3  L 


B  sinh [a  (d  -y)]+B,  cosh[a  (d  -y) ] 
W/J3C3  [  3A  y3  3  3B  y3  3 


H  =  j 
X3  J 


] 


B^cosh [a^3 (d3~y) ] +B3Bsinh  tay3 <d3~y)  J 


REGION  4: 


„TM 

"x4 


_2  _2  2  2 

8  +  P  -  a  =  w  u  e„ 

'x  1  z  y3  j  3 


P  a  f  "1 

-  ■  x  y^  B  sinh[a  (y-d  ) ] +B  cosh[a  (y-d  )  ] 

W(i^C^  [  4A  y 4  3  4B  y4  3  J 

p  a  r  1 

_  y -  B  sinh [a  (y-d)]+B  cosh[a  (y-d  )  ] 

U)1  E  [  4A  y4  3  4  B  y4  3J 


,TM 

Jz4 


M 

X4 


K- 

REGION  5: 


=  j  —  [b  cosh  [a  (y-d  )]+B  sinh[a  (y-d  )]1 
J  (1^  [  4A  y  4  3  4  B  y4  3  J 

„2  _2  2  2 

8  +  8  -  a  =  «  u  e 

x  z  y4  *4  4 


P 


,TH 

JX5 

M 

xS 

2 


P  a  f  1 

=  -  X  ys  B  sinh [a  (y-d  )  ] +B  cosh[a  (y-d)] 

5 A  y5  4  5B  y5  4  J 

=  i  —  [b  cosh[a  (y-d  ) ] +Bc  sinh[a  (y-d  )]] 

J  L  5A  y5  4  5 B  y5  4  J 

2  2  2  2 

p  +  P  -  a  =  u  M  e 

x  z  y5  5  5 


REGION  6: 


^TM 

Jx6 


M 
X  6 


P  a 

B  x  sinh [a  (d  -y) ] 
6  0)fl  C  y6  6 

6  6 

0Z 

j  Bg  —  cosh[ay6(d6-y)] 


(A. 4a) 

(A.  4b) 
(A. 5) 

(A. 6a) 

(A. 6b) 
(A. 7) 

(A. 8a) 

(A. 8b) 

(A. 8c) 
(A. 9) 

(A. 10a) 

(A. 10b) 
(A. 11) 

(A. 12a) 

(A. 12b) 
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06  "  +  K  ~  V  =  U\Be  <A-13> 

Since  the  SDA  requires  both  the  x  and  z  directed  electric  fields 

,  TM  TM  , 

at  the  center  conductor  interface,  both  E  and  E  are  included.  Note 

x4  z4 

,  TM  TM 

that  using  E  and  E  could  have  been  used  instead  of  the  electric 
x3  z  3 

fields  in  region  4  because  the  center  conductors  are  assumed  to  be 
infinitely  thin. 

Now  that  the  fields  have  been  determined,  the  ten  independent 
boundary  conditions  (two  have  already  been  applied)  can  be  used  to 
eliminate  the  BN'S-  First  the  tangential  fields  are  matched  at  the 
first  dielectric  interface  from  the  bottom. 


1)  E 


TM 

Xl 


TM 

TM  _T  M 

—  E 

.  or  E  .  =  E 

y=h  x2 

J  1 

y=h  zl  y=h  z2 

r  1  j  j 

p  U.  MU.  r- 

x  yl  ,  ,  .  x  y2 

- — sinh  (a  h  )= - —  1 

i  utl1c1  y1  1  ^2^2  L 

a  r 


y=h. 


B  sinh (a  h  ) +B  cosh (a  h  ) 
2  A  y2  2  2B  y2  2 


u  c  a 
'rl  rl  ■ 

B  =-sech (a  h  ) - 

l  yl  l  U  e 


r2  r2  yl 


2)  H 


TM 

Xl 


TM 

TM  TM 

L  =  H 

,  or  H  =  H 

y=hi  x2 

y=h  zl  y=h  z2 

1  1  J  1 

B  sinh (a  h  ) +B  cosh (a  h  ) 
2  A  y2  2  2B  y2  2 


y=h. 


0Z  13  z  r  l 

:Bx  -  cosh(ayihi)  =  :-|^B2Acosh(ay2h2)+B2Bsinh(ay2h2)j 


Substituting  in  for  B^  using  (A. 15) 


e  coth(a  h  )  a 
rl  yl  ' 


a 


yi 


/  w- 

i _ y2  : 

c 

r2  L 


B  sinh (a  h  ) +B  cosh ( 
2A  y2  2  2B 


a  h  )  = 

y2  2  J 


-  B  cosh (a  h  )  -  B  sinh (a  h  ) 
2  A  y2  2  2B  y2  2 


c  coth (a  h  ) 
rl  yl  1 


yi 


c  c  coth  (a  h  ) 

r2  r2  y2  2 

B  -  +  B  - - - 

2A  a  2b  a 

y2  y2 


c 

—  B 
a  2  a 
y2 


-  B  (c  /a  ) 
2B  r2  y2 


(A. 14) 

(A. 15) 


(A. 16) 


(A. 17) 


e  coth  (a  h  ) 

r2  y2  2 


(A. 18) 
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c 

e  coth (a  h  ) 

r2 

rl  yl  l 

a 

a 

y2 

L  y1 

e  coth(a  h  )"|fe  coth(a  h  )"|  fc 
rl  yl  1  r2  y2  2  r2 


(A. 19) 


In  a  similar  manner,  the  tangential  fields  at  the  first  dielectric 
interface  from  the  top  are  now  matched. 


E  =  E  or  E  .  =  E 

x5  y=d  x 6  y=d  z5  y=d  z6  y=d 


P  a 

X  V6  . 

B  - — sinh (a  h  )  = 

6  W(1  €  y  6  6 

6  6 

£  a 
'  x  ys 


- B  sinh(a  h  ) +B  cosh(a  1 

Wjl  C  I  5A  y5  5  5B  y5 


(A. 20) 


u  e  a 
r6  r6  v5 

B  =-sech  (a  h  ) - — 

6  ye  6  Li  c  a 

r5  r5  y6 


B  sinh (a  h  )+B  cosh (a  h  )  (A. 21) 

5A  y5  5  5B  y5  5 


i  =  H  or  H  =  H 

x5  y=d  x6  y=d  z5  y=d  z6  y=d 


jB  —  cosh  (a  h  ) 
J  6  n  ye  6 

6 


Ar 

=  ’sri1 


B  cosh (a  h  ) +B  sinh (a  h  ) 
5 A  ys  5  5B  y5  5 


(A. 22) 


Substituting  in  for  B  using  (A. 21) 

6 


e  coth(a  h  )  a 
r6  y6  6  y5 


B  sinh  (a  h  ) +B  cosh  (a 
C  5  A  y5  5  SB  y5 


v]- 


-  B  cosh (a  h  )  -  B  sinh (a  h  ) 
5 A  y5  5  5B  y5  5 


(A. 23) 


c  coth (a  h  ) 
r6  ye  e 


c  c  coth  (a  h  ) 

B  -El  +  B  — _ 32_L 

5a  a  5b  a 

ys  y5 


e  e  coth  (a  h  ) 

rS  „  r5  y5  5  „  ,  .  .  2 

-  B  - - -  -  B  (e  /a  ) 

a  5A  a  5b  rs  ys 

ys  y5 


(A. 24) 


c 

c  coth (a  h  ) 

rs 

re  y6  e 

a 

a 

ys 

ye 

c  coth(a  h  )1  Tc  coth  (a  h  )1  fc 
re  ye  6  rs  y5  5  r5 


(A. 25) 
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To  simplify  the  notation,  let 


fTM 

si 


TM 
S  2 


fTM 

S3 


c  coth (a  h  ) 
ri  yl  l 

a 

yi 

e  coth (a  h  ) 


r2 


y2  2 


a 

y2 

c  coth (a  h  ) 
r3 _ y3  3 

a 


y3 


e  coth (a  h  ) 

rTM  _  r4  y4  4 

S  4 


TM 

S5 


TM 
S  6 


a 

y4 

e  coth (a  h  ) 
r5  y5  5 

a 

ys 

c  coth (a  h  ) 
r6  y6  6 

a 

y6 


(A. 26a) 

(A. 26b) 

(A. 26c) 

(A.  26d) 

(A. 26e) 

(A . 2  6f ) 


and 


TM  TM  2 

Y  Y  +(e  /a  ) 
TM  _  Si  S2  r2  y2 

2  yTM  +  yTM 

Si  S2 

Y™  Y™+(c  /a  )2 
TM  =  s  6  s5  r5  y5 

5  yTM  +  yTM 

S  6  S  5 


(A. 27a) 


(A. 27b) 


Using  these  expressions  in  equations  (A. 19)  and  (A. 25)  and  solving  them 
for  and  respectively, 


2  A 


-  B 


a  Y™  Y™+(c  /a  )2 
y2  si  s2  r2  y2 


2  B  C 


r2 


yTM  +  yTM 
sl  S2 


a 

B  Jti  y™ 

2  B  e  2 
r2 


TM  ,  TM  .  ,2 

a  Y  Y  +  c  /a 
y5  S6  S5  r5  y5 

B  =  -  B  — - - -  =  -  B 

5A  SBC 

r5 


TM  TM 

Y  +  Y 
S6  S  5 


a 

_y5  yTM 

5  B  C  5 
r5 


(A. 28a) 


(A. 28b) 


Proceeding  to  the  next  dielectric  interface  up  from  the  1-2 
interface,  the  fields  at  the  2-3  interface  are  forced  to  be  continuous. 


5)  E 


TM 

X2 


TM  TM 

TM 

,  =  E  ,  or  E 

^  =  E 

y=d  x3  y=d  z3 

J  2  2 

y=d  z3 

2 

/3  a  /3  a 
x  y2  _  x  y3 


2 B  W2C2  WM3e3L  3A 


y=d. 


B  sinh (a  h  ) +B  cosh ( 
"*  y3  3  3B 


a  h  ) 
y  3  3  J 


(A. 29) 
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6)  H 


2  B 


TM 


a  a  u  r 

V2  v3  r2  r 

=  -  B  sinh(a  h  ) +B  cosh(a  h  ) 

E  £  LI  3A  V3  3  3B  V3  3 

r2  r3^r3 L  J  J 


x2 


j  B 


2A 


TM 

TM 

TM 

,  =  H 

,  or  H 

,  =  H 

y=d  x3 

2 

y=d  z  2 

2 

y=d  z  3 

*  2 

—  -  j  — fl 

A  M2  L 

t,L 


y=d 


B  cosh (a  h  )+B  sinh(a  h  ) 
3  A  y3  3  3B  y3  3 


B  cosh (a  h  ) +B  sinh  (a 
3  A  y 3  3  3  B  y3 


h3’] 


Substituting  in 

for 

B 

2  A 

a 

-  B  — ^  Y 

TM 

1  T= 
h 

1  PO 

2  B  E 

r2 

2 

"r3 

a 

B  =  - 

2  B  E 

r2 

Mr2 

[b 

L  3  A 

Equating  (A. 34) 

and 

(A. 3 

y3  3  3B 


y3 


v][<1 


-i 


(A. 30) 


(A. 31) 

(A. 32) 


(A. 33) 

(A. 34) 


a  u 
y3M 


-  1 


r3  r3 


B  sinh (a  h  ) +B  cosh (a  h  ) 
3 A  y3  3  3B  y3  3 


r2 

- B  cosh (a  h  ) +B  sinh (a 

Mr3L  3A  y3  3  3B  y3 

TM 


vM1 


Multiplying  both  sides  by  [Y  e  u  ]/[a  1 1  sinh (a  h  )] 

2  r3  r3  y3  r2  y3  3 


(A. 35) 


B  +B  coth  (a  h  )  Y™  =  - 
[  3A  3B  y3  3  J  [  2  J 


B 


3A 


e  coth (a  h  ) 
yTM+  r3  y3  3 

2  a 

y3 


r3 


e  coth (a  h  ) 
r3  y3  3 

B  - - - +B 

3A  a  3B  a 

y3  y3 


-  B  Y‘"coth(a  h  ) +c  /a 
3B  |_  2  y3  3  r3  y3j 


TM 


Using  (A. 26c)  and  factoring  cr3/ay3  out  °f  the  left  hand  side, 


B 

3  A 


TM  TM  2 

a  „  Y™  YT”+(c  J  a  X 
y3  2  S3  r3  y3 


3  B  E 


r3 


yTM  +  yTM 
2  S3 


(A. 36) 


(A. 37) 


(A. 38) 


Since  this  expression  is  very  similar  to  (A. 27a),  let 


Y™  Y™+{e  /a  )2 
rTM  2  S3  r3  y3 

'3  TM  TM 

Y  +  Y 
2  S3 


(A. 39) 
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Then  (A. 38)  simplifies  to; 


B  -  -  B  Y™  JH 
3  A  3B  3  C 

r3 


Again  in  a  very  similar  manner,  the  tangential  fields  at 
dielectric  interface  are  equated. 


(A.  40) 

the  4-5 


7)  E 


TM 

X5 


TM 

_TM 

TM 

=  E 

.  or  E 

—  E 

y=d  x4 

J  4 

y=d  z5 

1  4 

y=d  z  4 

1  4 

y=d 


/3  a  /3  a  1 

B  — X  =  -  — — —  B  sinh  (a  h  ) +B  cosh(a  h  ) 

5B  W|J5C5  4A  y4  4  4B  yl  4  J 

a  a  u  r  -i 

y5  y4rr5 I 

-  =  -  B  sinh(a  h  ) +B  cosh(a  h  ) 

C  C  U  4A  y4  4  4  B  y4  4 


5B 


8)  H 


TM 

X5 


r5  r4  r4 
TM 


TM 


=  H  or  H 

y=d  x4  y=d  z5 

4  4 


=  H™ 
y=d  z  4 

J  4 


y=dA 


j  B 


P„  Prr  T 

-  =  j  - 

A  M5  l 


B  cosh (a  h  ) +B  sinh (a  h  ) 
4  A  y4  4  4  B  y4  4 


] 


B  cosh (a  h  ) +B  sinh (a 
5A  U [  4A  y 4  4  4B  y4 

Substituting  in  for  B  using  (A. 28b) 

5  A 


v] 


a  p 

B  — ^  Y™  =  — ^|B  cosh  (a  h  ) +B  sinh  (a  h  ) 

5B  e  5  li  4A  y4  4  4b  y4  4 

r5  ^r4 1  1 


a 


B 


y5  . 


5B  C 


f’ 

•v.L 


r5  r  4 

Equating  (A. 42)  and  (A. 46) 


B  cosh (a  h  ) +B  sinh (a  h  ) 
4A  y  4  4  4  B  y  4  4 


a  u 
y* 


Ka  L 


r4  r4 


B  sinh (a  h  ) +B  cosh  (a 
4  A  y4  4  4  B  y4 


v]  = 


!!«[; 

'‘r.L 


B  cosh (a  h  ) +B  sinh (a  h  )  Y 
4  A  y4  4  4B  y4  4  5 


TM 


TM 


,-l 


Multiplying  both  sides  by  [Y  c  u  ]  /  [a  u  sinh (a  h  )] 

5  r4  r4  y4  r5  y4  4 


B  +B  coth (a  h  )  Y*M  =  - 
L  4A  4  B  y4  4  JL  5  J 


r4 


c  coth (a  h  ) 
r4  y4  4 

B  - - - l-B 

4  A  CC  4B  a 

y4  y4 


4A 


c  coth (a  h  ) 
tm  r4  y4  4 

Y  +  - - - 

5  a 

y  4 


B  Y‘'coth(a  h  ) +c  /a 
4B  5  y4  4  r4  y4 


TM 


(A. 41) 


(A. 42) 


(A. 43) 


(A. 44) 


(A. 45) 


(A. 46) 


(A. 47) 


(A. 48) 


(A. 49) 
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Using  (A.26d)  and  factoring  e^/a  out  of  the  left  hand  side, 

tm  tm  ,  .  .  2 

a  Y  Y  +(e  /a  ) 

B  =  -  b  5  sA  y* 

4A  4B  C  TM  ,  TM 

r4  Y  +  Y 

5  S  4 


(A. 50) 


Again  this  expression  is  very  similar  to  (A. 27a),  so  let 


tm  tm  ,  ,  ,2 

Y  Y  +(e  /a  ) 
rTM  _  5  s4  r4  y4 

'4  =  yTM  yTM 

5  S4 


(A. 51) 


Then  (A. 50)  simplifies  to; 


a 

TM  y4 


B  =  -  B  Y 
4  A  4B  4  C 


(A. 52) 


r4 

Finally,  the  fields  at  the  center  conductor  interface  are  matched. 
Initially,  the  electric  fields  at  the  interface  are  simply  set  equal  to 
each  other.  The  condition  that  the  fields  are  zero  over  each  of  the 


center  conductors  is  applied  implicitly  through  the  current  or  electric 
field  expansion  functions  when  Galerkins  method  is  used  in  the  final 


step. 


9)  E 


TM 

X3 


TM 

TM 

TM 

=  E 

,  or  E 

—  E 

X 

O 

II 

> 

y=d  z3 

1  3 

y=d  z4 

J  3 

(3  a 

X  v3 

B  - - —  cosh(0)  =  -  B 

3  B  wM3c3 

a  e  u 
y4  r3  r3 

g  =  _  g  J _ 

3B  4B  a  E  Lt 

y3  r4  r4 

Using  this  relation  in  (A. 40) 


/3  a. 
x  y4 

4  B  LOU  E 
4  4 


y=do 


cosh (0) 


(A. 53) 


(A. 54) 


B 


3  A 


a  u 

B  yTM  y4^r3 

4B  3  E  U 

r4^r4 


(A. 55) 


Finally,  the  discontinuity  of  the  magnetic  fields  due  to  the 


current  densities  on  the  center  conductor  interface  is  used. 
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10)  H 


TM 

X3 


-  H™ 
y=d3  x4 


=  j™ 

v=d  z 

3  3 


0  0 

.  2  .  z  TM 

n  —  B  cosh(0)  -  j  —  B  cosh(0)  =  J 
J  U  3A  J  ii  4A  Z 

^3  ^4 

Using  (A. 52)  and  (A. 55) 

a  a  u 

B  Y™  -yl  +B  Y™  -£i  =  -j  -i  J™ 

4  B  3  C  4B4C  0  Z 

r4  r4  1  z 

u  c 

„  ,  _TH  4  r4  1 

B  =  -j  J  - -  - 

4  B  z  0  a 


(A. 56) 


TM  TM 
z  y4  Y  +Y 
3  4 


(A. 57) 


(A. 58) 


To  verify  that  the  TMy  mode  satisfies  the  boundary  conditions  in 
the  spectral  domain,  it  is  neccesary  to  show  that  using  the  z  directed 


magnetic  fields  at  the  interface  gives  the  same  result  as  using  the  x 
directed  fields. 


(TM 
z  3 


-  H™ 
y=d  z  4 

3  3 


=  _j™ 

y=d3  x 


0  0 

-j  —  B  cosh(0)  +  j  —  B  cosh(0)  =  -J™ 
J  U.  3A  J  U  4A  X 

3  ^4 

Using  (A. 52)  and  (A. 55) 


a  a  u 

B  y™  +B  y™  -yi  =  -j  _i  j™ 

4B  3  E  4B  4  C  J  f3  X 

r4  r4  x 

ME 

„  .  TM  ^4  r4  1 

B  =  -j  J  Jr -  - 

4  b  x  0  a 


-  ,,TM ,  TM 
x  y4  Y  +Y 
3  4 


Equation  (A. 61)  is  equivalent  to  (A. 58)  if 


J™/0  =  J™/0 

X  X  z  z 


(A. 59) 


(A.  60) 

(A. 61) 


(A. 62) 


or 


0  J™  =  0  J™ 

Z  X  X  z 


(A.  63) 


which  was  shown  to  be  true  for  the  TMy  case  in  Section  2.3.2.  Thus  the 


TMy  mode  does  satisfy  the  boundary  conditions  in  the  spectral  domain 
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and  can  be  used  independently  to  solve  the  boundary  value  problem.  To 

y 

complete  the  TM  solution,  substitute  (A. 58)  into  (A. 52) 


LL  Y 

.  TM  4  4 

B  =  j  J  -r—  - 

4A  Z  Y™+y™ 

3  4 


(A. 64) 


Using  (A. 58)  and  (A. 64)  for  B  and  B  ,  substituting  back  into 

4B  4A 

the  equations  for  the  electric  field  in  region  4,  equations  (A.8a-b), 


and  setting  y  =  d^,  the  fields  at  the  center  conductor  interface  are; 


J™  0 

TM  .  Z  'x  TM 

E  =  j -  — —  Z 

x4  we  0 
0  z 


(A. 65b) 


TM  .  Z  TM 

E  =  j -  Z 

z4  we 

0 


(A. 65a) 


where 


TM  TM  TM 

’j  =  Y  +Y 
4  3 


(A. 66) 


TM  TM 

Y^  and  Y^  can  be  thought  of  as  the  modal  admittances  seen  by  the 

current  densities  looking  outward  in  each  direction  from  the  center 

■  TM  .  , 

conductor  interface.  Thus  Z  is  the  parallel  combination  of  these 

admittances . 

.  V  Y 

Having  completed  the  solution  for  the  TM  mode,  the  TE  modal 

configuration  is  considered.  The  vector  potentials  for  this  mode, 

F  (0  ,y,0  ),  are  chosen  as  follows  for  each  of  the  six  regions: 
yi  x  z 

REGION  1 :  0  s  y  <  h 


F  =A  (0  ,0  ) sinh (a  y) 
yl  1  x  z  yl 

REGION  2 :  h  <y<d,d  =h+h 
-  1  2  2  12 

F  =A  (0  ,0  ) cosh [a  (d  -y)]+A  (0  ,0  ) sinh [a  (d  -y) ] 

yl  2A  '  x "  z  y2  2  2B  x  '  z  y2  2  J 

REGION  3 :  d  iysd,d  =h+h+h 
-  2  J  3  3  12  3 

Fyl=A3A  Pz>  COSh  [<Xy  3  (d3_y)  1  +A3B(Px'Pz)sinh  [“y3  (d3'y)  1 


(A. 67a) 


(A. 67b) 


(A. 67c) 
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REGION  4 :  d  SySd,d  =h+h+h+h 
-  3  J  4  4  1234 

F  —A  (p  ,P  )  cosh  [a  (y-d  )]+A  (P  ,p  )sinh[a  (y-d  )  ]  (A.67d) 

yl  4A  x  z  y4  3  4  B  x  z  y4  3 

REGION  5:  d  i  y  s  d  ,  d  =  h  +h  +h  +h  +h 
-  4  1  5  5  12345 

F  =A  (P  ,13  )cosh[a  (y-d  )]+A  (p  ,f3  )sinh[a  (y-d  ) ]  (A.67e) 

yl  5A  x  z  y5  4  5B  x  z  y5  4 

REGION  6:  d  ^  y  £  d  ,  d  =  h  +h  +h  +h  +h  +h 
-  5  J  6  6  123456 

F  -AO  ,p  ) sinh [a  (d  -y) ]  (A.67f) 

yl  6  x  z  y6  6 

Where  the  A  (8  ,  (3  )'s  are  unknown  functions  of  the  transform 
N  x  z 

variables,  which  are  to  be  determined  from  the  boundary  conditions. 

They  are  abbreviated  as  A^  throughout.  FyN,s  are  chosen  based  on  1) 

the  tangential  E  fields  vanishing  at  y  =  0  and  2)  the  tangential  E 

fields  vanishing  at  y  =  d  .  Again,  only  the  x  directed  fields  are  used 

6 

,  y 

because  of  the  independence  of  the  TE  mode.  The  fields  in  each  region 
are  determined  using  the  six  vector  potentials  (A.67a-f)  along  with 
equations  (2.38a-f). 


REGION  1: 


T  E  Z 

E  =  -j  A  —  sinh (a  y) 
xl  J  l  c  ylJ 

p  a 

T  E  X  Vl 

H  =  -  A  - —  cosh  (a  y) 

xl  l  yl 

2  2  2  2  2 

p  -  13  +  P  -  a  =  u  n  c 

1  x  z  yl  ^11 


REGION  2: 


P 

T  E  .  Z 

E  =  -j  - 

X2  J  C 

2 


.TE 


[* 

«.  r 


cosh [a  (d  -y)]+A  sinh [a  (d  -y) ] 
2  A  y2  2  1  2B  y2  2  J 


sinh[a  (d  -y)]+A  cosh[a  (d  -y) ] 
2 A  y2  2  3  2B  y2  2  * 


2  2  2  2  2 

p  -  p  +  p  -  a  =  U  fl  E 
2  x  z  y2  ^2  2 


(A. 68a) 

(A. 68b) 
(A. 69) 

(A. 70a) 

(A. 70b) 
(A. 71) 
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REGION  3: 

TE 

E 

X3 

3  1 

hte 

X3 

^3e3L 

^3  = 

2  2 

■  0  +  0 

X  z 

REGION  4: 

ete 

x4 

4  L 

TE 

EZ4 

■ s*  r» 

hte 

x4 

1 3  a 
x  y4 

LOLL  E 

4  4 

2  2 
K  +  K 

REGION  5: 

ete 

x5 

■  ^  r 

hte 

X5 

^5C5 

i  cosh [a  (d  -y)]+A  sinh[a  (d  -y) ] 
3  A  y3  3  1  3B  y3  3  1 


\  sinh[a  (d  -y)]+A  cosh[a  (d  -y)  ] 
3  A  y3  3  3B  y3  3  1 


2 

X 

y3 


3  3 


A  cosh [a  (y-d  ) ] +A  sinh[a  (y-d  )] 
4A  y4  3  4B  y4  3 


4  A 


yA 


A  B 


yA 


A  sinh [a  (y-d  ) ] +A  cosh [a  (y-d  )] 
4  A  y4  J  3  4  B  y4  3 

2  2 

-a  =  u)  u  c 
y  4  4  4 


5  A 


y5 


5B 


ys 


h 


sinh [a  (y-d  ) ] +A  cosh [a  (y-d  )  ] 
k  y5  4  5B  y5  4 


_  2  Jl  2  2  2 

0  =  /3  +  0  -  a  =  to  n  c 

'5  '  X  “ z  y5  ^5  5 


REGION  6: 


,TE 

sinh [a 

=  “  j  A  - 

X  6 

6  e 

6 

ye 

T  E 

0  a 

-A  X  y6 

cosh [a 

X  6 

6  LOU  E 
^6  6 

y6 

2  2  2  2  2 
/3  =  0  +  0  -  a  =  u  V  c 

6  '  x  z  y6  *6  6 


Now  the  boundary  conditions  are  applied  to  eliminate 


constants,  beginning  with  the  continuity  of  the  tangential 
the  first  dielectric  interface  above  the  ground  plane. 


(A. 72a) 

(A. 72b) 
(A. 73) 

(A. 74a) 

(A. 74b) 

(A. 74c) 

(A. 75) 

(A. 76a) 

(A. 76b) 
(A. 77) 

(A. 78a) 

(A. 78b) 

(A. 79) 
the  unknown 
fields  at 
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te  te  te  te 

1)  E  =  E  or  E  =  E 

xl  y=h  x2  y=h  zl  y=h  z2  y=h 

1  1  3  1  3  1  3  1 


z  , 

nA  — sinh (a  h 
lc  yi 

l  3 


h  > .  -i  ^r» 
y1  1  c2L 

Crl  f 

) -  A  cosh(a 

lc  2  A 
r?  I- 


cosh(a  h  ) +A  sinh(a  h  )  (A. 80) 

2  A  y2  2  2B  y2  2  I 


A  =  sech (a  h  ) - A  cosh (a  h  ) +A  sinh(a  ) 

1  yl  1  cr2 L  2 A  y2  2  2B  y2 


(A. 81) 


TE  TE  TE  TE 

2 )  H  =  H  or  H  =  H 

xl  y=h  x2  y=h  zl  y=h  z2  y=h 

3  1  3  1  1  J  1 

3  a 

-  A  — X  y^cosh (a  h  )  = 
l  w/iici  yl  l 


j  a  r 
x  y2 

- —  A  smh(a  h  )  +  A  cosh  (a  h  ) 

2 C 2  |_  2 A  y2  2  2B  y2  2 


(A. 82) 


Substituting  in  for  A^  using  (2.3-81) 


a  coth (a  h  ) 
yl  yl  1 


i  cosh (a  h  )  +  A  sinh(a  h  ) 
2 A  y2  2  2B  y2  2 


A  sinh  (a  h  )  -  A  cosh (a  h  ) 
I  2A  y2  2  2B  y2  2 


a  coth  (a  h  )  a  coth  (a  h  )  a 

yi  yi  i  A  y2  y2  2  +A  _y2 

II  2  A  U  2B  U 

r 1  r2  r2 


a  a  coth  (a  h  ) 

„  /  /  ,2  y2  R  y2  y2  2 

A  (a  /u  )  — A  — - - - 

2A  y2  ^r2  U  2B  li 

^  r2  r2 


a 

a  coth (a  h  ) 

y2 

yi  yi  l 

H 

H  , 

r2 

rl 

a  coth (a  h  )  a  coth  (a  h  )  a 

-a  .  yi  yi  1  y2  v2  2  +  _y£ 

2  a"  u  nu 

r 1  r2  r2 

L  J  J  L  J 

Now  the  continuity  of  the  tangential  fields  at  the  first 
dielectric  interface  below  the  cover  sheet  is  enforced. 


(A. 83) 


(A. 84) 


(A. 85) 


TE  TE  TE  TE 

3)E  =  E  ,  or  E  ,  =  E 

x5  y=d  x 6  y=d  z5  y=d  z6  y=d 

3  5  J  5  ^  5  J  5 


-  n  A  — sinh  (a  h  )  = 
J  6  e  y  6  6 

6  3 


3  r A 


cosh (a  h  )  +  A  sinh (a  h  ) 
5  A  y 5  5  5B  y5  5 


(A. 86) 
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A  =  sech(a  h  ) -  A  cosh  (a  h  ) +A  sinh(a  h  )  (A. 81) 

6  y 6  6  e  |_  5 a  y5  5  5b  y5  5  J 


4)  H 


TE 

X5 


TE 

.  =  H 

y=d5  X6 


or  hTE 
y=d5  25 


=  hte 
y=d5  26 


y=dc 


P  a 

■  x  \r  6 

-A  - — cosh  (a  h  )  = 

6  U>/d  c  y  6  6 


6  6 


P  a  r 

V5  I 

- —  A  sinh (a  h  )  +  A  cosh(a  h  ) 

w/j  e  [  5A  y5  5  SB  y5  5 


(A. 88) 


Substituting  in  for  A  using  (A. 87) 

6 


a  coth (a  h  ) 
y6  y6  6 


r  6 


A 


cosh (a  h  ) +A  sinh (a  h  ) 
5  A  y5  5  5B  y5  5 


]' 


a  r 

y5 

A 


r5 

a  coth (a  h  ) 
y6  y6  6 


sinh (a  h  )  -  A  cosh (a  h  ) 
5  A  y5  5  5  B  y5  5 


r  6 


a  coth (a  h  )  a 

A  - — — - —  +A  Jfi 

5  A  U  5B  11 

rr5  r5 


a 


ys 


5B  /i 


r5 


a  a  coth  (a  h  ) 

,  .  .  2  y5  y5  y5  5 

-  A  (a  /u  )  -  A  — - - - 

5A  y5  'rS  U  5B  u 

~r5  'rS 

a  coth  (a  h  )  a  coth (a  h  ) 
y6  y6  6  y5  y5  5 


r  6 


r5 


5  A 


a  coth (a  h  ) 
y6  y6  6 

r  J 

a  coth (a  h  ) 
y5  y5  5 

a 

ys 

2] 

^rs 

1 

(A. 89) 


(A. 90) 


(A. 91) 


V  .  .  i 

Just  as  in  the  TM  case,  the  notation  is  simplified  by  letting 


TE 

Si 


TE 
S  2 


rTE 

S3 


TE 
S  4 


TE 
S  5 


^T  E 
36 


a  coth (a  h  ) 

yi _ yi  i 

»*rl 

a  coth (a  h  ) 
y 2  y2  2 

u 

r  2 

a  coth (a  h  ) 
y3 _ y3  3 

3 

a  coth (a  h  ) 
y 4  y4  4 

a  coth (a  h  ) 
ys _ ys  5 

^rS 

a  coth (a  h  ) 
y6  y6  6 


(A. 92a) 

(A. 92b) 

(A. 92c) 

(A . 92d) 

(A. 92e) 

(A. 92f ) 
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and 


YTEYTE+(a  /„  )2 

si  S2  y2  ^r2 

yTE  +  yTE 

SI  S2 

YTE  YTE+  (a  /f!  )2 

S6  sS  y5  'rs 

yTE  +  yTE 

S  6  S  5 


(A. 93a) 


(A. 93b) 


Using  these  expressions  in  equations  (A. 85)  and  (A. 91)  and  solving  them 


for  A  and  A  respectively; 


2B 


fi  YTE  YTE+(a  /fi  ) 2  /i 

_  r2  si  s2  y2  r2  r2  te 

_  _  —  _  A  _  Y 

2  A  a  TE  T  E  2A  a  2 

y2  Y  +  Y  y2 


Si 


S  2 


H  Y  Y  +(a  /p  ) 

*r5  s 6  s5  y5  ^r5 

A  =  -  A  -  - - -  =  -  A 

5A  SB  a 

y5 


yTE  +  yTE 
S  6  S  5 


r5  yTE 
SB  a  5 

ys 


(A. 94a) 


(A. 94b) 


The  next  boundary  of  concern  is  the  2-3  dielectric  interface,  and 
so  the  tangential  fields  are  forced  to  be  continuous. 


5)  E 


TE 

X2 


TE 

TE 

TE 

^  =  E 

,  or  E 

=  E 

y=d  x3 

1  2 

y=d  z3 

1  2 

y=d 

1  2 

Z  3 

y=d 


A  cosh (a  h  ) +A  sinh (a  h  ) 
3  A  y3  3  3B  y3  3 


'  z  pz  r 

-j  A  —  =  -j  -  i 

2a  e  e 

2  3  L 

er2r  1 

A  =  -  A  cosh  (a  h  ) +A  sinh  (a  h  ) 

2A  e  [  3A  y3  3  3B  y3  3  J 


6)  H 


TE 

X2 


=  H 

y=d  x3 

T  2 


TE 

,  or  H 
y=d  z  2 

*  2 


=  hte 

y=d  z3  y=d 

2  2 


pa  Pa 
x  y2  _  x  y3 

2B  toil  e  wii  e 

2  2  3  3 


a  a  e 
y2  _  y3  r2 


2B  u  a  c  3A 

r2  ^r3  r3  L 


A 


[A  sinh (a  h  ) +A  cosh (a  h  ) 

3  A  y  3  3  3B  y  3  3  J 

a  h  ) 
y3  3  J 


sinh  (a  h  ) +A  cosh(c 
y3  3  3B 


Substituting  in  for  A^  using  (A. 94a) 


-AY 


a  e 
te  _  y3  r2 


2  A  2  U  C 

r3  r3 


2  A 


A  sinh (a  h  ) +A  cosh (a  h  ) 

L  3  A  y  3  3  3  B  y  3  3  J 


a  e  r 
y3  r2 

-  A  sinh (a  h  ) +A  cosh ( 

u  e  3 a  y3  3  3B 

^r3  r3  J 


(A. 95) 

(A. 96) 


(A. 97) 

(A. 98) 

(A. 99) 

(A. 100) 
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Equating  (A. 96)  and  (A. 100) 

e  r 

r2  I 

-  A  cosh  (a  h  ) +A  smh(c 

e  I  3A  y3  3  3b 


r2  r  .  1 

- A  cosh (a  h  ) +A  sinh (a  h  )  = 

er3 [  3 a  y3  3  3b  y3  3  J 

oc  c  r  q  r  -j ' 

_  y3  r2  a  sinh(a  h  ) +A  cosh  (a  h  )  YTE 
/i  C  I  3 A  y3  3  3 B  y3  3  J  [  2  J 


(A. 101) 


Multiplying  both  sides  by  [Y  c  ]/[e  sinh(a  h  )] 

1  2  r3  r2  y3  3 

[a  coth(a  h  ) +A  1  FyTE"|  =  -  JEirA  +A  coth  (a  h  ) 
[  3A  y3  3  3BJ  [  2  J  Id  l  3A  3B  y3  3 


r  a  coth  (a  h  ) 
yTE+  y3  y3  3 

2  U 

L  rr3 


]-*4 


Y  coth (a  h  ) +a  /u 
2  y3  3  y3  rr3 


(A. 102) 


(A. 103) 


Using  (A.  92c)  and  factoring  ji^/a  out  of  the  left  hand  side, 

li  YTE  YTE+(a  //d  ) 2 

^r3  2  S3  y3  ^r3 

A  =  -  A  -  - - - 

3  B  3  A  a  TE  TE 

y3  Y  +  Y 

2  S3 

Since  the  above  equation  is  very  similar  to  (A. 93a).  let 

YTE  YTE+  (a  /id  )2 
yte  2  S3  y3  ^r3 

3  yTE  +  yTE 

2  S3 

Then  (A. 104)  simplifies  to; 

u 

.  „  T  E  r3 

A  =  -  A  Y  - 

3B  3A  3  a 

y3 


(A. 104) 


(A. 105) 


(A. 106) 


As  before,  the  tangential  fields  at  the  4-5  dielectric  interface 
are  now  equated. 

TE  I  TE!  TE  I  _TE  I 


TE 

TE 

TE 

,  =  E 

or  E 

=  E 

y=d  xA  y=d 

A  A 

Z5  y=d 

A 

z4  y=d 

1  4 

-j  A  -  =  -j  -  A 

5A  C  C  4A 

5  4  L 


cosh (a  h  ) +A  sinh ( 
y  4  4  4  B 


a  h  ) 
y4  4 


A  =  — 
5  a  e 


cosh (a  h  )+A  sinh ( 
*  y4  4  4  B 


a  h  ) 
y  4  4  J 


(A. 107) 


(A. 108) 


TE 

TE 

TE 

^  =  H 

or  H 

=  H 

y=d  x  A  y=d 

A  A 

zs  y=d4 

Z  4 

0  a  £  a  •, 

A  — x  yS  =  -X - A  sinh (a  h  ) +A  cosh(a  h  ) 

5B  4  A  y4  4  4  B  yl  4  J 


(A. 109) 
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a  a  e 

ys 


5B  H 


U  L  r 

y4  r5 
u  c 

r5  r4  r4  L 


A  sinh(a  h  ) +A  cosh (a 
4  A  y4  4  4  B  y4 


v] 


Substituting  in  for  A  using  (A. 94b) 

5  B 


TE 

-  A  Y 

5  A  5 


a  C  r 

a  c  r 
y4  r — [a 


5  A 


sinh (a  h  ) +A  cosh  (a 
4 A  y4  4  4B  y4 


sinh (a  h  ) +A  cosh  ( 
4  A  y  4  4  4  B 


v] 

w][*r] 


-i 


Equating  (A. 108)  and  (A. 112) 


C  r 

—  A 
Cr4  L 


cosh (a  h  )+A  sinh ( 
4  A  y4  4  4B 


a  h  )  = 

y4  4  J 

-  — - — — —  [a  sinh(a  h  ) +A  cosh  (a  h  )  1  |"yT  E1 

e  ,L  4A  y*  4  48  y*  4  J  L  5  J 


-1 


r4  r4 


Multiplying  both  sides  by  [Y  c  ]/[e  sinh (a  h  )] 

5  r4  r5  y4  4 

[a  coth  (a  h  ) +A  1  ["yT E~1  =  -  — —  [a  +A  coth  (a  h  )1 
|_  4A  y 4  4  4bJ  L  5  J  fl  L  4 A  4B  y4  4  J 

[a  coth  (a  h  )n  r 

yTE+  A  1  =  -  A  [YTEcoth(a  h  )+a  / [i  1 

5  Mri,  J  4a[_  5  y4  4  y4  ^r4j 

Using  (A.92d)  and  factoring  u  /a.  out  of  the  left  hand  side, 

r4  y4 

H  YTE  YTE+  (a  //i  ) 2 

^r4  5  s  4  y 4  r4 

ft  =  •  ft  _____  _ 

4  B  4  A  a  TE  ,  TE 

y4  Y  +  Y 
J  5  S  4 

Again  this  expression  is  very  similar  to  (A. 93a) 

YTE  YTE+  (a  /n  )2 
yte  5  s 4  y 4  ^r4 

4  TE  TE 

Y  +  Y 
5  S  4 

Then  (A. 116)  is 


A  =  -  A  j”  ^ 
4  B  4A  4  a 


(A. 110) 


(A. Ill) 

(A. 112) 


(A. 113) 


(A. 114) 

(A. 115) 


(A.  116) 


(A. 117) 


(A. 118) 


y4 

Having  completed  all  the  other  boundary  conditions,  the  fields  at 

,  y 

the  center  conductor  interface  are  matched.  As  with  the  TM  mode,  the 
electric  fields  at  the  interface  are  intitially  set  equal  to  each 
other.  The  condition  that  the  electic  fields  are  zero  over  a  finite 
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distance  in  the  spatial  domain  is  applied  implicitly  through  the 
current  or  electric  field  expansion  fuctions  used  in  Galerkin's  method. 

,TE 


9)  E 


TE 

TE 

TE 

—  E 

,  or  E 

^  =  E 

y=d  x4 

3 

y=d3  z  3 

y=d  z  4 

J  3 

X3 


/3Z  (3z 

-j  —  A  cosh(O)  =  -  j  —  A  cosh(O) 
€  3A  J  C  4  A 

3  4 


y=d 


A  =  A 
3A  4A  E 


r3 


(A. 119) 


(A. 120) 


r4 


Using  this  relation  in  (A. 106) 


u  e 

TE  r3  r3 


A  =  —  A 
3B  4  A  3  a  e 


(A. 121) 


y3  r4 

Finally,  the  discontinuity  of  the  magnetic  fields  due  to  the 
current  densities  on  the  center  conductor  interface  is  used. 


10)  H 


TE 

X3 


y=d 


^TE 
X  4 


TE 

.  —  J 

y=d  z 
3 


pa  Pa 

A  ^2  cosh(0)  +  A  -iLVl  cosh  (0)  -  JTE 
3B  <4!  e  4  B  10 fU  E  Z 


(A. 122) 


3  3 


4  4 


Using  (A. 118)  and  (A. 121) 


u>u  c 

-  A  y”  -  A  YTE  =  -±±  J?E 
4  A  3  4A  4  p  z 


x 


,  Uli  E 

1  04  jTE 

^4A  “  "  yTE+yTE  0x  Z 

3  4 


(A. 123) 


(A. 124) 


y  y 

As  with  the  TM  case,  it  is  neccesary  to  verify  that  the  TE  field 

configuration  satisfies  the  boundary  conditions  independently  in  the 

spectral  domain.  For  this  to  be  true,  using  the  z  directed  magnetic 

fields  at  the  center  conductor  interface  should  yield  the  same  answer 

as  using  the  x  directed  fields. 


188 


x3|y=d3  x4|y=d3  x 

pa  Pa 

A  Z  y3  cosh(O)  +  A  -  Z  y/|  cosh(O) 
3  B  wM3C3  4B 

Using  (A. 118)  and  (A. 121) 


(A. 125) 


TE  TE 

A  Y  —  A  Y 
4A  3  4A  4 


i  wu  e 
1  0  4  TE 

A  —  -  "  ■  - -r -  J 

4  A  YTE+yTE  Pz  X 

3  4 


0)1 1  C 

0  4  TE 

-p7  x 


(A. 126) 


(A. 127) 


Equation  (A. 127)  is  equivalent  to  (A. 124)  if 


-  JTE/p  =  JTE/p 

X  Z  Z  X 


(A. 128) 


p  JTE  =  -  p  JTE 
XX  z  z 


(A. 129) 


•  y  . 

which  was  shown  to  be  true  for  the  TE  case  in  section  2.3.3.  Thus, 

y  y 

both  the  TE  and  TM  modes  satisfy  the  boundary  conditions  in  the 
spectral  domain  independently.  To  complete  the  solution,  substitute 
(A. 124)  into  (A. 118) . 


Y  U  0)11  c 
4  r4  'O  4  TE 

A  —  -  -  - rr -  J 

4  B  TE  TE  a  p  Z 

Y  +Y  y4  1  x 
3  4 


(A. 130) 


Using  (A. 124)  and  (A. 130)  for  A  and  A  ,  and  substituting  back 

4B  4A 

into  the  equations  for  the  electric  field  in  region  4,  equations 
(A.74a-b),  and  setting  y  =  d3,  the  fields  at  the  center  conductor 
interface  are; 


P 

TE  .  TE  TE  Z 

E  =  10)11  J  Z  — 
X4  J  Z  P 

X 


(A. 131a) 


TE  .  TE  TE 

E  =  -  10)11  J  Z 
z4  0  z 


(A. 131b) 
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Looking  over  the  derivation,  it  is  evident  that  there  is  a  pattern 

to  the  determination  of  the  constants.  Observing  the  equations 

TM  TE 

defining  the  Y  ;  (A.27a-b),  (A. 39),  and  (A. 51),  and  Y  ;  (A.93a-b), 

i  i 

(A. 105),  and  (A. 117),  it  is  evident  that  there  is  a  pattern  that  can  be 
used  to  generalize  the  derivation  to  an  arbitrary  number  of  layers. 
Using  this  pattern,  with  the  geometry  in  Fig.  2.2,  the  modal  impedances 
can  be  found  using  the  equations  in  Section  2.3.8. 


